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4 DIMENSIONAL MAGNETIC RESONANCE IMAGING 

This application claims priority to U.S. provisionai application serial No, 
60/226,667, filed on August 22, 2000, the complete disclosure of which is 
5 incorporated herein by reference. 



1. Field of the Invention : 

10 This invention relates to method and apparatus for imaging an object, such as a 

body. More particularly, the mvention provides high resolution, real-lime images in 4 
dimensions with little or no deterioration torn motion artifact. 

15 2. Background of the Invention : 

Nuclear Magnetic Imaging (NMR) which is commonly called magnetic 
resonance imaging (MRI) entails I .) magnctizijig a volume with a constant primary 
magnetic field in a z-direction, 2.) providing u gradient along the axis of the z-directed 

20 field to select a slice in the xy-plane, the plane perpendicular to the direction of the 
primary field, 3.) providing electromagnetic radiation resonant with the Larmor 
frequency of protons in die slice, 4.) providing a puJse of resonant electromagnetic 
radiation to flip the magnetization vector into the transverse plane or plane of the slice, 
and 5.) applying a magnetic field gradient along an axis in the xy-planc of the z-dircctcd 

25 field with excitation at the Lannor frequency to provide phase dispersion of the NMR 
signal along tlie axis to encode spatial information, and 6.) recording the free induction 
decay (FID) radio emission signals following excitation, 7.) recording a plurality of such 
FIDs, each recorded following an excitation with a rotated direction of the gradient in tlie 
xy-planc, and 8.) reconstructing the image firom the plurality of the FIDs. An integer n 

30 of FBDs each having a phase gradient that corresponds to the magnetic field gradient that 
was rotated to n unique directions in the xy-plane comprise a set along two orthogonal 
axes in phase or k-space. A two dimensional Fourier transfonn of the data set is used to 
reconstruct an n by n pixel image. 

MRI is of primary utility in assessing brain anatomy and pathology. But Jong 

35 NMR relaxation times, a parameter based on how rapidly excited nuclei relax, have 

prevented NMR from being of utility as a high resolution body imager. The most severe 
limitation of NMR technology is that for spin echo imaging n, the number of free 
induction decays ("FIDs"), a nuclear radio fxequeacy energy emitting process, must 
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equal the number of Unes in the image. A single FID occurs over ^proximately 0. 1 
seconds. Not coiisidcnag the spin/lattice relaxation time, the time for the nuclei to 
reestablish equilibrium following an RF pulse, which may be seconds, requires an 
irreducible imaging time of n times 0. 1 seconds, which for 5 1 2 x 5 1 2 resolution rejquires 
5 approximately one minute per each two dimensional slice. Hiis represents a multiple of 
1 500 times longer that the time that would freeze organ movements and avoid image 
deterioration by motion artifoct. For example, to avoid deterioration of cardiac images, 
the imaging time must not exceed 30 msec. A method for speeding NMR imaging flips 
the magnetization vector of the nuclei by less than 90 degrees onto the xy-plane, and 
10 records less FIDs. Such a method, known as the flash method, can obtain a 128 x 128 
resolution in approximately 40 seconds. Another techm'quc used lo decrease imaging 
time is to use a field gradient and dynamic phase dispersion, corresponding to rotation of 
the lield gradient, during a single FID to produce imaging times typically of 50 msec. 
Both methods produce a decreased sigiial-to-noise ratio ("SNR") relative to spin echo 
1 5 methods. The magnitude of the magnetization vector which links Uie coil is less for the 
flash case because the vector is flipped only a few degrees into the xy-plane. The echo- 
planar technique requires shorter recording times with a concomitant increase in 
bandwidth and noise. Both methods compensate for decreased SNR by increasing the 
voxel size with a concomitant decrease in image quality. Physical limitations of these 
20 techniques render obtaining high resolution, liigh contrast vascular images impractical. 

SUMN4ARY OF THE INVENTION 

It is an object of the invention to provide high resolution multi-dimensional 
25 images of an object, such as a body, tissue, or working cardiopuhnonary system. 

It is a fwther object of the invention to rapidly acquire the data to provide 
magnetic resonance images of a body with reduced motion artifact<5- 

These and other objects of the invention are attained by providing an apparatus 
for obtaining a magnetic resonance image ofa body using data acquired over the tluec 
30 spatial dimensions plus time, rather than acquiring data only in time at one receiving 
antenna. A preferred embodiment of the apparatus of the invention includes a radiation 
source for applying a first radiation field having a magnetic component to tlie body, to 
magnetize the body. 'J*he apparatus further includes a source for applying a second 
radiation fidd to the body, to elicit a radiation field fh>m the body. A detector senses 
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this radiation field, and produces a signal that a reconstruction processor employs to 
create tlic magnetic resonance image of Ihe body. 

A NMR image is obtained of a magnetized body from a tJuee-dimensional map 
of the intensity variation of the NMR signal produced by each voxel of the magnetized 
5 body, and detected over a three-dimensional volume of space external to the body, herein 
referred to as the "sample space " The data is acquired over three spatial directions plus 
time. In an embodiraent, the NMR signals are detected over a three dimensional detector 
array as a function of time. The NMR signals may be sampled at least at the Nyquist 
rate, i.e., at a rate tliat is twice the highest temporally frequency of the NMR signal and 

10 twice the highest spatial frequency the Fourier tiansforai of the NMR image of the 
phantom. Sampling at the Nyquist rate or higher allows the spatial variations of the 
external NMR signal to be acquired. 

hi an embodiment, the NMR signal at each detector as a function of time is 
processed by a method such as a Fourier transform operation to give a plurahty of 

1 5 Fourier components each having the same frequency, an intensity and a phase angle. 
The NMR signal of each voxel at any given detector gives rise to a Fourier component 
with a unique phase angle relative to tlie Fourier component of any other voxel of the 
phantom at that detector, llic s^ of Fourier components that correspond to the NMR 
signal of a ^ven voxel over the detectors is determined. This may be achieved by using 

20 a first component having a phase angle and calculating the phase angle as a function of 
• spatial position of tlie first detector relative to any other detector and identifying the 
component at each detector having the calculated phase angle. The sets are detennined 
for all of the voxels. A NMR image is obtained from the sets wherein each set of Fourier 
components comprises a three-dimensional map of the intensity variation of the NMR 

25 signal produced by each voxel of the magnetized body. This practice of the invention 
preferably employs a Fourier Iransfomi algorithm, described in Fourier Transform 
Reconshuction Algorithm Section, to determine the spatial location of each voxel from 
the corresponding set of Fourier components comprises a three-dimensional map of the 
intensity variation of the NMR signal produced by each voxel of the magnetized body 

30 over the sample space. Iliis is repeated for each set to fonn the NMR image of the 
object. 

NMR images produced according to prior art methods and systems rely on 
applying an additional magnetic field in die direction of tlic primary field having a 
gradient along an axis in the transverse plane to cause a phase variation of the NMR 
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signal along the axis in the transverse plane. The direction axis of the gradient is varied 
a plurality of tinies to gives rise to an equivalent number of lines in the reconstructed 
image. In the present invention, the unique phase variation of the NMR signal is 
provided by the combination of 1 .) the angle 9 suspended between the direction of the 
5 detcx^tor and the radial vector, the vector from the dipolc to the detector, and 2.) the angle 
^ due to a separation distance r between a voxel and a detector given by the 
wavenumber of the RF field k times r . 

My prior inventions disclosed in U.S. Patent No. 5,073,858 and US. Patent 
Application No. 09/191,454, the complete disclosures of both which are incorporated 

1 0 herein by reference, are in part based on the realization that matter having a permeability 
different from that of free space distorts a magnetic flux applied thereto. This property is 
called magnetic susceptibility. An object, herein called a phantom, can be considered as 
a collection of small volume elements, herein referred to as voxels. When a magnetic 
field is applied to the phantom, each voxel generates a secondary magnetic field at the 

15 position of the voxel as well as external to the phantom. The strength of the secondary 
magnetic field varies according to the strength of the applied field, the magnetic 
susceptibility of the material within the voxel, and tlie distance of the exterual location 
relative to the voxel. For example, my U.S. Patent No. 5,073,858 and U.S. application 
No. 09/191,454 (cache that the net magnetic flux at a point exhinsic to a phantom to 

20 which a magnetic field is applied, is a sum of the applied field and the external 
contributions from each of the voxels. The '858 patent fiirther teaches sampling the 
external flux point by point and anploying a reconstruction algoritlun, to obtain the 
magnetic susceptibility of each voxel from the sampled external flux. 

Unlike the '858 patent that relies on a static response from a magnetized body to 

25 detemiinc the magnetic susceptibility of the body, my invention disclosed in U.S, Patent 
Application No. 09/191,454 elicits a radiative response from a magnetized body by 
subjecting the body to a resonant radiation field. One embodiment of the '454 invention 
generates a three-dimensional magnetic susceptibiUty imageof an object including a 
patient placed in a magnetic field from a three-dimensional map of a radio frequency 

30 (RF) magnetic field external to the patient, induced by subjecting selected nuclei of the 
body to a resonant RP field. Application of an RF pulse to the body causes the body to 
emit the RF magnetic flux external to the body. A Fourier transform of this external flux 
produces its frequency components ("Larmor fitxjucncies"). Each Lannor frequency is 
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used to dctemime the magnetic susceptibility of the voxels of the body producing that 
Lannor frequency. Further, the intensity variation of the external RF field over a three- 
dimensional volume of space is used to determine the coordinate location of each voxel. 
One practice of the inventions disclosed in my U.S. Patent No. 5,073,858 and 
5 U.S. Patent Application No. 09/19! ,454 obtains a thrc^^mensional magnetic 

susceptibility map of a magnetized body from a three-dimensional map of a secondary 
magnetic flux produced by the magnetized body, and detected over a three-dimensional 
volume of space external to the body, herein referred to as the "sample space." The 
extrinsic magnetic flux is samipled at least at the Nyquist rate, i.e., at twice the spatial 

10 frequency of the highest frequency of the Fourier transform of the magnetic 

susceptibility map of tlie phantom, to allow adequate sampling of spatial variations of the 
external magnetic flux. This practice of the inventions preferably employs a Fourier 
traiTsfonn algorithm, described in Fourier Transform Reconstruction Algorithm Section, 
to form the magnetic susccplibility map of the object. 

15 The present invention relates to systems for providing images of distributions of a 

quantity, in a chosen region of tlie body, by gyromagnetic resonance, particularly nuclear 
magnetic resonance (NMR) techniques. Such techniques may be used for examining 
bodies of different kinds. A particularly beneficial application is the cxanrinalion of 
patients for medical piuposes. Unlike my US. Patent No. 5,073,858 and U.S. Patent 

20 Application No. 09/1 91,454, the present invention employs nuclear magnetic resonance 
(NMR) to induce a magnetized phantom of essentially constant magnetic susceptibility 
to emit an cxtemal radiation having a mag^ietic field component In particular, 
^plication of an RF piilse, resonant with selected nuclei of a magnetized body, can 
polarize the nuclei through rotation of their magnetic moments. Tlie polarized nuclei 

25 within a voxel precess about the local magnetic field Ln the voxel at a Lannor frequency 
determined by the applied magnetic fieJd at position of the voxel. The superposition of 
external RF fields produced by all the voxels of the body creates the total external RF 
field at each detector that is time dependent. The external RF field recorded at tlie 
detectors as a function of time contains components each having a unique phase angle 

30 relative to other components. Each component corresponds to an emitting voxel of the 
phantom. The time dqjcndent signal at each detector may be transformed into a series of 
coinponeuts having intensity and phase data. Each set of components of the >MR signal 
over the sample space due to a given voxel is detennined from the phase data and the 
detector positions. The spatial variation of the NMR signal over tlie sample space is 
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usod to dctennine the location of the voxel in the phantom. Tliis is repeated for all sets 
of components, each corresponding to a voxel to reconstruct the NMR image. 

The radiation source for magnetizing a body to be imaged can be a direct current 
("DC") magnet, including a superconducting magnet. The radiation sources and 
5 amp]} fiers for applying an RF pulse to the magnetized body are well known in the art, 
and include, but arc not limited to, klystrons, backward wave oscillators, Gunn diodes. 
Traveling Wave Tube amplifiers, A preferred embodiment of the invention employs a 
three dim^siona! array of antennas as detectors for sensing the external RF field. 

One practice of the invention dctex:l5 the external RF field in the near field region 
1 0 where the distance of a detector sensing radiation from a voxel at a distance r from the 
detector is much smaller than the wavelength X of the radiation emitted by the voxel, 
i.e., r «X (or kr « 1 ). The near fields are quasi-stationary, that is they oscillate 
harmonically as e"'*" y but are otherwise static in character, Hius, the transverse RF 
magnetic field of each voxel is tiiat of a dipole. In one embodiment, an array of 
1 5 miniature RF antennas sample the external RF field over a three-dimensional volume of 
space that can be either above or below the object to be imaged. The distance r of the 
detector from the voxel gives rise to the phase term c"*' of the component of the 
detected RF signal where A is the wavcnumber of the NMR signal The harmonic 
oscillation of each RF dipole is equivalent to the dipole rotating in the transverse plane. 
20 The detector is responsive to a component in this plane. At a point in time, each 

dipole IS directed at an angle 0 relative to tlie direction of detection of the detector. The 
phase angle 0 of tlie RF dipole relative to the direction of detection axis of the detector 
gives rise lo a phase angle term e''° . In a preferred embodiment, the sum of the phase 
angles, kr and 0 , are unique for each voxel at each detector. The position of eadi 
25 detector relative to a different detector may be used to calculated the phase angle of the 
second relative (o tlae first. This may be repeated over all of the detectors to give the set 
of intensities of tfie NMR signal over the sample space due to a voxel. The location of 
each voxel is detennined through the spatial variations of tiic intensity of the NMR field 
of the set of components associated by phase. Thus, the phase of the components of the 
30 external RF radiation, and the intensity variations of the external RF radiation provide 
the necessary information for providing a NMR image of the magnetized phantom, such 
as a human body. Such a NMR image can be employed to obtain anatomical images of a 
human body based on selected physiological parameters. 
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Iji an embodiincnt of the present invention, Uic NMR image of an object 
including a patient placed in a magnetic field is generated from a three-dimensional map 
of the transverse resonant radio frequency (RF) magnetic fliix external to the patient. 
The external RF field recorded at the detectors as a function of time contains components 
5 each having a unique phase angle relative to other components. Each component 
corresponds to an emitting voxel of the phantom. The time dependent signal at each 
detector may be transformed into a series of components having intensity and phase data. 
Bach set of components of the NIVIR signal over the sample space due to a given voxel is 
detennined from the phase data and the detector positions. The intensity variation of the 
1 0 transverse RF field over the sample space is used to delennine the coordinate location of 
each voxel. The RF field is the near field wliich is a dipole that serves as a basis clement 
to form a unique reconstruction. The geometric system function corresponding to a 
dipole which determines the spatial intensity variations of the RF field is a band-pass for 
K'^K' Preferably, each volume clement is reconstructed independently in parallel with 
1 5 all other volume elements such that the scan time is no greater tlian tJie nuclear free 
induction decay (FID) time. 



Secondary Magnetic Field 

20 The magnetic moment wi, of each voxel is a magnetic dipole. And the phantom 

can be considered to be a tlireo-dimcnsional array of magnetic dipoles. At any point 
extrinsic to the phantom, the z-componcnt of the secondary flux, b\ from any single 
voxel is 

25 where x, y, and z are the distances from the center of the voxel to the sampling point. Jt 
is shown in APPENDICES I-IV that no geometric distribution of raagoeUc dipoles can 
give rise to Eq. (1). Therefore, the flux of each magnetic dipole (voxel contribution) 
forms a basis set for the flux of the array of dipoies which comprise the NMR image of 
the phantom. 

30 Eq, ( I) is a system function which gives the magnetic flux output in response to a 

magnetic dipole input at the origin. The phantom is an array of spatially advanced and 
delayed dipoles weighted according to the magnetic moment of each voxel; tliis is the 
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input function. The secondary flux is the sfupeiposition of spatially advanced and 
delayed flux, according to Eq. (1); this is the output function. Thus, the response of 
space to a magnetized phantom is given by the convolution of Eq. (1) with tiie series of 
weighted, spatially advanced and delayed dipoles representing the NMR image of the 
5 phantom. 

In Fourier space, the output function is the product of the Fourier transform (FT) 
of the system function and the FT of the input function. Thus, tlie system function filters 
the input fuiiction. The output function is the flux over all space. However, virtually all 
of the spectrum (information needed to reconstruct the NMR image) of the phantom 
1 0 exists in the space outside of the phantom because the system function is essentially a 
band-pass filter. This can be appreciated by considering the FT, H[k^yk^], o(Eq. (I): 

»,* 

where is the spatial frequency in the xy-planc or -plane and is the spatial 
frequency along the z-axis. H[k^fkJ is a constant for k^ and k^ essentially equal as 
1 5 demonstrated graphically in FIGURE 1 c. 



Band-Pass Filter 

When a static magneHc field Ho with lines in the direction of the z-axis is applied 
20 to an object comprising a materia] containing nuclei such as protons ihai possess 
magnetic moments, the field magnetizes the material. As a result a secondary field 
superposes the appL'ed field as shown in FIGURE 9. In the applied magnetic field, the 
magnetic moments of each nuclei precesses about the applied magnetic field. However, 
the magnetization of any one nucleus is not observed from the macroscopic, sample. 
25 Rather the vector sum of the dipole moments fix)m all magnetic nuclei in the sample is 
observed This bulk magnetization is denoted by the vector M . In thermal equilibrium 
with the primary field Bq, the bulk magnetization M is parallel to . The 
magnetization vector then comprises magnetic dipole m . The secondary field outside of 
the object (phantom) and detected at a detector 301 is tliat of a series of magnetic dipoles 
30 centered on volume elements 302 of the magnetized material. In Cartesian coordinates, 
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the secondary magnetic flux, B\ at the point (jf j',2) due to a magnetic dipole h;iving a 
magnetic dipote moment at the position (x^^yf^y z©) is 



D 



. (22' -x'-/) 

where i, is the unit vector along the z-axis. The field is the convolution of the system 
function, h(x,y^z) orh(p,^,z) (the left-handed part of Eq. (4)), with the delta fimction 
(the right-hand part of Eq. (4)), at the position (x^jy^^z^). A veiy important theorem of 
Fourier analysis states that tlic Fourier transform of a convolution is the product of the 
1 0 individual Fourier transforms [2]. The Fourier transform of the system function, 
/j(r,y,z) orh(p^(p,z), is given in APPENDIX V. 

The z-coniponent of a magnetic dipole oriented in tlie z-direction has the system 
function^ h(x,y,z)^ which has the Fourier transfomi, ^K^^y^K]* w^^^^^ is shown in 
FIGURE Ic. 

4;rfL'l-A:/l 

V 

The output function, the secondary magnetic field, is the convolution of the 
system function, ^(x,y,2) -lhe geometric transfer function for the z-component of a z- 
oriented magnetic dipole with the mput fimetion-a periodic array of delta functions each 
20 at the position of a magnetic dipole corresponding to a magnetized volume element. 

[x -^y +z j 

The Fourier transform of a periodic array of delta functions (the right-hand side of Eq. 
(7)) is also a periodic array of delta functions in k-space; 
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By (he Fourier Theorem, the Fourier transform of the spatial output function, Eq. (7),' is 
the product of tlic Fourier transform of the system fiuiction given by Eq. (6), and the 
Fourier transform of the input function given by Bq. (8). 



(9) 

\ y, z./ 

In the special case that 



the Fourier transform of the system function (the left-hand side of Eq. (9)) is giveu by 

H=4n {)!) 
Thus, tlie Fourier transfonn of tlie system runctioii band-passes the Fourier trans fonn of 

10 the input fiinction. Both the input function (the right-hand pait of Eq. (7)) and its Fourier 
transfonn (the right-hand part of Eq. (9)) are a periodic array of delta functions. No 
frequencies of the Fourier transform of tho input fimction arc attenuated; thus, no 
infonnation is lost in the case where Eq, (10) holds. 

In an embodiment of the present invention, the magnetization vector is rotated 

1 5 into the transverse plane by an additional RF field H, , The magnetization vector then 
comprises a rotating magnetic dipolc m in the transverse plane. The NMR image may 
be reconstnicfed by sampling the external field from a series of RP dipolcs rather lhan 
thai from a series of static dipolcs. in this case, the Fourier transform of the system 
function also band-passes the Fourier transform of the input function. Thus, tlie 

20 resolution of the reconstructed NMR image is limited by the spatial sampling rate of the 
secondary RF magnetic field according to the Nyquist Sampling Theorem. 

Reconstruction 

25 The NMR image may be reconstituted using a Fourier transfonn algorithm. The 

algorithm is based on a closed-form solution of tlie inverse problem — solving the spatial 
distribution of an array of magnetic dipolcs from the measured extrinsic secondary (RF) 
field that is transverse to the magnetic flux that magnetizes the voxels. The Uansvcrsc 
Rl^" magnetic field of each voxel is that of a dipole, tlie maximum amphtudc is given by 

30 Hq. (I) wherein the Lannor frequency of each voxel is essentially the same, and , the 
magnetic moment along the z-axis» of Eq. (1) corre^q^onds to the bulk maguctization M 
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of each voxel. In terms of the coordinates of Eq. (I), an array of miniature RF antennas 
point samples the maximum dipoie compon^t of the RP signal over the sample space 
such as the half space above (below) the object to be imaged wherein each RF signal is 
has a unique phase siiift due to the relative spatial relationship of each voxel and each 
5 detector. The phase shi ft of each component and the relative q>at!al relationship of the 
detectors is used to assign a component from each detector to a set. Each set of 
components associated by phase comprises the spatial variations of the intensity of the 
transverse RF field of any given voxel, Tiie intensity variation over the sample space is 
used to determine the coordinate location of each voxel. In the limit, each volume 
10 element is reconstructed independently in parallel with all other volume elements such 
that the scan time is no greater than the nuclear free induction decay (FID) time. 

The NMR scan performed on the object to be imaged including a human 
comprises the following steps: 

15 • The magnetic moments of nuclei including protons of tlie object to be imaged that are 
aligned by the primary field are further aligned by a radio frequency (RF) pulse or 
series of pulses. 

• The strength and dination of the rotating H, (RF) field that is resonant with the prototis 
20 of (he magnetized volume and is oriented perpendicularly to tlie direction of the 

magneti'/ang field is applied such that the final precession angle of the magnetization 
is 90 (^„^ = 90° ) such tliat the Rl^" dipolc is transverse to the primary magnetizing 
field and perpendicular to the RI^ magnetic field detector. 

25 • NMR pulse sequences which piovide the signals for a ?; or image may be applied. 
For example a 90 pulse may be followed by a series of 1 80 pulses. One sequence 
is the Carr-Purcell-Meiboom-Gill (CPMG) sequence [3]. 

• The fi"ee induction decay signals are recorded. 

30 

• The time dependent signals are Fourier transformed to give the intensity and phase of 

each component. The NMR signal of each voxel at any given detector gives rise to a 
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Fourier component with a unique phase angle relative to the Fourier component of 
any other voxel of the phantom at that detector. 

• The matrix of Fourier components that correspond to tlie NMR signal of a given voxel 
5 over the detectors is dctemiincd. ITiis may be achieved by using a first component 

having a phase angle and calculating the phase angle as a function of spatial position 
of the first detector relative to any other detector and identifying the component at 
each detector having the calculated phase angle. The matrices are determined for all 
of the voxels. The measurements of the spatial variations of the transverse RF field 
10 of a given matrix is used to determine the coordinate location of each voxel. Thus, 

each matrix of components associated by phase comprises the intensity variation over 
the sample space of the RF field of tlie bulk magnetization tVl of each voxel. 

• The Fourier transfonn algorithm is performed on each set of components over the 
15 detector array to map each bulk magnetization M corresponding to a voxel to a 

spatial location over the image space. The bulk magnetization map (NMR image; 
.also the input function) is given as follows. With respect to the coordinate system of 
Eq. (1), (x, y, and z are the Cartesian coordinates, , the magnetic moment along 
the z-axis, of Eq. (1) corresponds to the bulk magnetization M of each voxel, and b' 

20 is the magnetic flux due to the magnetic moment shown in FIGURE 9; the 

relationship to the NMR coordinate system is given in the Reconstruction Algorithm 
Section) the origin of the coordinate system, (0,0,0), is tlic center of the upper edge of 
the phantom. The phantom occupies tlie space below the plane x, y, e=0{?<Oin 
the phantom space)> and the sampling points lie abovcthcplane (z> 0 in the 

25 sampling space). The magnetic flux iii the sampling space is given by multiplying 

the convolution of the input function with the system function by the unitary fimction 
(one forz > 0 and zero for z < 0). The input function can be solved in closed-form 
via the foilovviug operations: 

30 U Record the RF NMR signal at discrete points in tlie sampling space. Bach point is 
designated (x, y, z, RF) and each RF value is an clement in matrix A. The time 
depeiKleut signals are Fourier transformed to give the intensity and phase of each 
component. The NMR signal of each voxel at any given detector gives rise to a Fourier 
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component with a unique phase angle relative to the Fourier component of any otlier 

voxel of the phantom at that detector. The matrix of Fourier components that correspond v 
to theNMR signal of a given voxel over the detectors is detennined. This may be 
achieved by usbg a first component having a phase angle and calculating the phase angle 
5 as a function of spatial position of the first detector relative to any other detector and 
identifying the component at each detector having the calculated phase angle. The 
matrices are determined for all of the voxels. The measurements of tlie spatial 
variations of the transverse RF field of a given matrix is used to determine the coordmatc 
location of each voxel. Thus, each matrix of components associated by phase comprises 
10 the intensity variation over the sample space of the RF field of the bulk magnetization M 
of each voxel. 

2. Discrete Fourier transform each matrix A„ to obtain each matrix . 

1 5 3. Muhiply each clement of each matrix by tlie corresponding inverse (reciprocal) 
value of the Fourier transform of the system function, Eq. (1), evaluated at the same 
frequency as the element of the matrix . This is matrix C,. 

4. Generate matrix by taking the discrete inverse Fourier transform of matrix C„. 

20 

5. Multiply each clement of each matrix by tlie distance squared along the z-axis to 
which the element corresponds to generate the position of tfie bulk magnetization M of 
voxel n. (This corrects the limitation of the sample space to z> 0). 

25 6. In one embodiment, the point spread of Oie reconstructed voxel is corrected by 

assigning one voxel above a certain threshold \vilh Uic bulk magnetization M . The other 
voxels are assigned a zero value. 

7. This procedure is repeated for all matrices A^, In the limit with suiBcient phase 
30 resolution, each volume element is reconstructed iiidependently in parallel wiUi all oilier 
volume elements such tjiat the scan time is no greater tha;: the nuclear free induction 
decay (FID) time. 
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. BRIEF DESCRIPTION OF THE DRAWINGS 



FIGURE la is a schematic diagram of an apparatus according to the teachings of 
the invention in which a magnet magnetizes a body to be imaged, an RF generator 
5 excites the body, and an array of detectors detect radiation emitted by the body in 
response to the RF excitation radiation; 

FIGURE 1 b is a schematic diagram of an apparatus according to the teacliings of 
the invention that employs a three-dimensional array of detectors to detect radiation 
emitted by nuclei of a body to be imaged in response to excitation of the nuclei by a 
1 0 radi ation source; 

FIGURE Ic is a plot of tlie Fourier Uransfonn H[k,,ky,k^] of the system function 

h(x,y,2) (Eq. (1)) corresponding to the z-componcnt of a magnetic dipole oriented in the 
z-dircction in accordance with the invention; 

FIGURE 2 is tl^e plot of the field of a ring of dipoles of radius R and magnetic 
15 moment m = 10"* Gcm^ given by Eq. (114) as a function of radius R where the position 
of the center of (he ring relative to the detector is the point (0,0, 1^) accordance wiOi 
the invention; 

FIGURE 3 is the plot of the field of a ring of dipoles of radius R = 0.2 cm and 
magnetic moment m-W* Gcm^ given by Eq. (1.14) as a function of the distance 
20 between the detector at the origin and the center of the ring at the points (0,0, 
z^Acmtoz-\San)\i\ accordance with the invention; 

FIGURE 4 is the plot of the field of a shell of dipoles of radius R and magnetic 
moment m = 10^ Gctn given by Eq. (II. 1 7) as a fimction of radius R where the position 
of the centOT of the shell relative to the detector is the point (0,0,10) in accordance with 
25 the invention; 

FIGURE 5 is Uie plot of the field of a shell of dipoles of radius R = 0.2 cm and 
magnetic moment w = 10^ Gan given by Eq. (A. 1 7) as a function of the distance 
between the detector at the origin and the center of the shell at the points (0,0, 
z= 4 cm to 2 = \ S cm) in accordance with the invention; 
30 FIGURE 6 is the plot of the field of a sphere of dipoles of radius R and magnetic 

moment /?i = 10* Gcm^ given by Eq, (IV. 16) as a function of radius R where the 
position of the center of the sphere relative to the detector is tlie point (0,0,10) in 
accordance widi the invention; 
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FfGURE 7 is the plot of the field of a sphere of dipolcs of radius R = 0.2 cm and 
magnetic moment m = 10^ Gcm^ given by Eq. (TV.! 6) as a function of the distance 
between the detector at the origin and the cent^ of the sphere at the points (0,0, 
z= 4 cmtoz = i5 €m)in accordance with the invention; 
5 FIGURE 8 shows a typical the nuclear magnetic resonance (4D-MRI) apparatus 

in accordance with the invention; 

FIGURE 9 shows the coordinate system (x, y, z) of Eq. (I) with a primary field 
Hg and the corresponding magnetic dipoic both oriented parallel to the z-axis wherein 
the zrcomponent of the flux due to the z-oriented dipole is measured at a detector 
10 according to Eq. (1 ) and shows the distances from the voxel to the detector in accordance 
with tlic invention; 

FIGURE 10 shows the stationary coordinate system of the nuclear magnetic 
resonance (4D-MRI) appai atus of FIGURE 8 corresponding to the coordinate system of 
FIGURE 9 in accordance with the invention; 
1 5 FIGURE 1 1 shows the rotating NMR coordinate system ( x ^ , y ^ , ) and the 

stationary coordinate system (x, y, z) of the NMR detector corresponding to the 
coordinate system of FIGURE 9 and FIGURE 10 of a primary field oriented parallel 
to the -axis and the z-axis and the corresponding transverse RF magnetic dipole 
oriented in the x^yj^ -plane and periodically parallel to the y-axis wherein the spatial 
variation of the RF y-component of the flux due to the RF dipole is measured at a 
detector according to Eq. (1) and shows the distances from the voxel to the detector in 
accordance witii the invention; 

FIGURE 12 is a schematic of the three dimensional detector (antennae) array 
witli respect to the stationary NMR coordinate syst^ of FIGURE 10 which coTresponds 
to the coordinate systems shown in FIGURES 9 and 1 1 in accordance with the invention; 
FfGURE 13 shows the general process of reconstruction by reiteration, and 
FIGURE 14 shows the stationary coordinate system (x, y, z) of the NMR detector 
corresponding to the coordinate system of FIGURE 9 and FIGURE 1 0 of a primary field 

oriented parallel to tlie z,, -axis and the z-axis and the corresponding transverse RF 
magnetic dipole oriented in tiie Xny^ -plane and periodically parallel to the y-axis 
wherein the spatial variation of the RF y-component of the flux duo to the RJF dipole is 
measured at a detector according to Eq, (I) and shows the distances and angles between a 
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voxel linear to a first detector, a second nonaligned voxel, and a second detector in 
accordance witb the invention; 

Further details regarding specific derivations, and calculations are provided in the 
attached s^pendices, wherein: 
5 APPENDIX I is the field produced by a ring of dipoles according to the present 

invention; 

APPENDIX n is the derivation of the field produced by a shell of dipoles 
according to (he present invention; 

APPENDIX in is the niathematical proof that the field produced by a shell of 
1 0 magnetic dipoles is different &om that of a single dipole according to the present 
invention; 

APPENDIX IV is the derivation of the field produced by a sphere of dipoles 
according to the present invention; 

APPENDIX V is the derivation of the Fourier transform of the system ftmction of 
15 the z-component of the magnetic field from a dipole oriented in the direction of the z- 
axis used in the reconstiiiction process according to the present invention; 

APPENDIX VI is the derivation of S' - //F® U(k^) convolution from Eq, (55) 
used in a reconstruction process according to the present invention, and 

APPENDIX VU is the derivation of tlic Inverse Transform of Eq. (69) to Give 
20 Inverse Transform I , Eq. (69), used in a reconstruction process according to the present 
invention. 

DETAILED DESCRIPTION OF THE INVENTION 

25 An exemplaiy embodiment of a nuclear magnetic resonance (4D-MR1) apparatus 

1 10 according to the teachings of the present invention is shown in FIGURE la. The 
apparatus 1 1 0 includes a magnet 1 12, such as a superconducting magnet, that provides a 
primary or magnetizing field, to magnetize a body 114 to be imaged, A magnetometer 
1 16 detemiines the primary or magnetizing field in tlie volume to be occupied by the 

30 body, i.e., in the image space, in the absence of the body. One practice of the invention 
utilizes a magnetometer that employs NTvlR of protons in water for determining the 
primary radiation field at multiple points in the image space. In an embodiment, tlie 
primary field is uniform as recorded by magnetometer 1 16, 
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A radiation source 1 18, such as a radio firequency generator, applies an RF pulse 
to the body in combination with RF antennae 1 18a when the body is placed in the 
magnetic field, to excite and thereby polarize selected nuclei of tlie body. The excited 
nuclei emit an RF radiation that a plurality of detectors 1 20» disposed in a plane above or 
5 below the object, can detect. The excitation pulse can be selected to rotate the 

magnetization of the nuclei, preferably by 90 degrees, with the respect to the primaiy 
field. In such a case, the RP radiation that the excited nuclei emit is primarily along a 
direction perpeudicuJar to tlie plane of the detectors, in an embodiment using techniques 
known by those sb'lled in the art, RF pulse sequences are applied to generate the data for 
10 Tj and/or 2\ NMR images. The detectors can be selected to respond only to components 
of a magnetic field perpendicular to the plane in which they reside. Thus, such detectors 
can detect the emitted RF field without interference from other components of the 
magnetic field to permit a unique reconstruction of the NMR image. In an embodiment, 
a drive mechanism 1 22 moves the detectors 120 in a direction perp^dicular to the plane 
1 5 of tlic detectors, to sample the external RF field over a three-dimensional volume. The 
separation of the detectors and the step size of the movement of the detectors along a 
direction perpendicular to the plane are selected such that the detectors sample the 
external RF field over a three-dimensional volume, i.e., the sample volume, at least at the 
Nyquist rate. Preferably, a tliree dimension array of detectors is used to sample the RF 

20 field at the Nyquisl rale. Such a sampling advantageously allows obtaining the NMR 
image of the body. An embodiment, employs an impedance-matched array of RF 
auteiuias that are time-multiplexed to reduce cross talk among them. The RF field may 
be sampled synchronously or the field may be sampled at known times so that llie phase 
at any given detector may be related to that at any other detector, 

25 The process of inagnetization of selected nuclei of the body can be better 

understood by referring to FIGURE 1 1 tliat shows a magnetizing field applied to a 
voxel 14a ofthe body 114. The RF excitation field is selected to be in a direction 
perpcndiculai' to the field Hp, and its magnitude is designated as H,. NMR active nuclei 
ofthe voxel, such as protons, possess botli angular momentum and a magnetic moment. 

30 Thus, the vector sum ofthe magnetic moments of all such NMR active nuclei present in 
tlK3 voxel I4a give rise to ihe bulk magnetization ofthe voxel I4a. The bulk 
magnetization vector M (not shown) designates the collective contribution of a selected 
type of NMR active nuclei to the magnetization ofthe voxel 1 4a which corresponds to a 
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RF parallel magnetic moment m . The RF excitation field is selected to be in resonance 
witli tlie selected type of nuclei, e.g., protons, of the voxel 14a, to rotate the bulk 
magnetization vector M . Jn a rotating fiame representation, designated in FIGURE 1 1 
as (Xr , , Zj^) and well known to those skilled in NMR, the magnetization vector M 
5 rotates about 11 , so long as the RF radiation is present The rate of rotation of the 
* magnetization vector about the applied field H, depends on the gyromagnetic ratio of 
the alTcctcd nuclei and the magnitude of the excitation field, H, . The duration of the RF 
radiation can thus be selected to cause a rotation of the magnetization vector, iniliaily 
aligned along H^,, onto the x^y,! -plane, for example by a 90*" rotation. After the RF 
10 excitation field is turned off, the rotated magnetization M processes about H,, at the 
position of tlie voxel 14a. 

The precession of the magnetization M of each voxel about H,, produces a 
radiating dipolar field corresponding to a magnetic moment m 'that has its maximum 
intensity along the y ^ direction. The precession fi-equcaicy of the magnetic moment 
1 5 called the Lamior fi-cquency is dependent on the magnetic flux of the primary field 

which is uniform in a preferred cmbodimrat A detector 20a can be, for example, an RF 
antenna pointing along the y direction to respond selectively to the radiating RF field 
from all of the voxels. This external RF field is recorded as a function of lime at each 
over tlic entire sample space. The NMR signal at each detector as a function of time is 
20 processed by a method such as a Fourier tiansfonn operation to give a plurality of 
Fourier components each having an iiuensity and a phase angle. Tlie NMR signal of 
each voxel at any given detector gives rise to a Fourier coinponent with a unique phase 
angle relative to tlie Fourier component of any other voxel of the phantom at that 
detector. Each set of components of the NMR signal over the sample space due lo a 
25 given voxel is d^ermined iiom the phase data and the detector positions. The spatial 
variation of the NMR signal over the sample space is used to determine the location of 
the voxel in tJie phantom. This is repeated for all sets of components^ each 
corresponding to a voxel to reconstruct the NMR image. Preferably, the position of each 
voxel is reconstructed independently in parallel with all other voxels such that the scan 
30 time is no greater the time for the excited nuclei to return to their uncxcited state called 
the nuclear fi^ee induction decay (FID) time. This may be achieved by using sufficiently 
higli spatial and temporal sampling rate such that each set of components associated by 
phase corresponds to a single voxel. 
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With further reference to FIGURE la, an analog to digital converter ("A/D 
converter") 124 convctls the analog signal outputs of the detectors 120 into digital 
signals. A preferred embodiment of the invention utili2C5 at least a 1 2 bit A/D converter 
to digitize the output signals. A reconstniction processor 1 26 receives the digital signals, 
5 and determines the NMR image of the body. The reconstniction processor 126 includes 
a Fourier transform processor 126a that obtains the phase components of the external RF 
field. In addition^ a spatial locator 126b, which is a part of the reconstruction processor 
1 26, employs the variation of the maximum intensity of the external RF field over the 
sample space (the tJiree dimensional space sampled by the detectors) to locate the 
10 positions of the voxels in the image space producing a particular set of components of 
the external RF field associated by phase, in a manner described in detail below (See the 
Reconstruction Algorithm Section). An display processor 128 displays a bvo- 
dimensional or a three-dimensional image corresponding to the NMR image of the body. 
The spatial locator 126b employs an algorithm, described in detail below in the 
1 5 Reconstruction AJgorithm Section, to determine the positions of the voxels of the 

magnetized body that produce each set of components of the external RF field associated 
by phase. It is demonstrated in the Uniqueness of Reccmstruction Section that a RF 
magnetic field produced by a geometric distribution of dipoles is unique. Therefore, a 
unique spatial distribution of magnetic dipoles, such as those corresponding to the buljc 
20 magnetization of each voxel due to prccessing nuclei, gives rise to a unique magnetic 
field. Thus, the measured extemal RF field can provide a unique solution for the spatial 
distribution of magnetic dipoles, i.e., magnetized voxels comprising excited nuclei, iii the 
body. 

Referring again to FIGURE 1 a, Uie spatial locator preferably employs a Fourier 
25 Transform Reconstruction Algorithm, described in detail in the Reconstruction 

Algorithm Section, to map a bulk magnetic moment M corresponding to a particular set 
of components associated by phase of the external RF field onto a spatial location or 
locations in the body. The latter case applies if more than one location in the body gives 
rise to a particular phase component of the extemal RF radiation. 
30 As shown in FIGURE J a, a preferred embodiment of the invention employs an 

open design magnet, such as a Helmholtz coil design, to allow positioning the array of 
detectors 120 close to the body 1 14 to be imaged. The orientation of the magnetic field 
with respect to the body can be selected to optimize the signal to noise ratio of the 
si|^als delected by the airay of detectors 120. For example, in case of imaging a patient 
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body, tJie primary magnetic field can be selected to be coaxial with the body, or it caii 
alternatively be perpendicular to the body axis. 

llic 4D-MR1 apparatus of FIGURE la provides a number of advantages. In 
particular, because all of the data is acquired at once rather tlian over hundreds of 
5 repeated sequences of pulses as is the case v/iih IMR systems of the prior art, the 
apparatus of tie invention is particularly suited for imaging cardiopulmonary and 
vascular systems. Further, since the apparatus may increase the number of detectors to 
increase the resolution, the apparatus of the invention may acfiieve a higher resolution 
such aslO^ cm^ ; thereby, permitting physicians to view human anatomy and pathology 
10 in a nianoer not available with conventional imaging techniques. Further, the present 
technique may provide a three duncnsional image that can be displayed from any 
perspective. 

An alternative embodiment of the nuclear magnetic resonance (4D-MR1) 
apparatus of the invention, shown in FIGURE lb, employs a three-dimensional array of 
1 5 detectors 230, spaced apart to detect spatial variations of the emitted RF radiation at least 
at the Nyquist frequency. A magnet 2 12 provides a magnetizing in a volume to be 
occupied by a body, i.e., image space. A magnetometer 216 measures the magnetizing 
field at a plurality of positions in the image space in the absence of the body to determine 
the uniformity of the primary field. As in embodiment of FIGURE la, the body 214 to 
20 be ijnaged is placed in a magnetizing field provided by the magnet 212. An RF 

generator 21 8 in combination with an RF antenna 218a apply an RF pulse or a sequence 
of RF pulses to the body to polarize selected nuclei of the body. The three-dimensional 
array of detectors 230 provide output signals in response to the RF radiation emitted by 
the body. The RF signal may be recorded synchroQOusly to pemiit the relative phases of 
25 the Fovnier components comprising the RF signal In an embodiment, the detection may 
be synchronized relative to the excitation. A digitizer 224 digitizes tlie output signals 
and sends the digital signals to a constmction processor 226 that determines variations of 
the NMR image of the body in a manner similar to that described in connection with the 
embodiment of FIGURE la. A display processor 228 receives the information regarding 
30 the spatial variations of Ihe intensity of the NMR signal from the construction processor 
226, and provides a two-dimensional or a three-dimensional NMR image of tlie body. 
Employing a three-dimensional airay of detectors advantageously decreases the 
acquisition time because the emitted RF signal over the ©itire three dimensional sample 
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space is detected at once. The shortening of (be acquisition time in turn reduces motion 
artifacts in the NMR image. 

Uniqueness of Reconstruction 

5 

The nature of the RF field can be determined from Maxwell's equations applied 
to a sinusoidal current. With a sinusoidal current J(x') confined to small region 
compared vWth a wavelength, the solution of (he vector potential A(3f) is [4] 

10 where Jt = is the wavenumber, and a sinusoidal time dependence is understood. The 
c 

magnetic induction is given by 

B = VxA 03) 
while, outside the source, the electric field is 

k 

1 5 For a source of dimension d , the iiclds in the far zone defined by d «X «r arc 

transverse to the radius vector and fall off as r'\ typical of radiation fields. For the near 
zone where r « A (or kr « I), Qie exponential in Eq. (12) can be replaced by unity. 
Then the vector potential is given by 

20 This shows that the near fields are quasi-stationary, oscillating harmonically as e''^' , but 
otherwise static in character. 

Nuclear magnetic resonance (NMR), which is commonly called magnetic 
resonance imaging (MRI), is a means to measure the primary and secondary magnetic 
fields to provide the input to the NMR reconstruction aigoritlim. The proton 
25 gyromagjietic ratio / 2;r is 

1 2% =42,57602 MHz T"' (16) 
The NMR frequency / is the product of tlie proton gyromagnctic ratio given by Eq. (16) 
and the magnetic flux B . 

/ = / 27cB 42.57602 MHz T'B (17) 
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A typical flux for a superconducting NMR imaging magnet is 0.25 T. According to Eq. 
(17) this corresponds to a radio frequency (RF) of 10.6 MHz whicb corresponds to a 
wavelength of 283 m . In the present invention, each RF antennae of an array is located 
at a distance of about 10 cm jrom the voxels within the image space. Thus, the RF field 
5 is detected in the near zone where r « A (or kr «\\ and the near fields according to 
Eq. (1 5) are quasi-stationary, oscillating harmonically as e'^^ , but otherwise static in 
character. The transverse RF magnetic field of each voxel is that of a RF dipole, the 
maximum amplitude is given by Eq. (1) wherein the Laroior frequency of each voxel is 
determine by the uniform primary field H^, and , the magnetic moment along the z- 

10 axis, of Eq. (1) corresponds to the bulk magnetization M of each voxel. 

An object containing nuclei with a magnetic moment, herein called a phantom, 
can be considered as a collection of small volume elements or voxels. When a static 
magnetic field Hq with lines in the direction of the z-axis is applied to an object 
comprising a material containing nuclei such as protons that possess magnetic moments, 

1 5 the field magnetizes the material. As a result a secondary field superposes the applied 
field as shown in FIGURE 9. In tiie applied magnetic field, tl)e magnetic moments of 
each nuclei precesses about the applied magnetic field. However, the magnetization of 
any one nucleus is not observed from the macroscopic sample. Rather tlie vector sum of 
the dipole moments from all magnetic nuclei in the sample is observed. ITiis bulk 

20 magnetization is denoted by the vector M . In thermal equilibrium with the primary field 
II, , the bulk magnetization M is parallel to IIq - In an embodiment of the present 
invention, the magnetization vector is rotated into the transverse plane by an additional 
RF field II, . The magnetization vector then comprises a rotating magnetic dipole m in 
the transverse plane. The NMR image may be reconstructed by sampling the external 

25 field from a series of RF dipoles. 

T he field streng^ of a magnetic dipole moment is a fimction of the ext^al 
position in space relative to the dipole. For convenience of analysis, the field of a scries 
of static dipoles m having the coordinates shown in FIGURE 9 is analyzed for 
uniqueness. (The uniqueness of the field of a set of static dipoles applies for the 

30 equivalent case of RF dipoles oriented in the transverse plane.) Considering FIGURE 9, 
the net magnetic field at a point extrinsic to the phantom is a sum of tlie applied field and 
the contributions of each of the voxels within the object, the secondary field. The field is 
point sampled over a tliree dimensional space. 
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The secondary magnetic field due to magnetized tissue has to be modeled as 
noninteracting dipoles aligned with Uie imposed field. It is demonstrated below that the 
field of any geometric distribution of dipoles is unique, and the superposition principle 
holds for magnetic fields; therefore, a uruque spatial distribution of dipoles gives rise to a 
5 unique secondary magnetic field, and it is fiirther demonstrated below that this secondary 
field can be used to solve for the NMR image in closed form. It follows that this map is 
a unique solution. To prove that any geometric distribution of dipoles has a unique field, 
it must be demonstrated that the field produced by a dipole can serve as a mathematical 
basis for any distribution of dipoles. This is equivalent to proving that no geometric 
10 distribution of dipoles can produce a field which is identical to the field of a dipole. 

By symmetry considerations, only three distributions of uniform dipoles need to 
be considered. A magnetic dipole has a field that is cylindrically s^onmetricaL A ring, a 
shell, a cylinder, and a sphere of dipoles are the only cases which have this symmetry. A 
cylinder is a linear combination of rings. Thus, the uniqueness of the dipole field is 
15 demonstrated by showing that it is different fi^om that of a ring, a shell, and a sphere. 
The uniqueness of the dipole from the cases of a ring, a shell, and a sphere of dipoles is 
demonstrated in APPENDIX I, APPENDIX H, and Appendix IV, respectively. The plot 
of the three cases of the field of a ring, shell, and a sphere of dipoles each of radius R 
and magnetic moment = 10* Gcm^ given by Eq. 0. 14), Eq. (II. 17), and Eq. 0V.16) of 
20 APPENDIX 1, APPENDK n, and APPENDIX iV as a function of radius R where the 
position of the center of each distribution relative to the detector is the point (0,0,10) is 
given in FIGURES 2, 4, and 6, respectively. Since the fiolds vary as a function of radius 
R , llie dipole field is not equivalent to these distributions of dipoles. Ft is further 
mathematically proven in APPENDIX HI that the field produced by a shell of magnetic 
25 dipoles is different from that of a single dipole. All other fields are a linear combination 
of dipoles. Thus, the dipole is a basis element for the reconstruction of a NMR image. 
Since each dipole to be mapped gives rise to a unique field and since the total field at a 
detector is the superposition of the individual unique dipole fields, linear independence is 
assured; therefore, the NMR image is unique. In other words, there is only one solution 
30 of the NMR image for a given set of detector values which spatially measure the 

superposition of the unique fields of the dipoles. This map can be reconstructed using 
the algorithms described in the Reconstruction Algoritlmi Section. 

Eq. (1) is a system fimction which gives the magnetic flux output in response to a 
magnetic dipole input at the origin. The phantom is an array of spatially advanced and 
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delayed dipoles weighted according to tlic bulk magnetization of each voxel; tliis is the 
input function. The secondary flux is the superposition of spatially advanced and 
delayed flux, according to Eq. (I); this is the output function. Thus, the response of 
space to a magnetized phantom is given by the convolution of Eq. (1) widi the serieis of 
5 weighted, spatially advanced and delayed dipoles representing the bulk magnetization 
map or NMR image of the phantom. The discrete si^aJs are recorded by a detector 
airay over the sample space comprising the xy-planc and the positive z-axis of FIGURE 
9. In an embodiment of the present invention, the magnetization vector is rotated into 
the transverse plane by an additional RP field H, . The magnetization vector then 
10 comprises a rotating magnetic dipole ra in the transverse plane. The NMR image may 
be reconstructed by sampling the external field from a series of RF dipoles. The discrete 
signals are recorded by a detector array over the sample space comprising the xz-p!ane 
and tlie positive y-axis of FIGURE 1 1 . 

15 r, and 7; NMR Images 

The NNfR active nuclei including protons posses both angular momentum and a 
magnetic moment. When nuclei are placed in a static magnetic field Hp, they process 
about the Held at a frequency proportional to tlie magnitude of 11^, The bulk 

20 magnetization M of each voxel comprises the vector sura of the magnetic moments from 
all of the nuclei in each voxel. If the processing nuclei are then subjected to an additional 
rotating (RF) field H,, which is synclwonous with the precession, tlieir magnetic 
moments and thus M will prccess about II ^ and rotate away from the primary field 
by an angle in a coordinate frame which rotates at the Larmor frequency [3]. ITie 

25 precession about H, continues as long as H, exists. Tlie final value of ^^^^ then depends 
on the strength of 11,, which dctennines the precession rate, and the time for wliich it is 
turned on. The nuclei absorb energy as they change their orientation. This is known as 
nuclear magnetic resonance (NMR). The temperature of the nuclei or nuclear spin 
system rises dming absorption of energy. When the H, field is removed, the spin 

30 system cools down until it is thermal equilibrium with its environment, llic exponential 
relaxation of the spin system temperature to that of the surrounding lattice is called span- 
lattice relaxation and has a time constant 7\ where a time constant is defined as the time 
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it takes for 63% of Uie relaxation to occur, llie NMR signal may also decay because the 
nuclei initially in phase following the U, pulse get out of alignment with each other or 
dephasc by local interactions with the magnetic fields of neighbor nuclei. The dcphasing 
of the NMR signal is due to differing precession rates effected by the local interactions 
5 and is described by an exponential time constant also known as spin/spin relaxation. 
The main source of NMR image ( also called magnetic resonance images (MRl)) 
contrast is r, and which depend on tissue types. 

In an embodiment of the present invention, a image is produced by a applying 
at least one pulse sequence that inverts the magnetization and records the relaxation, a 
10 technique called inversion recovery. For example, the RF receivers ate switched on to 
follow the decay following the nuclear excitation comprising a H, pulse. M, (0, the 
time dependent bulk magnetization in the direction of the primary field H^^ (coordinates 
of FIGURE 1 1) is examined at a time after an inverting pulse by applying another H, 
pulse equivalent to a rotation by 90° after waiting the time following the initial 
15 inversioa The 90-^ pulse puis the z magnetization mXO into the transverse plane for 
observation. Changing tfie waiting time allows for observation of MXO at different 
times during relaxation. 

In an embodiment of the present invention, a Tj image is produced by a appfying 
at least one pulse sequence that flips the magnetization vector into the transverse plane 
20 and records the transverse relaxation by producing at least one spin-echo. The 7; image 
depends on the NMR signal decaying because the nuclei initially in phase following the 
pulse get out of alignment with each other or dephase by local interactions with the 
magnetic fields of neiglibor nuclei. The dcphasing of the NMR signal is due to differing 
precession rates effected by the local interactions ajid is described by an ex|x>nential time 
25 constant T ^ . Another unwanted source of dephasing is due to inhomogeneities in the 
primary field across the image space that causes an addition contribution to the 
precession rates of the magnetic nuclei. The signal ones observes after flipping the 
magnetization into the transverse plane (hat includes dephasing from an inhomogeneous 
H, field is still knovm as Uie free induction decay QFIDl but the is called relaxation. 
30 The T , relaxation may be recovered fi-om a T* FID since Uie inliomogeneity is 
constant and may be reversed. The dephasing due to the static inhomogeneity of 11^ 
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may be canceled out by applying a 180° H , pulse at time along the y r -axis as shown 
in FIGURE 1 1 where the relaxation ocxurs in time t^. After an additional Hme 
the total time elapsed after the PO'' H, pulse is li^ which is the spin-echo tinie. At 
this time tlie dephasing due to the static inhomogeneity of is exactly canceled out; so» 
5 the relaxation is strictly due to those processes that create 2\ relaxation whicli is 
recorded. Ajfter each 180^ pulse, another spin echo is formed. The envelop of the 
maximum ampHtude of tlie spin echoes is the relaxation. In an embodiment, a pulse 
sequetice to give the data for a image known as the Carr>Purcell-Mciboora-Gill 
(CPMG) sequence (3] comprises applying a 90° pulse along the x^-axis followed by a 
10 series of 1 80'' pulses along the -axis at times + 2ni^ where n is an integer including 
zero. 

In other words, when nuclei are placed in a static magnetic field and then 
subjected to an additional rotating (RF) H,, which is synchronous witb tbeir 
precession, M will precess about H, and rotate away from the primary field by an 
15 angle . The magnitude of M is a maximum initially and decays with time. This 

occurs by emission of the same multipolarily radiation that it absorbed and by transfer of 
energy to the surrounding lattice. Thz intensity of the radiation is a function of M and 
the coordinate posinon relative lo the RF emitting voxel. In the present invention, the 
measurement of the intensity of the RF signal is perfomied over time and space 
20 following and/or encoding pulses. The signal as a function of time at a given 
detector position is Fourier transformed to give the components each having an 
amplitude and a unique phase. Each set of componenls of the NMR signal over tlie 
sample space due to a given voxel is determined from the phase data and the detector 
positions. The location of each voxel is determined through the spatial variations of the 
25 intensity of the transverse NMR field of the set of components associated by phase. 



4D-MRI System 



An NMR apparatus used to generate, and measure tl\e secondary field and 
30 reconstruct the image is shown in FIGURE 8, and the corrcsponduig coordinate system 
is shown in FIGURE 10. The apparatus comprises 1.) a magnet including a 
superconducting magnet lo magnetize a volume of an object or tissue to be imaged, 2.) a 
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means (magnetometer) to determine the primary or magnetizing flux over the image 
space in tlie absence of the object to be imaged. In one embodiment, NMI< of a proton 
containing homogeneous phantom such as water is determined on a point by point basis 
to map the primary field, or the magnetic flux at multiple points is obtained 
5 simultaneously using the reconstruction algorithm described herein, 3.) a radio frequency 
(RF) generator and transmitter including an antennae such as saddle coils to excite the 
protons of the magnetized volume^ 4,) a means including an antennae coil to sample the 
dipoie component (z-componcnt in terms of Eq. (I)) of the RF secondary magnetic field 
at tfie Nyquist rate in time over the proton free induction decay, 5.) a detector array of 

10 elements of the means to sample the time signals including an antennae array which is 
selectively responsive to the dipoie component (z-componcnt of Eq. ( 1 )) of die RF 
magnetic field of the magnetic moments of the protons which are aligned along the 
transverse axis to ideally point sample tfie secondary magnetic field at the Nyquisi rate 
over the spatial dimensions which uniquely determine the NMR image which is 

1 5 reconstructed from the measurements, 6.) an analog to digital converter to digitize the RF 
signals, 7.) a time Fourier transfoim processor to convert the signal at each detector over 
time into its Fourier components, 8.) a processor to associate the Fourier components due 
to each voxel by phase, 9.) a reconstruction algorithm processor including a Fourier 
transform processor to convert eadi set of con^onents into a voxel location of the bulk 
20 magnetization in the image space in parallel over all tlie voxels to form an NMR image, 
and 10.) image processors and a display such that the NMR image can be rotated in 
space to be displayed from any perspective as a three dimensional or two dimensional 
(tomographic) image. 

25 Magnetizing Field 

In an embodiment of the present invention, the applied magnetizing field which 
permeates the object to be imaged including tissue is confined only to that region which 
is to be imaged. The confined field limits the source of signal only to the volume of 
30 interest; thus, the volume to be reconstructed is limited to the magnetized volume which 
sets a limit to the computation required, and eliminates end effects of signal originating 
outside of tlie edges of die detector array. In the NMR case, the field having a steep 
gradient at the edges limits the imaged region by providing a range of Larmor 
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frequencies wherein the data is comprises a narrow frcqxicncy band at a desired Larraor 
frequency. 



Detector Array 

5 

An embodiment of the NMR imager of (he present invention comprises a detector 
array of multiple detector elements which are arranged in a plane. The array may be a 
two dimensional detector array which is translated over the third dimension during the 
scan, or it may be a three dimensional detector array. The individual detector of the 
1 0 array may respond to a single component of the secondary magnetic field which is 

produced by the magnetized object including tissue where tJie component of the field to 
which the detector is responsive determines the geometric system jfimction which is used 
in the reconstruction algorithm discussed in the Reconstruction Algoritlim Section. The 
detectors ideally point sample the secondary magnetic field at the Nyquist rate over the 
1 5 spatial dimensions which umquely determine the NMR image which is reconslnictod 
from the measurements. 

Small antennas may measure the RF signals as point samples without significant 
decrease in the signal to noise ratio relative to large antennas by using impedance 
matching wliile minimizing resistive losses by using superconducting reactance 
20 elements, for example. In an embodiment, cross talk between anteiuias is ameliorated or 
eliminated by time multiplexing the signal detection over Uie anuy of anteimas. The RF 
field may be sampled synchronously or the field may be sampled at known times so that 
the phase at any given detector may be felaleU to that at any other detector 

Micromaguctic field sensors that are used to detect the primary field in the 
25 absence of the object to be scanned include NMR detectors and superconducting 

quantum interference devices (SQUIDS). Additional devices have been developed that 
are based on galvanometric effects due to the Lorentz force on charge carriers. In 
specific device configurations and operating conditions, the various gal vanomagnetic 
effects (Hall voltage, Lorentz deflection, magnetoresistive, and magnctoconcentration) 
30 emerge. Semiconductor magnetic field sensors include MAGFETs, magnetotransitors 
(Ml), Van der ?auw devices, integrated bulk Hall devices including the vertical MT 
(VMT), and the lateral MT (IMT), silicon oa sapphire (SOS) and CMOS 
magnetodiodes, the magnetounijunction transistor (MUJT), and tfie carrier domain 
magnetometer (CDM). magnetic avalanche transistors (MAT), optoelectronic magnetic 
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field sensors, and magnetorcsistive magnetic field sensors. In llie case of NMR 
measurement of Ihe secondary field (and/or the primary field), the detector array 
comprises RF antennas described in the NMR Primary Magnet, Gradient Magnets, RF 
Generator, RF Transmitter, and RF Receiver Section. 

5 

Scanning Methods 

The NMR scan performed on the object to be imaged including a human 
comprises the following steps: 

10 

• The magnetic moments of nuclei including protons of the object to be imaged that arc 

aligned by the primary field are further aligned by a radio frequency (RF) pulse or 
series of pulses. 

1 5 • The strength and duration of the rotating H, (RF) field that is resonant with the protons 
of the magnetized volume and is oriented perpendicularly to the direction of the 
magnetizing field is applied such that the final precession angle of the magnetization 
is 90 (^„, = 90"* ) such tliat the RF dipole is transverse to the primary magnetizing 
field and perpendicular to the RF magnetic field detector. 

20 

• NMR pulse sequences which provide the signals for a 7J or 7^ image may be applied. 

For example a 90 pulse may be followed by a series of 1 80 pulses. One sequence 
is the Carr-Purceil-Mciboom-GHI (CPMG) sequence [3]. 

25 • The firee induction decay signals are recorded. 

• The time dependent signals are Fourier transformed to give tlie intensity and phase of 

each component. The NMR signal of each voxel at any given detector gives rise to a 
Fourier component with a unique phase angle relative to the Fourier component of 
30 any other voxel of the phantom at that detector. 

• The matrix of Fourier components that correspond to the NMR signal of a given voxel 

over the detectors is detennined. This may be achieved by using a first component 
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having a phase angle and calculating the phase angle as a fimction of spatial position 
of the first detector relative to any other detector and identifying the component at 
each detector having the caicuJated phase angle. The matrices are determined for all 
of the voxels. The measurements of the spatial variations of the transverse RF field 
5 of a given matrix is used to determine the coordinate location of each voxel Thus, 

each matrix of components associated by phase comprises the intensity variation over 
the sample space of the RF field of the bulk magnetization M of each voxel. 

- The Fourier transfomi algorithm given in the Fourier Transform Rcconstraction 
1 0 Algorithm Section is performed on each set of components over the detector array to 
map each bulk magnetization M corresponding to a voxel to a spatiaj location over 
the image space. 

• In one embodiment, the point spread of the reconstructed voxel is corrected by 
1 5 assigning one voxel above a certain threshold with the bulk magnetization M . The 
other voxels are assigned a zero value. In this case» the reconstruction is digital 
versus analog. ]n terms of the signal to noise ratio, the superiority of digital ov^ 
analog is generally known to tl)Ose sidlied in the art of signal processing. 

20 • In the limit vifiih sufficient phase resolution, each volume element is reconstructed 

independently in parallel with all other volume elements such that the scan time is no 
greater than the nuclear free induction decay (FID) lime. 

NMR Primary Magnet, Gradient Magnets. RF Generator, RF Transmitter, and RF 

25 

Receiver 

The primary magnet is that typical of a magnetic resonance imaging apparatus. 
The primary magnetizmg field is in the 2-dircction as shown iji FIGURE 10, usually 
30 coaxial with the patient's body. However, in an another embodiment the primary 
magnetizing field is in the direction which is perpendicular to the patient's body. An 
open primary maguet design is preferable such as a Hchnholtz coil design to 
accommodate tlie three dimensional RF antennae an^ay. A further field is applied to have 
a gradient in the z-direction (FIGURE 10). ITiis provides a unique field in a chosen 
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volume of the patient including a cross-sectional slice of the patient. In another 
embodiment, the NMR (Lannor) frequency from each voxel is determined by the 
magnitude of the magnetizing flux at the voxel position. The primary magnetic field has 
a shaqi gradient at the edges of tlic image space wherein the Lannor frequencies outside 
5 of a selected range are rejected as arising from voxels outside of the image volume. 

A rotating RF field H, , of frequency chosen to cause resonance in the slice is 
then applied. Thus, only Ihc molecules in the slice resonate. The resonance signal from 
the shce can then be detected. NMR pulse sequences which provide the signals for a 7; 
or T; image may be applied. For example a 90 pulse may be followed by a series of 
10 1 80 pulses. One sequence is the Carr-Purcell-Mciboom-Gill (CfUG) sequence [3]. 

An embodiment of the NMR imager of the present invention comprises a RF 
generator 1 4 and 18 and RF transmiHer coils 9 shown in FIGURE 8 which provide the 
rotating U, (RF) field that is resonant with the protons of the magnetized volume and is 
oriented perpendicularly to the direction of the magnetizing field as shown in FIGURE 
15 11. Wlien the precessing nuclei are subjected to the additional rotating (RF) field H, . 
which is s>'nclironous with the precession, iheir magnetic moments and thus M 
precesses about H, and rotate away from the primary field by an angle in a 
coordinate frame which rotates at the Lannor frequency [3]. The precession about H, 
continues as long as H , exists. The final value of then depends on the slxiiRgfh of 
20 fl , , wliich determines the precession rate, and the time for which it is turned on. In an 
embodiment, tlie strength and duration of H, is such that = 90'' such that the dipole 
is oriented in die x^y^-planeof FIGURE I L A rotating coordinate frane is 
traditionally used to explain the physics of NMR [3]. ITius, is terms of the traditional 
NMR coordinate designation as described by Patz [3] and Hounsficld (7J sliown in 
25 FIGURE 1 1 , the X R -axis and y ^ -axis rotate about the primary field H,, oriented parallel 
to the z„ -axis al the Lamior frequency relative to the staUonary NMR coordinate system 
(x. y, z) shown in FIGURE 1 0. Thus, the additional rotating (RF) field H, and the 
transverse RF magnetic field are stationary in the rotating NMR coordinate system, but 
rotate at the Larmor frequency in the stationary NMR coordinate system. Both the 
30 rotating and stationary coordinate systems are shown in FIGURE 11. In an" embodiment, 
the rotating H, (RF) field is along the x„ -axis and rotates the magnetization vector by 
the angle =^90^ to into the x^y^ -plane. In terms of the reconstruction by the present 
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invention, the system functioa of the field corresponding lo the rotated magnetization 
vector is equivalent to that of the detection of the z-component of a z-onented dipole. 
llie y-axis is the unique axis of the NMR system shown in FIGURE 8 and the stationary 
NMR coordinates shown in FIGURBS 1 0 and 11. This axis corresponds to tlie unique 
5 axis of the Fourier Transform Reconstruction Algorithm Section, the z-axis shown in 
FIGURE 9f and the magnetization axis is the z-axis. Tlie magnetization axis of the NMR 
rotating coordinates [3,7] is the Zj^-axis shown in FIGURE 11, The measured transverse 
RF magnetic dipole oriented in the xy-plane is periodically parallel to the y-axis and 
rotates about the primary field oriented parallel to the z-axis at the Larraor frequency, 
10 The RF transmitter includes saddle coils. The RF receiver is a tlu*ee dimensional 

array, or two arrays ISO*' from each other, or two orthogonal pairs of arrays, wherein 
each member of a pair is 180'' from each other. The detector ideally point samples tlie 
RF field at the Nyquisl rate described in the Reconstiiiction Algorithm Section. In one 
embodiment shown in FIGURE 1 1, each detector 20a is a coil antennae perpendicular to 
15 the y-axis and selective to the y-componeut of the RF secondary ma^etic field arising 
from each voxel 14a magnetized in tlie2-axis direction. This two dimensional array is 
translated along ttie y-axis during a scan where readings of the secondary magnetic field 
are obtained as a fimction of the translation. In another embodiment shown in FIGURE 
12, the array 401 is three dimensional comprising multiple parallel two dimensional 
20 arrays 402 wherein each two dimensional array has a plurality of antennae coils 403 
which delect the RF field from each magnetized voxel 404. 

In another embodiment, the rotiting H, (JRF) field shown in FIGURE 1 1 is along 
the y R -axis and rotates the ma^etization vector by the angle ^j,, i= 90* to project into 
the X ^yf^ -plane. In teiins of the reconstruction by the present invention, the system 
25 function of the field corresponding to the rotated magnetization vector is equivalent to 
that of the detection of the z-componeni of a z-oriented dipole. The x-axis is the um'que 
axis of the NMR system shown in FIGURE 8 and the stationary NMR coordinates 
shown in FIGURES 10 and 11. This axis corresponds to the unique axis of the Fourier 
Transfoini ReconshiJction Algorithm Section, the z-axis shown in FIGURE 9, and the 
30 magnetization axis is the z-axis. The magnetization axis of the NMR rotating 

coordinates [3,7] is the z^ -axis shown in FIGURE 1 1 . The detector array receives 
similar signals as those produced by the transmitter, and both are perpendicular to tiie x- 
axis. Each detector is a coil antennae perpendicular to the x-axis and selective to the 
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component of the RF secondary magnetic field ansing frona voxels magnetized in the z- 
axis direction. This two dimensionai array is translated along the x-axis during a scan 
where readings of the secondary magnetic field arc obtained as a function of the 
translation. In another embodiment^ the array is three dimensional comprising multiple 
5 parallel two dimensional arrays 402 as shown in FIGURE 12 with the exception that the 
array 401 is perpaidicular to the x-axis. The measured transverse RF magnetic dipole 
oriented in the xy-plane is periodically parallel to the x-axis and rotates about the 
primary field oriented parallel to the z-axis at the Lanmor fi^equency. 

A linear combination of the cases of rotation of the RF magnetization along the 
10 XR-axis and the y ^ -axis is v^ithin the scope of the present invention and may be adopted 
in a manner straightforward to those skilled in theNMR art to apply the fields and 
detectors described for this invention. 

The method described herein may be performed on a suitable NMR examining 
apparatus such as that shown in simplified form in FIGURE 8, and the corresponding 
1 5 coordinate system is shown in FIGURE 10. Illustrated schematically only arc coils 6, 
which provide Bq , the steady primary field 7» which provide G^, the field gradient in tlic 
X-axis direction as shown, 8 which provide the , the field gradient in the y-axis 
direction as shown, 9 which provide tlie Rl^ field, and 10, which provide G,, the field 
gradient in the z-axis direction as shown. The coils are driven by B^, G^, G^, RF, and 
20 G, drive amplifiers 1 1 ,12, 13 J 4, and 15 respectively, controlled by B^^G^, RF, and 
(7, control circuits 16, 17, 18, and 19, respectively. These circuits caj^ take suitable 
forms which will be well known to those with experience of NMR equipment and other 
apparatus using coil induced magnetic fields. The circuits are controlled by a central 
processing and control unit 20 to achieve the desired primary field, field gradients, and 
25 RF field to rotate the magnetization vector such that it is perpendicular to coil 9. 

Tile KF coils may be two saddle shaped coils 9 which are driven in parallel to 
provide the rotating RF field. The Ftt) signals sensed arc received in tliis example by tlic 
three dimensional array of RF coils 30 and are ainpUlied by an RF amplifier 21 before 
being applied to signal handling circuits 22. The three dimensional detector array 30 is 
30 shown in more detail in FIGURE 12 as the tlvec dimensional detector anray 401 . The 
circuits 22 are arranged to make any appropriate calibrations and corrections but 
essentially transmit the signals to the processing circuits lo provide the required 
representation of the examined volume. Tliese circuits can conveniently be combined 
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with the circuits which control the primary field, field gradients, and RF field and thus 
are included in the ciicuils indicated at 20. The picture thus obtained is viewed on a 
display 23, such as a television monitor, and this may include inputs and other 
peripherals 24 for the provision of commands and instructions to the machine, or other 
5 forms of output. The display is not limited and includes any medium of conveing the 
image. Examples of displays, but not limited to, include printers, cathode ray tube 
displays, liquid crystal displays, plasma screens, three dimensional modelers, 
holographic displays, laser monitors, and projection monitors. 

The apparatus also includes field measurement and error signal circuits 25 which 
10 receive signals via amplifiers 26 fi-om field probes X^, X^, , and shown. 

The patient 27 is inserted in the tubular fonner of and coils 7, 8 and is 
supported there by a suitable couch or otlier supporting means. Such supports may be 
readily provided in any suitable form. 

The coils 7, S arc two sets of coils axlally displaced, each sei comprising two 
1 5 pairs of saddle coils the pair 7 being at 90** to Uie pair 8, These coils are themselves 
inserted into the central aperture in coils 6 which in an embodiment arc wound in 
four parts connected in series to provide an approximately circular configuration which 
is well known to be desirable for production of a uniform field. Further details of the • 
coil winding will not be given since suitable coils can readily be devised, by those willi 
20 the appropriate skills, to provide tJie fields required. 

The appropriate stores provide the amplitude and duration signals which are 
converted to analog forai in digital to analog converters (DACs) and applied to 
respective coil drive circuits x, y, z, RF. The respective drive circuits, which can take 
any form well known for driving field coils, provide the specified cunents to the 
25 appropriate coil for the specified duration. The apparatus and circuits described so far 
may be adopted to provide different gradients and RF fields, by appropriately adjusting 
the stored sequences and profile data. Similarly other known NMR apparatus which are 
capable of applying a steady magnetic field, a pulsed RF field, and , G,, and G, field 
gradients to a body, may be adopted in a manner straightforward to those skilled in the 
30 ^^MR art to apply the fields described for this invention. 

Gradients may be applied in any direction to further enhance the image 
reconshruction or image quality by methods known to those skilled in the art. hi an 
embodiment shown in FIGURE 8, the system componenls which provide imaging 
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enhancing gradients are 7, which provide C?,, the field gradient in the x-axis direction as 
shown, 8 wWch provide the G^, the field gradient in the y-axis direction as shown, and 
10, which provide G„ the field gradient in the z~axis direction as shown. For example, a 
magnetic field gradient along an axis in the xy-plane of the z-directcd primary field H^, 
5 is applied to produce a phase dispersion as a function of the distance along the axis, hi 
an embodiment, the phase due to tlie gradient may be linear and may be larger Uxm the 
phase of any Fourier component in the absence of a gradient During reconstruction, the 
phase infomiation due to the gradient may be applied to refine the assignment of the 
position of each voxel with respect to the distance along the gradient axis by methods 

1 0 known to those skilled in the art. 

With the basic signal handling system of the present invention, (lie FID signals 
from the signal sensing coils of the detector array 30 (shown in more detail in FIGURE 
12 as coils 403 of the tliree dimensional detector array 401) arc amplified in an Rl- 
amplifier and applied via an analog to digital converter (ADC) to a store such as a 

1 5 random access memory (RAM). The data is then processed according to the procedure 
given in the Reconstruction Algoritlun Section. 



Phase Angle and Associated Set of Foiuier Components 

20 

In an embodiment, tlie NMR signal at each detector as a function of time is 
processed by a method such as a Fourier transform operation to give a plurality of 
Fourier components each having an intensity, a phase angle, and the same frequency. 
The NMR signal of each voxel at any giv^ detector gives rise to a Fourier component 

25 with a unique phase angle relative to the Fourier component of any other voxel of the 
phantom at that detector. Tlie set of Fourier components that correspond to the NMR 
signal of a given voxel over the detectors is determined. This may be achieved by using 
a componttit with a defined standard phase relative to a first detector. Over the an:ay of 
sampled points, the phase angle of components at each detector may be converted to the 

30 corresponding phase at the position of Uic component having the standard phase. The 
components that have the standard phase are associated into a set of components 
comprising the components Irom a given voxel over the sample space. 
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In an embodiment, the phases of the components from the voxel are determined 
at the other detectors using ttie relative position of other detectors of the sample space 
relative to the first detector and the standard phase of the component at the first detector. 
In other words, the set of associated components may also be dctcnnincd by 1 .) 
5 identi^ing a first component having a characterizable phase angle, 2.) calculating tlie 
phase angle as a function of spatial position of the other detectors relative to tlie first 
using the phase angle of the component at the first detector, and 3,) identiiying the 
component at each detector having the calculated phase angle. These components 
identified from the data are associated to form the set. 
10 In the present invention, the phase variation of the NMR signal is provided by tlie 

combination of 1.) the angle 9 suspended between the direction of the detector and tJie 
radial vector, the vector fix)m the dipole to the detector, and 2.) the angle ^ due to a 
separation distance r between a voxel and a detector given by the wavenuinber of tlie RF 
field k times r . The distance r of the detector fi-om the voxel gives rise to the phase 
15 term e''*^ of the component of the detected RF signal where k is the wavenumber of ll\e 
NMR signal. Tlie harmonic oscillation of each RF dipole is equivalent to the dipole 
■ rotating in the transverse plane- The detector is responsive to a component in this plane. 
At a point in time, each RF dipole is directed at an angle 0 relative to the direction of 
detection of the detector. The angle $ of the RF dipole relative to the direction of 
20 detection axis of the detector gives rise to a phase angle term e'^'^ . In a preferred 

embodiment, the sum of the phase angles, Ar and d, are unique for each voxel at each 
detector. The position of each detector relative to a different detector may be used to 
calculated liie phase angle of tlie second relative to the first. Hiis may be repeated over 
all of the detectors to give the set of intensities of the NMR signal over the sample space 
25 due to each voxel. 

In an embodiment, tlie voxels that are on the same axis with a given detector or a 
plurality of detectors that align vertically with one or more voxels is detemiined. The 
phase angle is tlicn only a ftjnction of the distance r of the detector from (he voxel. Such 
components may be identified by the presence of at least one component of a phase 
30 given by ^ + where is the phase angle of the first component and Ar is the 
distance of the second detector relative to the first of the detector. This dctector{s) is 
defined as having a zero angle d . In this case, the phase angle ^2 ^ second component 
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at a second detector aligned on the detection axis corresponding to a first component at a 
first detector having a phase angle i^, is given by 



From the components identified as coaxial with a detector, the phases of the components 
5 from the voxel are determined at the other detectors using the relative position of other 
detectors of the sample space relative to tlic first detector and <j>^ of flic component at the 
first detector. 

FIGURE 14 shows the stationary coordinate system (x, y, z) of the NMR detector 
conresponding to the coordinate system of FIGURE 9 and FIGURE 10 of a primary field 

10 Hq oriented parallel to the z^^-axis and the z-axis and the coiresponding transverse RF 
magnetic dipole oriented in the x^y^ -plane and periodically parallel to the y-axis 
wherein the spatial variation of the RF y-componcnt of tlic flux due to the RF dipole is 
measured at a detector according to Eq, (I). FIGURE 14 further shows the distances and 
angles betv\'ecn a voxel 801 linear to a first detector 802, a second nonlinear voxel 800, 

15 and a second detector 803 in accordance with the invention. The NMll signal ofeach 
voxel at any given detector gives rise to a Fourier component v/ith a unique phase angle 
relative to the Fourier component of any other voxel of the phantom at that detector. In 
the present invention, the unique phase variation of the NMR signal is provided by the 
combination of I .) the angle 0 suspended between the direction of the detector and the 

20 radial vector, the vector fi:om the dipole to the detector^ and 2.) the angle ^ due to a 
separation distance r between a voxel and a detector. The pliase ^ due to a separation 
distance between a voxel ajid a detector of r is given by tlie wavenumber of the RF field 
k times r. 



In tenns of the NMR coordinates of the detector shown in FIGURE 1 1 and FIGURE 14, 
the phase angel (f> is given by 



25 



i^^kr (19) 
where the wavenumber k is given by 



(ao, 




(21) 
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where Ax , Ay^ and Az is the sqjaration distance between a voxel and a detector in the 
x,y,^nd z directions, respectively. In the case that a two dimensional slice is imaged 
in the transverse plane, the phase angel ^ is given by 

^=A^ = f VW + (A>'y (22) 

5 All voxels which have the same radial distance from a given detector will have 

the same phase given by Eq. (19), However, at each case that the radial distances r of 
two voxels to a given detector are equivalent, a unique angle 9 suspended between the 
direction of the detector an(i the vector along r exists. In an embodiment, the RF dipoles 
of all voxels are time synchronous. The ^MR pulsing and detection may be 

1 0 synchronized. At least one standard dipole may be used as reference phase to set the 
phases of the RF dipoles of the phantom. Al! voxels that are not collinear with the axis 
of a detector will posses an angle 6 suspended between the direction of the detector and 
the radial vector, the vector from the dipole to the detector. Since the RF dipoles are 
time synchronous, at each point in time» tliis angle corresponds to a phase angle 0 of the 

1 5 RF signal from each voxel at each detector. In the case that a two dimensional slice is 
imaged in the transverse plane, the phase angle 0 is given by 

0 Ax ^ sin-' f ..sin-^-— ^ (23) 

The total phase angle of the RF signal from each voxel at each detector is given as tlie 
sum of ^ and 6. 

20 0^=^ + 0 (24) 

The uniqueness of tlic phase angle a Fourier component of the NMR signal of 
each voxel at any given detector relative to a Fourier component of any other voxel of 
the pliantom at that detector is demonstrated by the following cases according to the 
coordinates shown in FIGURE 14: 

25 

Case I Ax ~ 0. 1 cm Ay = 0.1 cm 

0=O.78rad ^ = 5.9 X\0'' rod 

30 ^ = 0 + ^ = 0.78 + 5.9 X 10"" rad = 0.7806 rad 
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Case2 Ax^O.l cm ==^20 cm 

= 0 + ^ = 5 A- 10"' + 8.377 10'' ra^ = 8.88 A' 1 0'' ra J 
Case3Ax = 20cm A)/=O.Icm 

7(20)' + {0.1/ A ^ 

e=l.56rad ^ = 8.37 AT 10"' ra^/ 

= 0 + = 1 .56 + 8.377 Z J 0"' rarf - 1 ,643 ra</ 

Case 4 Ax = 20 cw Ajy^ = 20 cm 

15 

0 - sin-' , ^ = a , y 

0 =0.785 rnrf ^=0.1184 rae/ 

+ 0.785 +0.1 184 rarf = 0.9034 raJ 
Tlic condition for uniqueness of the phase #y ofeach voxel at a given detector 
20 depends on the unique sum given by Eq. (24). Since all voxels which have the same 
radial distance from a given detector will have the same phase ^ given by Eq. (19), but a 
unique angle 9 , each voxel will have a unique phase when the following condition is 
satisfied: 

the maximum angle 9 , is greater than <^,^ . the maximum angle for any 

25 r, 

^ou«>*m.x (25) 

^mw corresponds (o and , the minimum separation of the voxel and tlic detector 
in the x and y directions, respectively, and is given by Eq. (23). 
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corresponds to , the minimum separation of the voxel and the detector in the x 
direction and , the maximum separation of the voxel and the detector in the >r 
direction, and is given by Eq. (22), 

5 =yV(^7+o!Z7 (27) 

Substitution of Eq. (26) and Eq. (27) into Eq. (25) gives the condition for uniqueness of 
the phase angie ^f. 

>0 For ;'«x»*™,Eq. (29) gives 
X : 2k 



For = . Eq.(30) gives 

A>2;rV2j'„. (31) 

■^>8-9)'„», (32) 
1 5 In terms o f frequency / , Oq. (32) gives 

where c is the speed of ligbt. For = 20 c/w, the condition for uniqueness of the 
phase angle is / =< 168 A///z. 

Consider a component of a voxel at a first detector with a phase ^, defined as a 

20 standard phase. The corresponding phases of the components from the same voxel are 
determined at other detectors using the relative position of the other detectors of the 
sample space relative to first detector. In the case that a voxel is collinear with the first 
detector, the phase angle 0=0. ITie phase angle of tlie corresponding component at 
a second collinear detector can be calculated from tlie phase ajiglc ^, of a first using tlic 

25 detector separation along the collinear axis. In the case that a two dimensional slice is 
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imaged and the second detector is displaced Ay in the y direction of the sample space 
rciative to the first detector with a phase angle ^„ ttie angle is given by 

*2 = V(*Ayy + ikAyf = V(Myy Hi (34) 
Consider the phase angle ^„ of a corresponding component at a second detector 
5 that is nonlinear with the first detector. The phase angle ^ of the second component 
due to the same voxel may be calculated from In this case a contribution exists from 
I .) the angle Q suspended between the direction of the detector and the radial vector, the 
vector from the dipole to the detector, and 2.) the angle ^2 ^ separation distance r 
between a voxel and a detector given by the wavcnuraber of the RF field k times r , In 
1 0 the case that a two dimensional slice is imaged and the second detector is displaced Ax' 
and Ay in the x and y directions of the sample space, respectively, relative to the first 
detector with a phase angle the phase 0 is given by 

I Ax" . Ax* - »i ^Ax* 
— =sm '-T = sm ^ -T 

^ V(Ax'y +(Ayy V(^y + (wy 



(35) 



— sm 



The angle ^2 is given by 



15 = -^{kd^i + {kL/f + (Mj;/ = ^Ac*)' + (My / + ^\ (36) 

The phase angle of the second component due to the same voxel is given by the sum 
ofEq- (35) and Eq. (36). 

^„ =sin-' . . , ^^(k^f^ikAyf H^ (37) 

The component at the second detector which has phase angle is associated with the 
20 component ^, from which ^^.^ was calculated. AH phase angles at all other detectors 
are determined from each standard phase angle This is repeated for all using those 
detectors which are detcmiined to are collinear witli a voxel to give the sets of 
components each comprising the spatial variation of the RF field over the sample space 
due to each voxel. 
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Reconstruction Algorithm 



When a static magnetic field Bq with lines in the direction of the z-axis is appHed 
to an object comprising a materia! containing nuclei such as protons that possess 
5 magnetic moments, the field magnetizes the material. As a result a secondary field 
superposes the apphcd field as shown in FIGURE 9. In the applied magnetic field, the 
magnetic moments of each nuclei precesses about the applied magnetic field. However, 
the magnetization of any one nucleus is not observed from the macroscopic sample. 
Rather tlie vector sum of the dipole moments fix)m all magnetic nuclei in the sample is 

10 obser\'ed. This bulk magnetization is denoted by the vector M. In thennal equilibrium 
with the primary field H^^, the bulk magneU'zation M is parallel to H<^. The volume to 
be imaged is divided into volume elements called voxels, and the magnetized voxel 302 
shown in FIGURE 9 with bulk magnetization M is modeled as a magnetic dipole ra . 
Consider the case wherein data comprising the z-componcnt of the magnetic field of a 

1 5 dipole oriented in the z-direction is acquired by detectors 301 in the three dimensional 
sample space comprising the xy-planc and the positive z-axis as shown in FIGURE 9. 
The magnetic moment, , of each voxel within the phantom is a magnetic dipole. And, 
the phantom can be considered to be a three-dimensional array of magnetic dipoles. At 
any point extrinsic to Ihc phantom, tlic z-component of tlic secondary flux, B\ from any 
20 single voxel is 



where x, y, and z are the distances from the center of the voxel to the sampling point It 
is shown in APPENDICES MV that no geometric distribution of magnetic dipolcs can 
give rise to Eq. (i). Therefore, the flux of each magnetic dipole (voxel contribution) 
25 serves as a basis element to form a unique reconstruction of the array of dipoles which 
comprise tlie bulk magnetization map or NMR image of the phantom. 

Tiq, (1) is a system function which gives tlie magnetic flux output in response to a 
magnetic dipole input at the origin. The pliantom is an array of spatially advanced and 
delayed dipoles wcigliied according to the bulk magnetization of each voxel; this is the 
30 input function. The secondary flux is the superposition of spatially advanced and 
delayed flux, according to Hq. (1); this is the output function. Tlius, tJie response of 
space to a magnetized phantom is given by the convolution of Eq. (1) with the series of 
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weighted, spatially advanced aiid delayed dipoles representing the bulk magnetization 
map of the phantom. TTie discrete signals are recorded by a detector array over the 
sample space comprising (he xy-planc and the positive z-axis. 

In an embodiment of Uie present invention, the magnetization vector is rotated 
5 into the transverse plane by an additional RF field B , , The magnetization vector then 
comprises a rotating magnetic dipole m in the transverse plane. The NMR image may 
be reconstructed by samph'ng tlie external field from a series of RF dipolcs rather than 
(hat from a scries of static dipoles. In tl)is case, the geometric system function is givai 
also given by Eq. (1 ), the function of the z-component of the flux B' from a z-oriented 
1 0 dipole at any point extrinsic to the image space, from any single voxel. The geometric 
system fiuiction corresponding to a dipole which determines the spatial intensity 
• variations of the RF field is a band-pass for = as shown in AITENDIX V, In the 
h'mit, each volume element is reconstructed independently in parallel with all other 
volume elements such that the scan lime is no greater than the nuclear free induction 
1 5 decay (FID) time. 

In an embodiment, the strength and duration of the rotating H, (RF) field that is 
resonant with the protons of the magnetized volume and is oriented perpendicularly to 
the direction of the magnetizing field is applied such that the final precession angle of tJie 
magnetization is 90*^ (^^ = 90^ ) such that the RF dipole is transverse to the static 
20 magnetizing field and perpendicular to the RF magnetic field detector 20a as shown in 
FIGURE 1 1 . The rotating NMR coordinates and the stationary NMR cooidinates are 
both shown in FIGURE 1 1 . According to Eq. (1), the signal as a function of time which 
is Fourier transformed arises from each transverse RF magnetic dipole oriented in the 
XrYu -plane (rotating NMR coordinates) wliich is periodically paraHcl to the y-axis and 
25 rotates about the primary field H^, oriented parallel to the z-axis at the Larmor frequency 
(stationary NMR coordinates). For tlie stationary NMR case, the y-axis corresponds to 
the z-axis of Eq. (1), and , the magnetic moment along the z-axis, of Eq. (1) 
corresponds to the bulk magnetization M of each voxel. 

In addition^ each Fourier component comprises an additional part that 
30 corrcspomis to an RF magnetic dipole which is periodically pai-allcl to the x-axis. This 
part also rotates about tlie primary field Ho oriented parallel to the z-axis at tlic Larmor 
frequency wherein llic x-axis of the stationary NMR cooidinates corresponds to the x- 
axis of Eq. (1). The parts rotate at the same frequency but are orthogonal. The 
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ortliogonaKty corresponds to a phase angle of the Fouiier transfomi of 90**; thus, each 
Fourier component contains a real and an imaginary part. The magnitude of the 
amplitude of the signal due to the dipolc oriented parallel to the y-axis (z-axis of Eq. (1)) 
exceeds that of the dipole oriented parallel to tlie x-axis (x-axis of Eq. (1)) which 
5 identifies this part of a Fourier component. In an embodiment, the spatial variation over 
the detector array of the part of the Fourier component with the maximum amplitude (the 
signal due to li)c dipole oriented parallel to the y-axis) is used to determine the 
coordinate location of each voxel using the geometric system function of the detection of 
the z-component of a z-oriented dipole (Eq. (1)) as given in the Reconstruction 

10 Algorithm Section. 

In another embodiment, tlie spatial variation over the detector array of the lesser 
magnitude orthogonal part of each Fourier component is used to determine the 
coordinate location of each voxel using the corresponding geometric system function as 
described in the Reconstruction Algorithm Section. The system function for the case of 

15 the minor orthogonal part in terms of the coordinates given in FIGURE 9, conesponds to 
tfie detection of the x-component of a z-orientcd dipole. The geometric system function 
is given by the function of the x-component of the flux B* from a 2-oriented dipole at any 
point extrinsic to the image space, from any single voxel: 

20 where x, y, and z are the distances from the center of the voxel to Ihe sampling point and 
is the magnetic moment along the z-axis corresponding to the bulk magnetization M 
of each voxel. 

In tlie reconstruction process described herein, the secoudai^ field may be in the 
same or in a transverse orientation relative to the primary field. In each case, the 

25 orientation of the secondary dipole field and the measured secondary dipole field 

component are according to Eq. (1). The relationship of tlie coordinate systems of the 
present invention are shown in FIGURES !>, 10, and 1 1. The primary and secondary 
fields arc parallel and stationary in FIGURE 9 versus transverse in the rotating and 
stationary NMR coordinates shovm m FIGURES 10, 1 1, and 12. The relationship of the 

30 coordinate system of the reconstruction method of the NMR system used herein to that 
of tlic coordinate system of Eq. (i) shown in FIGURE 9 is according to FIGURE 1 1. 
ITic magnetization axis of Eq (I) is the z-axis; whereas, the magnetization axis of the 
stationary NMR coordinates is the z-axis shown in FIGURE 1 0, and the magnetization 
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axis of the NMR rotating coordinates [3,7] is the -axis shovvoi in FIGtBEUS J 1 . 
Regardiog Eq, (1) and (he reconstruction algorithjn in the NMR case, the z^component of 
Eq, (1) is substituted with the stationary NMR y^component, the y^component is 
substituted with the stationary NMR z-component, and the x-coraponent is substituted 
5 with the stationary NMR x-component. 

The reconstruction algorithm can be a reiterative, a matrix invereion, or a Fourier 
Transfonn algorithm. For all reconstrucHon algorithms, the voJurae to be imaged is 
divided into volume elements caUed voxels, and the magnetized voxel with magnetic 
moment is modeled as a magnetic dipole. In an embodiment, the matrix of Fourier 
10 components that correspond to the NMR signal of a given voxel over the detectors is 
determined. The matrices are determined for all of the voxels. The measurements of the 
spatial variations of the transverse RF field of a given matrix is used to determine the 
coordinate locaUon of each voxel. Thus, each matrix of components associated by phase 
comprises the intensity variation over the sample space of the RF field of the bulk 
15 magnetization M of each voxel. The matrices are the input for the reconstrucHon 

algorilhni. This procedure is performed either in parallel or series for each matrix. Tho 
bulk magnetization map (NMR image) is the superposition of the independent images, 
each of which corresponds to a given voxel. The supeqjosition of images is plotted and 
displayed. 

20 Aft embodiment of a matrix inversion reconstruction algorithm comprises the 

steps of 1.) using the geometric system function to determine the spatial intensity 
variation of the transverse RF field over the detector array, 2.) inverting the 
conrcsponding matrix, and 3.) multiplying the signal over the detector array by the 
mverted matrix to give the voxel sources. For example, a matrix inversion 
25 reconstruction algorithm is to detcmiine a coefficient for each voxel raathemalically (Eq. 
(41)) or by calibration which when multiplied by the bulk magnetization M of each 
voxel is that voxeFs contribution to the signal at a given detector with the corresponding 
um'que phase at each detector. Tliis is repeated for every detector, and those coefficients 
are used to det^ine a matrix which, when multiplied by a column vector of the bulk 
30 magnetization M values of the voxels, gives the voltage signals at the detectors. This 
matrix is inverted and stored in memory. Voltages as a function of time arc recorded 
over the detector anay. The signal as a function of time is Fourier transfomied to give 
the Fourier components each having an amplitude and a unique phase. In an 
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embodiment, of a matrix inversion reconstruction algorillun, the components are 
rauItipKed by the inverse matrix, to generate the bulk magnetization M map. 

In another embodiment of a matrix inversion reconstruction aJgorithm, the matrix 
of Fourier components that correspond to Ihe NMR signaJ of a given voxel over the 
5 detectors is determined. The matrices are determined for all of the voxels. The 
racasuremcats of the spatial variations of the transverse RF field of a given matrix is 
used to determine the coordinate location of each voxel. Thus, each matrix of 
components associated by phase comprises the intensity variation over the sample space 
of the RF field of the bulk magnetization M of each voxel. Each matrix of components 
10 are multiplied by the inverse matrix, to generate the bulk magnetization M map. In one 
embodiment, the point spread of the reconstructed voxel is corrected by assigning one 
voxel above a certain threshold with the bulk magnetization M . The other voxels are 
assigned a zero value. This procedure may be repeated for all voxels. In the limit with 
sufficient phase resolution, each volume element is reconstructed independently in 
15 parallel with all other volume elements such that the scan time is no greater than the 
nuclear free induction decay (FID) time. The total bulk raagjietization map (NMR 
image) is the superposition of the separate maps for eadi magnetic moment which is 
plotted and displayed. 

In an embodiment of a reiterative reconstruction algoritlun, the geometric system 
20 function.is used to determine the system of linear equations wliich gives the intensity, 
spatial variation, and phase of the RF field over the sample space. The signal as a 
function of time is Fourier tran&fomicd to give the Fourier components each having an 
amplitude and a unique phase. The system of linear equations gives the voltage and 
phase fi-om each voxel at each sensor based on the bulk magnetization M value of each 
25 voxel and the position of the voxel relative to the sensor. In an embodiment, weighting 
coefficients are determined based on these equations. The coefficients may be 
detennined mathematically (Eq. (41)). alternatively, they may be determined by 
calibi-ation. Tlic bulk magnetization M for each voxel is estimated, and Uie signals at 
each detector are calculated at each phase. In an embodiment, the bulk magnetization M 
30 value of each voxel times its weighting coefficient and its calculated phase at a given 
detector is compared to the measured voltage and pliase. A correction is made to IVf of 
each voxel. This gives rise to a second, or recomputed, estimate for M of each voxel. 
The signal value from this second estimate is computed and corrections arc made as 
previously described. This is repeated until the correction for each reiteration 
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approaches a predefined Kmit which serves lo indicate that the reconstruction is within 
reasonable limits of error, lliis procedure is repeated for all saisors. The final bulk 
magnetization map is plotted and displayed. 

The general process of reconstruction by reiteration is shown according lo Uie 
5 steps of FIGURE 1 3 (and is implemented in processor 20 in FIGURE 8). The image 
displayed according to the process 200 is merely a mapping of the buHc magnetization 
M of voxel sections of the object examined. The signal as a function of time is Fourier 
transformed to give the Fourier components each liavmg an ainphtude and a unique 
phase. The bulk magnetization M at a given location determines the intensity and 

10 spatial variation of the RF field over the sample space. Accordingly, signals produced 
by the RF magnetic sensors 1 10, in terms of volts as a function of phase, are a direct 
result of the bulk magnetization M of the voxel elements. Therefore, a reference voltage 
is generated at 210 from which the actual or measured sensor voltage at a given phase is 
subtracted at 220. The reference voltages arc modeled by assuming a voltage with a 

1 5 miique phase from each voxel at each sensor. Tlierefore, tlie voltage and phase from 
each voxel is determined at a sensor according to a weighting based on the position of 
the voxel relative to the sensor. The weighting may be given by Eq. (1) in terms of the 
coordinates x, y, z of FIGURE 9. The phase may be determined as given in the Phase 
Angle and Associated Set of Fourier Components Section. The resulting modeled or 

20 calculated voltage signals are compared at step 220, providing a difference or A(x,y,z) 
signal, weighted at step 230 to produce a weighted difference signal, which is then added 
to the previously estimated bulk magnetization M value for each voxel clement al step 
240. The resulting level, available in three dimensions corresponding to the axes x, y, 
and z, is selectively displayed on the display at step 250. Having adjusted the estimated 

25 bulk magnetization M for each voxel, the calculated bulk magnetization M is 

recalculated at step 260, the resulting estimated sensor voltage at each pliase is then 
compared to the actual sensor voltage at each phase at step 220, tlie process 200 being 
repeated until the difference is reduced to a tolerable value. This procedure is repeated 
for all sensors. The final bulk magnetization map is plotted and displayed. 

30 In terms of the coordinates of Eq. ( I ) and FIGURE 9» Uic reconstruction 

algorithm using Fourier Transforms (FFT) involves exploiting the FIT to solve Eq. (4 1 ) 
given below. For the case tliat follows, data is acquired in tlie x, y, and z-dircctious, but 
in geuCTal, data is acquhed over the dimensions which uniquely dctennine the bulk 
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magnetization map (NMR image). Also, the present analysis is for measuring the z- 
coniponent of the magnetic field of adipole oriented in the z-direction; however, Uie 
analysis applies to the other two orthogonal components where the geometric system 
function for the z-component of the z-oriented dipole is rq)laced by the geometric 
5 system ftmction for the x or y-componenls of the magnetic field produced by the dipole 
where the geometric system function is defined below as the impulse response of the 
detector to the given component of the field of a dipole of given orientation. The sample 
space, or space over which the secondary RF field is measuied, is defined m the present 
example as the three-dimensional space comprising the entire xy-plane and the positive 
10 2-axis, as shown in FIGURE 9. Other sample spaces are valid and each requires special 
consideration during the reconstruction as described below. The discrete voltages 
recorded fiom an infinite detector array in the xy-plane which is translated to infinity 
along the z-axis starting fi-om the origin where the detector array is responsive to the 
component of tlie secondary magnetic field is given by Bq. (1), where Uic voltage at any 
15 point in space produced by dipoles advanced m the z-direclion and advanced or delayed 
in the x and y-directions relative to the origin is given by the following Eq. (40), where 
the voltage response is Q times the secondary magnetic flux strength in the case shown 
in FIGURE 9. 

hi a preferred embodiment, the data comprises the RF field in the transverse 
20 plane over the sample space shown in FIGURE 1 1 . The relationship of Ihe coordinate 
designations of Eq. (1) and FIGURE 9 versus the NMR system and the relationsliip 
between the magnetic moment of a static secondary field and the bulk magnetization M 
of the RF field is given in the Reconstruction AlgoriUun Section and in tlie Fourier 
Transform Reconstruction Algorithm which follows. 

25 

Fourier Transform Reconstruction Algorithm 

In terms of the coordinates of Eq. (I) and FIGURE 9, tl)e volume to be imaged is 
30 divided into volume elements called voxels and the magnetized voxel is modeled as a 
magnetic dipole . Data comprising the z-component of fiie magnetic field of a dipole 
oriented in the z-direction is acquired in tlic three dimensional sample space comprising 
the xy-plane and the positive z-axis. 
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The phantom can be considered to be a three-dimensional array of magnetic 
dipoles. At any point extrinsic to the phantom, the z-componcnt of the secondary flux. 
B , from any single voxel is 

5 where x, y, and z are the distances from the center of ttie voxel to the sampling point. It 
is shown in APPEmiCES MV that no geometric distribution of magnetic dipoles can 
give rise to Eq. (1). Therefore, the flux of each magnetic dipolc (voxel contribution) 
fomis a basis set for Ihc flux of the array of dipoles which comprise the bulk 
magnetization map of the phantom. 
10 Eq, (1) is a system function which gives the magnetic flux output in response to a 

magnetic dipoJe input at the origin. Tlie phantom is an array of spatially advanced and 
delayed dipoles weiglited according to the bulk njagnctiz^ation of each voxel; this is the 
input function. The secondary flux is the superposition of spatially advanced and 
delayed flux> according to Eq, (I); this is the output function. Thus, the response of 
15 space to a magnetized phantom is given by the convolution of Eq. (I ) with tlie series of 
weighted, spatially advanced and delayed dipoles representing the bulk magnetization 
map of the phantom. The discrete voltages recorded by a detector array over the sample 
space comprising the xy-plane and the positive z-axis arc 

(41) 

where the flux magnetizing each voxel is unity, tlie volume element is made unity, and 
the voltage response of each detector is Q times the secondary magnetic jQux stiength. 
llie variables of Eq. (40) ai e defined as follows: 

Xr^ ^,.«, the bulk magnetization of the voxel located at the 

position defined by the Dirac delta function. 



20 



k^y k^, ky 

30 



dipole spacing in the x. y, and z-dircctions, 
respectively 
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the dimensions of the image space, Uie total distance in the x, y, and z-directions, 
respectively, for which the bulk magnetization of the voxels is nonzero 

the detector spacing in the x, y, and z-dircctions, respectively 

5 

The voltage signals recorded at the detector array over tlie sanaple space is 

at at CO 

(42) 

Eq. (41) can be represented symbolically as follows 
10 s = Q[gx[h®/\xuiz)] (43) 

where Q is the proportionality constant between the signal voltage and the output flux 
strength, s is the discrete function of voltage signals recorded of the secondary flux over 
the sample space, g is the secondary magnetic flux sampling function given as follows: 

eo «0 <o 

Z 'VpJ' - w,^ ,2 - 1 (44) 

15 In Eq. (42), h is the system function which is also defined as the geometric system 
function given as follows: 

The system function is the impulse response of the detector array whore tlic external 
magnetizing flux is set equal to unity. Fn Hq. (42), / is the bulk magnetization fimction 
20 or NMR image function given as 

In Eq. (42), w(z) is the unitary z function which is one for positive z and zero otherwise. 
'J'he function g dtscretizes the continuous voltage function V given by Eq. (40) which is 
the function h convolved with the fimction / then multiplied by the function u{z), Tlie 
25 discrete voltages recorded over the sample sj^ace are used in a Fourier transfonn 
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algorithm derived herein to reconstruct the NMR image. Consider the fiinction s of Eq. 
(41) which is giv<ai as follows: 

s^{h®f)xu(z) (47) 
Eq. (46) is equivalent to the function h convolved with the function / then multiplied 
5 by the function m(z). The function S which is the Fourier transform of ^ is given as 
follows: 

S^(H>cF)<E>U(k^) (48) 
'Vht fiinction S is equivalent to the resultant function of the function Hy the Fourier 
transform of the system function- A , multiplied by the function F, the Fourier transform 
10 of the bulk magnetization function- /, convolved with the function U(kX Fourier 
transform of the unitary z function- m(z) . The Fourier transform of the bulk 
magnetization function— / (Eq. (45)) 



(49) 



is 



15 ^= £ E Z -''^ (50) 

where jc„ = njk^\ y„ ^ nj(^\ z„ --nyK • Fourier Irarisfonn of 
m(z) = ifor z > 0 and u(z)^ Iforz <0 is [S] 

uiO-{6khj (51) 

The Fourier transfonn of tlie system function (Sec APPENDIX V) 
20 ,^^i^£zJ^^J^ (52) 



IS 



25 where 
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X 

A:=2«/;=2k- (55) 

y 

K='2nf,=2n- 

Z 

The resultant function of the product of the functions // (Eq. (52)) and F (Eq. (49)) tJicn 
convolved with the function (Dq. (50)) is 

•i, , , 2 /-I Zj An, 

5 Hie result is given as follows, and tlic derivation appears in AI*PEND1X VL 

(57) 

Substitution of the Fourier transform of the system function, (Hq. (52)), the Fourier 
transform of the bulk magnetization function> F (Eq. (49)), and factoring out in the 
10 second term gives 

£- y y V g -/iM- J 

f 1 



(58) 

Multiplication by the complex conjugates gives 
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^ '.-0 J-.— /,/2 JI,=-I,/2J: 



Art) /», .nj*' 



f 









n-/ 



I 

J J 



(59) 



Factoring out the common terms gives 

5 Combining exponential terms gives 
S 



(60) 



(61) 



Thefimction S divided by the fimclion // is 



O Q 



+ c**^'"' - 



2. <0 



(62) 



1 0 The inverse Fourier transform of the iiinction S divided by the function // is ^iven as 
follows, where the symbol B'''(0 is defined as the inverse Fourier transform of the 
function Q. 



The inverse Fourier transform of the function F is 



}; z.<0 (63) 



J 



15 i^'(/=-)=/=Z E E ;t...../[^-^..y-^..--^J 



(64) 
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Specifically, the inverse Fourier transform of the function F in terms of the dipole 
spacing is 



(65) 



Consider the general Fourier transform pair in cylindrical coordinates: 

Under circular symmetry, lhat is when g is independent of 0 (and hence G independent 
of ^ ) the inverse Fourier transform is [9] 

A(r,z) = 2«:JJ//(j,<s)X(2;wry^'^ja:y^G) (67) 

The inverse Fourier transform of Eq. (62): 



given by Eq. (64) and Eq. (66) is 



J 



• „^ A 



j+l,f* W,ftt I i J •IP^'J P P 



15 



(69) 

The solution of Eq. (68) appears as follows, and the derivation appears in APPENDIX 
VII. 
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.^i.-O ,.--6/2 Jt.-4i2 L 



5[z} 



(70) 



Inverse Transfonn 



5 Thus, the inverse Fourio- transform of — (Eq. (62)) is given by Eq. (69). 

H 

-^-ii Z Z Z X,.,,,e-^^'-''''-'''4e^^^^^ 

L^'«-0 >',i^-/,/2 x.-/,/2 L -fj. 



1 0 The convolution replaces each Goordinate with a spatially shifted coordinate. 

.r^i n ^ z ^ 

The result of the evaluation of Eq. (71) at each coordinate x^,y„,z„ is 

L//4,,.,_,./^—[2^.'J 

Solving for the bulk magnetization of eadi dipole gives 



(71) 



(72) 



(73) 



15 



L//J, 



(74) 



2^* 



llie solution of the bulk magnetization of each dipolc follows fiom Eq. (73). 
Discrete values of (he voltages produced at the detector array due to the secondary 
magnetic field are recorded during a scan which represent discrete values of function 4- 
(Eqs. (41-42)); tlius, in the reconstruction algorithm that follows, discrete Fourier and 
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Inverse Fourier transforms replace the corresponding continuous fimclions of Bq. (73) of 
the previous analysis. 

Discrete values of // of Bqs. (52-54), the Fourier transform of the system 
fiinction, replace the values of tlie continuous function. Furthennore, each sample 
5 voltage of an actual scan is not truly a point sample, but is equivalent to a sample and 
hold wliich is obtained by inverting the grid antennae matrices or which is read directly 
from a micro-antennae as described in the Finite Detector Length Section. The spectrum 
of a function discretely recorded as values, each of which is equivalent to a sample and 
hold, can be converted to the spectrum of the function discretely recorded as point 
10 samples by dividing the former spectrum by an appropriate sine function. This operation 
is performed and is described in detail in the Finite Detector Length Section. From these 
calculated point samples, the bulk magnetization Amction (NMR Image) is obtained 
following the operations of Eq. (73) as given below. 

15 Reconstruction Algorithm 

1 ) Record the RF NMR signal at discrete points over the sample space. Each point is 
designated (x, y, z» RF) and each RF value is an element in matrix A. 

2) 20 hivert the detector grid matrices defined by the noncoinmon areas of the overiapping 

elements of the detector airay described in the Finite Detector Length Section to obtain 
the sample and hold voltages which form Matrix A* (if micro-antennas are used, form 
Matrix A* of the recorded voltages directly). 

3) 25 Fourier transfonn the time dependent signals to give the intensity and phase of each 

component, 'Vhc NMR signal of each voxel at any given detector gives rise to a Fourier 
component with a unique phase angle relative to the Fourier component of any other 
voxel of the pbantooi at that detector. The matrix of Fourier components that correspond 
to the NMR signal of a given voxel over the detectors is determined. This may be 
30 achieved by using a first component having a phase angle and calculating the phase angle 
as a fujiction of spatial position of the first detector relative to aaiy otiier detector and 
identifying the component at each detector having the calculated phase angle. The 
matrices A are determincti for all of the voxels. The measurements of the spatial 
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variations of the transverse RF field of a given matrix is used to determine the coordinate 
location of each voxel. Thus, each matrix of components associated by phase comprises 
the intensity variation over the sample space of the RF field of the bulk magnetization M 
of each voxel. 

5 

4) Three-dimcnsionaiiy Fourier transfonn each matrix , using a discrete Fourier 
transform formula such as that which appears in McC. Siebert [10] to obtaiji each matrix 

of elements at frequencies coiresponding lo the spatial recording interval in each 
direction. 

5) Multiply each element of each matrix by a value which is the inverse of the Fourier 
transfonn of a square wave evaluated at the same fi*equency as the element where the 
square wave corresponds to a sample and hold operation performed on the continuous 

1 5 voltage function produced at the detector array by the secondary field. This 

multiplication forms each matrix * which is the discrete spectrum of the continuous 
voltage function discretely point sampled (See the Finite Detector Lengtli Section for 
details of this operation). 

smc(A:,,A^.,A:j 

20 

6) Multiply each element of each matrix * by the value which is the inverse (reciprocal) 
of the Fourier transform of the system function evaluated at the same firequcncy as tlie 
clement to form each matrix 

25 

7) Inverse three-dimensionally Fourier transform each matrix ^„ * * using the discrete 
inverse Fourier transform formula such as that which appears in McC. Siebert (10] to 
form each matrix whose elements correspond to the bulk magnetization of tlie dipoles 
at the points of the image space of spatial interval appropriate for the frequency spacing 

30 of points of matrix ^„ * * . 
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8) Divide each clement of each matrix C by — r to correct for the restriction that the 

sample space is defined as z greater than zero ( z > 0). This operation creates each matnx 
/), which is the bulk magnetization M map (NMR image). 

5 Q^^oSx.y.z) 



22 



In the NMR case shown in FIGURES 8, 10, 11, and 12, Step 7 of the 
Reconstruction Algorithm of the Fourier Transfonn Reconstruction Algorithm Section is 
10 equivalent to tlie general case except for the coordinate designations and (he relationship 
between the magnetic moment of a static secondaiy Held and the bulk magnetization M 
of tlie RF field whereui , the magnetic moment along the z-axis, of Bq. (1) 
corresponds to the bulk magnetization M of each voxel: 

8*55 Divide each element of each matrix C, by to correct for the restriction that the 

2yn 

sample space is defined as y greater than zero, y > 0. Tliis operation creates each matrix 
Z) which is the bulk magnetization M map. 

In one embodiment, the point spread of the reconstructed voxel is corrected by assigning 
20 one voxel above a certain threshold with tlie bulk magnetization M . The other voxels 
are assigned a zero value. This procedure may be repeated for all voxels. In the limit 
with sufficient phase resolution, each volume element is reconstructed independently in 
parallel v.nth all other volume elements such that the scan time is no greater than tlie 
nuclear free induction decay (FID) time. 

25 

9) Superimpose the separate bulk magnetization M maps corresponding to each matrix A,^ 
Plot the bulk magnetization M with tlie same spatial intoval in each direction as the 
sampling interval in each direction over all matrices A„ (i.e..tbc total bulk magnetization 
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M map which is plotted and displayed is the superposition of the separate maps of each 
magnetic moment corresponding a matrix A J, 

(The above steps for each matrix relate generally to the program implementation shown 
5 in the listing of the Exemplary Reconstruction ?ro^[^m as follows. The above Steps 1 
and 2 relate to the Data Statements beginning at lines 50; and Step 4 relates to the X, 2 
and Y FFT operations of lines 254, 455 and 972, respectively. Steps 5 and 6 are 
implemented by the processes of lines 2050, 2155 and 2340; and Step 7 relates to tlie X, 
Y and Z inverse transform of lines 3 170, 3422 and 3580, respectively. Step 8 and 8* . 
10 relates to the correction and normalization process of line 4677,) 

Exemplary Reconstruction Program 

10! 4D-MRI ALC30RITHM 

25 OPEN "4D-MRIPR0TOTyPEI.LIS" FOR OUTPUT AS #1 
15 26 C«.05 

27 DIM X(9,9,9) 

28 FOR W:=^l TO 9 STEP 1 

29 FOR T=l TO 9 STEP 1 

30 FOR S=l TO 9 STEP 1 
20 31 X{W,T,S}=0 

32 NEXT S 

33 NEXT T 

34 KEXT W 

35 X(5,5,6)=l 
25 36 DIM DI(9,9) 

37 PRINT »1, "DIPOLE PHANTOM" 

40 FOR Q«l TO 9 STEP 1 

41 FOR R=l TO 9 STEP 1 

42 FOR a=l TO 9 STEP 1 
30 43 LET H=X(U,R,Q) 

44 LET DI(U,R)=H 

45 NEXT U 

46 NEXT R 

47 t4AT PRINT 81, DI, 
35 48 NEXT Q 

SO DATA 0,0,0,0,0,0,0.0,0,0,0,0,0 

52 DATA 0,0,0,0,0,-.6, -.9,-1.1,-1.2, -1.8,5,7,3,5, .4,-.4,-.S,-. 7, -.6, .7 
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54 DATA 4.2,19.7,7.8, .1, - .9, -.3, -.4, -.7, ,7,13.9,25.65.10.1, .9, - .8,-. 3, - 
.6 

56 DATA - .4,2.7,10.7,12.4,5.6,2.6, -.5, -.3,-, 5, -.6, -.6.. 3, 3.7,. 9,-. 5,- 
1.0 
5 58 DATA 

0,0, 0,0, 0,0, 0,0, 0,0,0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0.0, 0,0, 0,0. 0.0, 0,0 
60 DATA 0, 0,^.5, -.8, -.8, -1.0, -1.2, -.8, 6, 9, .2,-. 2,-. 4,-, 6, - 
,6,1.2,4.1,12.7,0, .6 

62 DATA --6,-.3,-.5,-.7, -.3,3.7,9.0S,7.6,1.3,-.S, -.3, -.4, - 
10 .2,2.7,4.3,10.6, .8 

6« DATA 1.7,-.4,''.3,",5,-.4,-.l,l.l,l,l,l,7, .2,- 
.6,0,0,0,0,0,0,0,0,0,0,0,0,0 

66 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-.5, -.6, - .6, - .5, - 

.5 

15 68 DATA .3,2.5, .1,-.!, -.3, -.4, -.3,1.3,3.9,0.1,1 .7, .9, -.3, -.3, - .4. -.6, - 
.8,1.5 

70 DATA 1-65, 4. 7, 1.3, -.3, -.2, -.2,, 1,2. 1,3. 5, 7.1,. 7,. 7,-. 3,-. 2,-. 4,- 
.2, .2,1.2 

72 DATA 1.1, 1.4,. 4,- 
20 .3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 

74 DATA 0,0,0,0,0,0,0,0,0,0,-.4,-.5,0,-.l.-.l,1.7;i.5,.2,-.l,-.2,- 
.2,0,1.3,3.1 

76 DATA 4.1, 3.4, .9,0, -.2. -.3, - .4, -.4,1.1, .65,2. 9, 1. 3. -.2, - .1, - 
.1, .2,1,2,1.6 

25 78 DATA 4 .9, .1, -5, -.3, -.2, -.3, -.1,0, .3,1.3, .6, .1, - 
,1,0,0,0,0,0, 0,0,0,0,0,0,0 

80 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,^.3,- 
.2, .1, .9,1.6,1.6 

82 DATA . 4, 0, 0, 1, -. 1, .2,1.0,2.2,2.7,3.7, .5,0, -.2,-. 2,-. 2,- 
30 .1, .6, .25,1.6, .7 

84 DATA -.2,-.l,-.l,.l,.7, .9,2.8,0,.3,-.2,-.3,-.2,-.l,0, .2,.3,.3,O,- 
.l/O,0,0 

B5 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
86 DATA 0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, -.2,- 
35 .1, .1,1.0,1,4,1,3, .2,0,0 

88 DATA -.1,0,.2,.7,1.5,2.0,2.9,.3,0,-.1,-.1,.1,0,.7,.25,.9,.4,- 
.1,0,0, .1, .5 

90 DATA .5, 1.3. 0, .2, - .1, 0, - .1, 0, . 1, .2, . 1, .2,0. - 
.1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
40 92 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,- 
.1,0, .1, .8, . 9, 1.0, .1.0,0 
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94 DATA 0,O,.2,.5,l.a,1.3,2.O,.3,O,0,-.l,0.O, .3,.15,,6,.2,- 
.l,0,0, .1, .4, .3, .9 
$6 DATA 0, .2, -.1^0, - 

.1, .1, .1, .1, .1, .2, 0,0, 0,0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
5 98 DATA 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,- 

.1,0,0, ,7, .7, .1,0, .1 

100 DATA 0, .1, .2, .3, .6, .8,1,0, . 2,0,0, - 
.1,0,0, .1, .15, .3, .1,0,0,0,0, .2, -2, .5 
102 DATA 0, .1, - 

10 .1,0,0,0, .1, .1, .1, .1,0,0,0,0.0,0,0,0,0,0,0,0,0,0.0,0,0,0,0,0 
104 DATA 

0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, .5, .4, .4, .1,0,0,0,0, .1, .1 
106 DATA .4, .5, .6, .1,0, 0,0, 0,0,0,. 05, .1,0,0,0,0,0, ,1, .1, .3,0, .1, - 
.1,0,0,0.0,0 

15 108 DATA .1,0,0,0,0,0,0,0,0,0,0,0.0,0,0,0,0,0,0,0,0,0 
120 DIM V(9, 9, 9) 
130 FOR Z=l TO 9 STRP 1 
140 FOR 0=1 TO 9 STEP 1 
14 5 FOR U=l TO 9 STEP 1 
20 147 REAP V(Z,Q,U) 

222 NEXT U 

223 NEXT Q 

224 NEXT Z 

225 PRINT ttl. "VOXiTAGE DATA" 
25 226 DIM V0{9,9) 

227 FOR Q=l TO 9 STEP 1 

228 FOR R=l TO 9 STEP 1 

229 FOR U=l TO 9 STEP 1 

230 LET Ii»V(U,R,Q) 
30 231 LET VO(U,R)«H 

232 NEXT U 

233 NEXT R 

234 KAT PRINT Si, VO , 

235 NEXT Q 

35 254 IFFT THE ROWS OF SAMPIrED VOIiTAGBS IN THE X DIRECTION 

255 DIM MR (9) 

256 DIM MI (9) 

257 DIM R(9) 

250 DIM RV{9,9,9) 
40 259 DIM IV(9,9,9) 

260 FOR V=0 TO 9 STEP 1 
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270 FOR M-1 TO 9 STEP 1 
2 BO FOR N=l TO 9 STEP 1 
290 LET y-V(M,N,V) 
300 LET R(»)=y 
5 305 NEXT H 

310 EXTERNAL SUB FFT(DIM() ,DIM() ,DIM() } 
320 CALL FPT(R(>,MR(),MI(>) 
330 FOR N=l TO 9 STEP 1 
340 LET Y=MR(N) 
10 350 LET RV(M,N,V)=Y 
360 NBXT N 

370 FOR 1J=1 TO 9 STEP 1 

380 LET Q=MT(N) 
•390 LET IV(M,N,V)=G 
15 400 NEXT N 

430 NEXT M 

435 KfBXT V 

417 DIM RVO(9,9) 

41'9 DIM IV0(9,9) 
20 420 PRINT "RV" 

421 FOR Q^l TO 9 STEP 1 

422 FOR R«l TO 9 STEP 1 
4 23 FOR U=l TO 9 STEP 1 
424 LET K=»RV(a,R,Q) 

25. 425 LET RVO{U,R)-H 
426 UEKT U 
4 27 NEXT R 

428 MAT PRINT »1, RVO, 

4 29 NEXT 0 
30 430 PRINTT "IV 

431 FOR Q«l TO 9 STEP 1 

4 32 FOR R=l TO 9 STEP 1 

4 33 FOR U-1 TO 9 STEP 1 

434 LET H«IV(U,R,Q) 
35 435 LET XVO(U,R)=.H 

4 3« NEXT U 

437 NEXT R 

4 3S MAT PRINT itl, IVO, 
439 NEXT Q 

40 4 55 ?FPT THE COLWTNS OF THE SAMPLED VOLTAGE IN THE Z DIRECTION 
460 DIM RRV(9,9, 9) 
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470 DIM IRV(9,9,9) 
475 FOR V=l TO 9 STEP 1 
480 FOR W-1 TO 9 STEP 1 
490 FOR N=l TO 9 STEP 1 
5 500 LET Y=RV(M,H,V) 
510 liET R(M)=y 
520 NEXT M 

530 CALL PFT(R() ,KR() .MIO > 
540 FOR M=l TO 9 STEP 1 
10 SSO VST H=MR<M) 

560 LET RRV(M,N,V}=H 
570 NEXT M 

580 FOR M=l TO 9 STEP 1 
590 LET G=MX<M) 
15 600 LBT IRV(M,N.V)=G 
610 NBXT M 
620 NBXT N 

625 NBXT V 

626 DIM RRVO(9,S) 
20 627 DIM IRV0(9,9) 

630 PRINT #1,"RRV« 

631 FOR Q=^l TO 9 STEP 1 

632 FOR R=l TO 9 STEP 1 
63 3 FOR U =1 TO 9 STEP 1 

25 634 LET H=RRV{U,R,Q) 
637 LET RRVO{U,R)=H 
636 NKXT U 

639 NEXT R 

640 MAT PRINT #1, RRVO, 
30 641 NEXT Q 

650 PRINT #1, "IRV" 

651 FOR Q=l TO 9 STEP 1 

653 FOR R=l TO 9 STEP 1 

654 FOR U=:l TO 9 STEP 1 
35 655 LET K« IRV(U,R,0> 

656 LET IRVO{U,R>-H 

657 NBXT U 

658 NBXT R 

659 MAT PRINT #1, IRVO, 
40 660 NEXT 0 

661 DIM RIV(9,9,9) 
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662 DIM IIV{9,9,9> 
68$ FOR V«l TO 9 STEP 1 
690 FOR N=l TO 9 STEP 1 
700 FOR M«l TO 9 STEP 1 
5 710 LET y=IV(M.N.V) 
760 LET R(M)-y 
770 NE;tT M 

830 CALL PFT(R(),MR(),MI()) 
B40 FOR M=l TO 9 STOP 1 
10 850 LBT H=»m(M) 

860 LET RIV(M,N,V)=H 
870 NEXT M 

672 FOR M= 1 TO 9 STEP 1 
873 LBT G<»MI(M) 
15 874 LKT 3IV{M,N,V)=G 

875 NEXT M 

876 NEXT N 

877 NEXT V 

870 DIM RlVO(9,3) 
20 879 DIM IIVO(9,9) 

880 PRINT Si, "RIV" 

881 FOR Q=l TO 9 STEP 1 

882 FOR Rrl TO 9 STEP 1 

883 FOR Ur-.i TO 9 STEP 1 
25 884 LET H=RIV(U,R.Q) 

BBS LBT RlVO{U,R)«H 

886 NEXT U 

887 NEXT R 

888 MAT PRINT ill, RIVO, 
30 889 NEXT Q 

890 PRINT ffl, "I IV" 

891 FOR Q»X TO 9 STEP 1 

892 FOR Rs=l TO 9 STEP 1 

893 FOR U=»l TO 9 STEP 1 
35 894 LBT H«IIV(U, R,0) 

895 LET IIVO(U,R)«H 

896 NEXT U 

897 NEXT R 

898 MAT PRINT ftl, IIVO, 
40 8 99 NEXT Q 

900 DIM RVA(9,9,9) 
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901 DIM IVA(9.9,9) 
904 DIM RVAO(9,9) 
906 DIM IVAO(9,9) 
908 FOR Q-1 TO 9 STEP 1 
5 910 FOR R*^l TO 9 STEP 1 

911 FOR U»l TO 9 STEP 1 

912 LET H=IIV(U,R,Q) 

913 H=<-1)*H 

914 LET G=«RRV(0,R,Q) 
10 915 LET L«G+H 

916 LET RVA(U,R,Q)==L 

917 NEXT U 

918 NEKT R 

919 NEXT Q 

15 920 FOR Q^l TO 9 STEP 1 

921 FOR R«»l TO 9 STEP 1 

922 FOR U=l TO 9 STEP 1 

923 LET H»IRV(U,R,Q) 

924 LET L=RIV(U,R,Q) 
20 925 LET IVA(U,R,Q)«H+L 

927 MKXT U 

928 NEXT R 
930 NEXT 0 

932 PRINT #1, "RVA« 
25 934 FOR 0=1 TO 9 STEP 1 
93G FOR R=l TO 9 STEP 1 
938 FOR U=l TO 9 STEP 1 
940 LET H*RVA(U,R,Q) 
942 LET RVAO(U,R)=H 
30 944 NEXT U 
946 NEXT R 

94 8 MAT PRlin* #1, RVAO, 

950 NEXT 0 

952 PRINT #1, "IVA" 

35 954 FOR 0»1 TO 9 STEP 1 
956 FOR R=l TO 9 STEP 1 
958 FOR U=l TO 9 STEP 1 
960 LET H»IVA(U,R,Q} 
962 LET IVAOCU,R)i=H 

40 964 NEXT U 
966 NEXT R 
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968 MAT PRINT «1. IVAO, 

970 NEXT Q 

972! PPT THE ROWS OF THE SAMPLED VODTAGBS IN THU Y DIRECTION 
98a DIM RVAy{9,9.9> 
5 990 DIM IRVAYO,9,9) 
992 DIM RVAYO(9,9) 
994 DIM IRVAY0(9,9) 

1012 FOR M«l TO 9 STEP 1 

1013 FOR N=l TO 9 STEP 1 
10 1014 FOR V=l TO 9 STBP 1 

1015 LET Y=RVA(M,N,V) 

1016 LET R(V)»Y 

1017 NEXT V 

1018 CALL FFT(R(),MR(),MI()> 
15 1019 FOR V=l TO 9 STBP 1 

1020 LET H^MR(V) 

1021 LKT RVAY <H,N,V)=H 

1022 NEXT V 

1023 FOR V=^l TO 9 STEP 1 
20 1024 LBT G=MI (V) 

1025 LBT IRVAY(M,N, V)-G 

1026 NEXT V 

1027 NBXT N 

1028 NEXT M 

25 1030 PRINT Si. "RVAY" 

1032 FOR Q=l TO 9 STEP 1 

1034 FOR TO 9 STEP 1 

1036 FOR U^^l TO 9 STBP 1 

103a LBT H=RVAY(U,R,0> 
30 1040 LET RVAYO(U,R)«H 

1042 NEXT U 

1044 NEXT R 

1046 MAT PRINT Sl.RVAYO, 

1040 NEXT Q 
35 1050 PRINT #1. -IRVAY" 

1052 FOR Q=l TO 9 STEP 1 

1054 FOR R=l TO 9 STEP 1 

1056 FOR U^l TO 9 STEP 1 

1058 LET H=1RVAYCU,R.Q) 
40 1060 LET IRVAYO|U,R)«H 

1062 NEXT U 
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1064 NEXT R 

1066 KAT PRINT #l,IRVAYO, 
1060 NEXT Q 
1070 DIM RlVy(9,9,9) 
5 1000 DIM IIVy{9r9,$) 

1085 DIM RIVY0{9, 9) 

1086 DIM IIVYO(9,9) 
1090 FOR M=l TO 9 STEP 1 
1100 FOR N=l TO 9 STBP 1 

10 113 B FOR V«l TO 9 STEP 1 

1139 LET y=IVA(M,N,V) 

1140 LET R(V)*Y 

1141 NEXT V 

1142 CALL FFT(R() ,MR() ,MI()) 

15 1143 ron v=i to 9 step 1 

1144 LET H=MR{V) 

1145 LET RIVY{M,W,V) =H 

1146 NEXT V 

114 7 FOR V=l TO 9 STBP 1 
20 n4 8 LET G=MI (V) 

1149 LET IIVY(M,N,V) =G 

1150 NEXT V 

1151 NEXT N 

1152 NEXT M 

25 1153 PRINT #1, "RIVY" 

1160 FOR 0=1 TO 9 STEP 1 

1162 FOR R=l TO 9 STBP 1 

1164 FOR U=:l TO 9 STEP 1 

1166 LEl^ K=RIVY(U,R,0> 
30 1170 LET RIVYO(U,R)=H 

1172 NEXT a 

1174 NEXT R 

1178 MAT PRINT #1. RIVYO, 

1180 NEXT 0 
35 1185 PRINT #1. "IIVY" 

1190 FOR O^l TO 9 STEP 1 

1200 FOR R^l TO 9 STEP 1 

1210 FOR U-1 TO 9 STBP 1 

1212 LET llVy(U,R,Q) 
40 1214 LET IIVYO{U,R)s«H 

1216 NEXT U 
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1218 NEXT R 

1220 PRINT ||1,IIVY0, 

1222 NEXT Q 
1230 DIM YRVA(9, 9,9) 
5 1240 DIM yiVA(9,9.9) 

1241 FOR Q=l TO 9 STEP 1 
1250 FOR R=l TO 9 STBP 1 
1260 FOR U=l TO 9 STEP 1 
1270 LET L-.IIVY{U,R,Q) 

10 1260 LET B=RVAy(U.R,Q) 

1290 LET yRVA(U,R,0)»B-L 
1300 WBXT U 
1310 NEXT R 
1320 MEXT Q 

15 1330 FOR Q-1 TO 9 STEP 1 
1340 FOR R-1 TO 9 STEP 1 
1345 FOR U«l TO 9 STEP 1 
3350 LET H=RIVY{U,R,Q) 
1360 LET L«IRVAy(U,R,Q} 

20 1370 LET YIVA(U, R,QJ=L+H 
1380 NEXT U 
1390 NEXT R 
1400 NEXT Q 
1410 PRINT ^1. "YRVA" 

25 1412 DIM YRVAO(9,9) 
1414 DIM YIVAO(9,9) 
1420 FOR 0=1 TO 9 STEP 1 
3430 FOR R=l TO 9 STEP 1 
1440 FOR U=l TO 9 STEP 1 

30 1450 LET H=:yRVA(U,R,Q) 
1460 LET YRVAO(U,R)=H 
14 70 NEXT U 
14 80 NEXT R 

3490 MAT PRINT #1,YRVA0, 

35 1500 NEXT 0 

1510 PRINT #1, -yiVA" 
1520 FOR Q=l TO 9 STEP 1 
1530 FOR R=X TO 9 STEP 1 
1540 FOR U=l TO 9 STBP 1 

40 1545 LRT H=YIVA(U,R,Q) 
ISSO LET YIVA0(U,R)=H 
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1550 NEXT U 
1570 NEXT R 

1580 MAT PRINT #l,yiVAO, 
1590 NEXT Q 

5 2050 »GENERATE TOE DISCRETE SPECTRW OP THE SYSTEM FUNCTION AND THE 
SINC 

2055 '.FUNCTION OF THE SAMPIJ3 AND HOLD CORRESPONDING TO THE FINITE 
DETECTOR 

2057 1 DIMENSIONS 
10 2060 DIM SPH(9,9,9) 

2061 DIM SFH0(9,9) 

2062 DIM SINC (9,9,9i 

2063 DIM SINCO(9,9) 
2065 LET SFH(4,4,4)«=4*PI 

J5 2066 FOR M==-4 TO 4 STEP 2 
2070 FOR Vc=-4 TO <i STEP 1 

2080 FOR TO 4 STEP 1 

2081 H=.04 

2082 J=ABS(M}+ABS(N}-i-ABS(V) 
20 2085 IF »J=0 THEN GO TO 209B 

2090 T=4*PI* ( (2*PI*N/9*l/.l) "24 ( (2*PI *V/9*l/ . I) "2 ) ) 
2092 B« {2*PI*N/9 * 1/ . 1) *2+ <2 *PT*M/9*l/C) ^2^ (2 *PI*V/9*l/ . 1) ''2 
2095 LET SFH(M+5,N+5,V+5)=T/B 
2098 G:=ABS (N) 4-ABS(V) 
25 2100 IF 0=0 THEN GO TO 2114 
2101 A«l 

2105 IF N^O THEN GO TO 2107 

2106 A=SIN(2*PI*10*N/9*.1}/(PI*N/9*10) 

2107 B=l 

30 2108 IF V=0 THEN GO TO 2110 

2109 B=rSIN(2*Pr*10*V/9*.l)/(PI*V/9*10) 

2110 H=^A*B 

2111 IF N=0 THEN H-.2*H 

2112 IF V=0 THEN H=,2*K 

35 2114 LET SrNC(M+5,N+5,V+5)=H 
2130 NEXT N 
2133 NEXT V 
2132 NEXT M 

2136 PRINT #1, "SFH" 

40 2136 FOR Q=l TO 9 STEP 1 

2137 FOR R=l TO 9 STEP 1 



wo 02n6956 PCT/USOl/25954 

70 

2138 FOR U=l TO 9 STEP 1 

2139 LET S=SFH(U,R,Q} 
214 0 liBT SF110(U,R)^S 
2141 NEXT U 

5 2142 NEKT R 

2143 MAT PRINT tfl, SPHO, 

2144 NEXT Q 

2145 PRIOT #1, '»SINC" 

2146 FOR Q=l TO 9 STEP 1 
10 2147 FOR R=3 TO 9 STBP 1 

2148 FOR U=l TO 9 STEP 1 

2149 LET S=SINC{U,R,Q) 

2150 IiET SXNC01U,R)=S 

2151 NEXT U 
15 2152 NEXT R 

2153 MAT PRINT #1, SIMCO, 

2154 NEXT Q 

2155 IIWVERSE THE DISCRETE SPECTRUT4 OF THE SYSTEM FUNCTION fiND THE SINC 
2157 ! FUNCTION 

20 2160 DIM HR(9,9,9) 

2170 DIM KR0<9, 9) 

2171 DIM SINCR(9,9,9) 
2175 FOR V=l TO 9 STBP 1 
2180 FOR M-l TO 9 STEP 1 

25 2190 FOR N=l TO 9 STEP i 

2200 LET Yc^SFH(M,N,V) 

2201 LET H=SINC(M,N, V) 

2210 IF Y=0 THEN GO TO 2212 

2211 v«i/y 

30 2212 IF H=0 THEN GO TO 2230 

2221 H=l/H 

22 3 0 LET HR(M,N,V)=Y 

2235 LET SINCR (M,N, V) =H 

224 0 NEXT N 
35 2250 NEXT M 

2260 NEXT V 

2310 PRINT #1, "HR" 

2311 FOR Q=l TO 9 STBP 1 

2312 FOR R«l TO 9 STBP 1 
40 2313 FOR U=l TO 9 STBP 1 

2314 LET S-lfR(U,R,Q) 
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2315 LET HRO{U,R)=S 

2316 NEXT U 

2317 NEXT R 

2320 MAT PRINT ftl,HRO, 
5 2321 NBXT Q 

2340 ! DIVIDE THE TRANSFORMED DATA BY THE TRANSFORM OP THB SYSTEM 
FUNCTION 

2345 I AND TUB SINC FUNCTION 
3030 DIM FR(S,9,9) 
10 3050 DIM FI(9,9,9) 
3052 DIM FRO (9, 9) 
3054 DIM FIO(9,9) 

3056 FOR V=l TO 9 STEP 1 

3057 FOR M=l TO 9 STEP 1 
!5 3058 FOR N=l TO 9 STEP 1 

3059 T=yRVA{M/W,V) 

3060 S=RR(M,N,V} 
3062 ri=SINCR(M,M, V) 
3065 Kt.S*T*L 

20 3066 hKT FR{M,)W,V)=K 

3067 NEXT N 

3068 NEXT M 
3070 NEXT V 

3080 FOR V=i TO 5 STEP 1 
25 3081 FOR M=l TO 9 STEP 2 

3082 FOR N»l TO 9 STEP X 

3083 H«SINCR(M,N,V) 
3093 K= HR<M,N,V) 
309*1 Ii= YJVA(M,N,V) 

30 3095 E«K*L*H 

3096 LET F1(M,N,V)=B 

3097 KEXT N 

3098 NEXT M 
3100 NEXT V 

35 3130 PRINT #1, "FR" 

3131 FOR Q=l TO 9 STEP 1 

333 2 FOR R=l TO 9 STEP 1 

313 3 FOR U«l TO 9 STEP 1 

3134 LET H=: FR(U,R,Q) 
40 3135 LET FRO(U,R)=R 

3136 NEXT U 
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3137 NEXT R 

3130 MAT PRINT #1, PRO, 
314 0 NEXT Q 
3141 PRINT #1, ™FI" 
5 3142 FOR 0=1 TO 9 STEP 1 

3143 FOR Ral TO 9 STEP 1 

3144 FOR U=l TO 9 STEP 1 

3145 LET H=FI(U,R,Q) 

3146 liBT PI0(U,R1=H 
10 3156 NEXT U 

3157 NEXT R 

3158 MAT PRINT #1. PIO, 
3160 NEXT Q 

3170 IINVERSB TRANSFORM THE SOWS IN THE X DIRECTIOT 
15 3180 DIM RP(9,9,9> 

3X85 DIM RFO(9,9) 

3187 DIM WlPO{9,9) 

3190 DIM MIF(9,9,9) 

3195 DIM 1(9) 
20 3196 FOR V=l TO 9 STEP 1 

3200 FOR M«l TO 9 STEP 1 

3210 FOR N-1 TO 9 STEP 1 

3220 LET Y-FR(M,N,V) 

3230 LET R{N) sY 
25 3240 NEXT N 

3250 FOR N:?l TO 9 STEP 1 

3260 LET Y=FI(M,N,V) 

3270 LET I(W)=Y 

3280 NEXT N 

30 3285 BXTERNAIi SUB IFT(DIM (> ,DIM{) ,DIM () , DIM () ) 

3290 CALL IFT(R() ,10 ,MR<) ,MI()) 

3300 FOR N=l TO 9 STEP 1 

3310 LET Y=HR(N) 

3320 LET RF(M,N,V)sY 
35 3330 NEXT N 

3340 FOR N=l TO 9 STEP 1 

3350 LET Y=MI(N) 

3360 LET MIF(M,N,V) '=Y 

3370 NEXT N 
40 33 BO NEXT M 

3385 NEXT V 
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3390 PRINT #1, "RF" 

3391 FOR Q*l TO 9 STEP 1 
3395 FOR R-1 TO 9 STEP 1 
3393 FOR U=l TO 9 STEP 3 

5 3394 LBT H«RF(U,R,0) 

3395 bET RFO(U,R)=H 

3396 NEXT a 

3397 NEXT R 

3400 MAT PRINT #1, RFO, 
10 3405 NEXT Q 

3410 PRINT jfl, •» MIP " 

3411 FOR 0=1 TO 9 STEP 1 

3412 FOR R=.l TO 9 STEP I 

3413 FOR U«l TO 9 STBP 1 
15 3414 LET H= KIF(U,R,Q) 

3415 LEI' MIFO(U,R)«H 

3416 NEXT U 

3417 NEXT R 

3420 MAT PRINT #X, MIPO, 
20 3421 NEXT Q 

3422 1 INVERSE TRANSFORM THB ROWS IN THB Y DIRECTION 
3430 DIM RFY(9,9,9) 

3432 DIM MIFY(9, 9,9) 

3433 DIM RFYO(9,9) 
25 3434 DIM MIFYO(9,9) 

3435 FOR M=l TO 9 STEP 1 
3440 FOR M=a TO 9 STKP 1 
34 50 FOR V=l TO 9 STEP 1 
3460 LET y«RF(M,N,V) 
30 3470 LET R(V)=Y 
3480 NEXT V 

3490 FOR V*l TO 9 STBP 1 
3500 LKT Y==MIF(M,N,V) 
3510 LET I(V)«Y 
35 3520 NEXT V 

3525 EXTERNAI, SUB IFT (DIM ( ) , DIM ( ) , DIM () , DIM ( ) ) 

3526 CALL IFT(R() , I () ,MR() ,MI () ) 

3527 FOR V=l TO 9 STEP 1 

3528 LET Y«MR (V) 

40 3529 LET RPY(M,N.V)=Y 
3530 NEXT V 
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3531 FOR V=l TO 3 STEP 1 

353 2 LET Y=KI (V) 

3533 LET MlPy (M,N, V)=y 

3534 NEXT V 
5 3535 MEXT N 

3536 NEXT M 

3537 PRINT #1, "RVY" 
353 D FOR Q^l TO 9 STEP 1 
3539 FOR R=l TO 9 STEP 1 

10 3540 FOR U=l TO 9 STEP 1 

3541 LET n^nFY{V,^,Q) 

3542 LET RPYO{U,R)=-H 

3543 MEXT 0 

3544 NEXT R 

15 3545 r-lAT PRINT #1, RPYO, 

354 6 NEXT 0 

3547 PRINT 1*1, •'MIPY" 
3550 FOR Q»l TO 9 STEP 1 
3555 FOR R=^l TO 9 STEP 1 
20 3560 FOR U=l TO 9 STEP 1 

3565 LET H«MIFY (U, R,0) 

3566 LET MIFYO(U,R)-H 

3567 NEXT U 
3566 NEXT R 

25 3570 MAT PRINT ftl^MIFYO, 
3575 NEXT Q 

3 580 i INVERSE TRANSFORM THE COLUMNS IN THE Z DIRECTION 
3581 DIM F(9,9,9) 

3590 DIM FO(9.9) 
30 3592 FOR V«l TO 9 STEP 1 

3593 FOR H=l TO 9 STEP 1 

3594 FOR W=l TO 9 STEP 1 
3600 LET Y«RFY(M,N.V) 
4546 LET R (M) =Y 

35 4 547 NEXT M 

4548 FOR M=^l TO 9 STEP 1 

4 54 9 LET Y:=MIFY {M,N, V) 
4 550 LET I(M)'=Y 

4 552 NEXT M 

40 4553 EXTERNAL SUB XFT2 (DIM () , DIM () , DIM () , DIMO ) 
4554 CALL IFT2 {R() , I () ,m () ,MI () ) 
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4556 FOR M^l TO 9 STEP 1 

4557 LET Y=MR<M) 
4559 LBT F(M,N,V)=y 
4600 NEXT K 

5 4602 NEXT N 

4604 NEXT V 

4605 PRIMT "F« 

4610 FOR Q-a TO 9 STEP 1 

4620 FOR R=l TO 9 STEP 1 
10 4630 FOR U-rl TO 9 STEP 1 

4635 LET H=P(U,R,Q) 

4640 LET PO(U,R)n=H 

4650 NEXT U 

4660 NEXT R 
15 4666 MAT PRINT tfl,FO, 

4670 NEXT Q 

4677 I CORRECT FOR THE U(Z> COMVOLUTIOW AND NORMALIZE THE RECOMSTRaCTlON 

4678 DIM CF(9,9,9) 
4780 Drr4 CF0(9, 9) 

20 4 809 FOR V=l TO 9 STEP 1 

4810 FOR N=l TO 9 STEP 1 

4871 LET K=Pfl,N,V) 

4B73 K^K/IOEB 

4974 LET CF{1,N,V)=K 
25 5915 NEXT N 

6616 NEXT V 

6617 FOR V=l TO 9 STEP 1 
6620 FOR M=2 TO 9 STEP 1 
6630 FOR N:=l TO 9 STEP 1 

30 6640 LET K=F(M,N,V) 

6642 P==0.5*PX/{M-1)"2/C"2 
6650 K«=K/P 

6660 LET CF{M,N,V)=xK 
6670 NEXT N 
35 6680 NEXT M 
6681 NEXT V 

6690 PRINT #1, "RECONSTRUCTION" 

6691 FOR Q«l TO 9 STEP 1 

6692 FOR R=>1 TO 9 STEP 1 
40 6693 FOR U«l TO 9 STEP 1 

6694 LET H=CP(U,R,Q) 
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€695 Hs=H/CF<5,5,6) 

6696 LET CPO(U, R>«H 

6697 NEXT U 

6698 NBXT R 

5 6700 MAT PRINT #1,CP0, 
6705 NEXT Q 
6710 END 

6722 SUB FFT(R() ,MR() ,MI() ) 
6726 FOR M=l TO 9 STEP 1 
10 6730 A«0 

6740 FOR N==l TO 9 STEP 1 
6750 LET H=R{N) 

6760 B=H*CX)S(2«PI*lM-5)*(M-5)/9) 
6770 A-A+B 

15 6780 NEXT N 
6790 A-A/9 
6000 LET MR(M):^A 
6810 NEXT M 

6820 FOR H=l TO 9 STEP 1 
20 6830 A=0 

664 0 FOR N=l TO 9 STEP 1 
6880 LET K=R(N) 
6890 H::^-H 

6900 B-H*SIN(2*PI* (M-S)*(M-5)/9) 
25 6910 A^A+B 

6920 NEXT N 

6930 A=A/9 

6935 LET MI (M) =:A 

6940 NEXT M 
30 6950 END SUB 

6960 SUD IFTmO ,I().NR() ,MI()} 

6970 DIM MRR(9) 

6980 DIM MRI (9) 

6990 DIM MIR(9> 
35 7000 DIM Mir (9) 

7010 FOR N«l TO 9 STEP 1 

7020 A=0 

703 0 FOR M=:l TO 9 STEP 1 

704 0 LET GsR(M) 

40 7050 B=G*COS(2*PI*(M-5)*(W-S)/9) 
7060 A=A+B 
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7070 NEXT M 

7 0 80 LET MRR(NJ«A 

7090 NEXT M 

7100 FOR N=l TO 9 STEP 1 
5 7130 A«0 

7120 FOR M=:l TO 9 STEP 1 
7130 LET G^R(M) 

7140 B==G*SIN(2*PI* (M-5) * {N-5) /9\ 
7150 A=A+B 
10 7160 NEXT M 

7170 LET MRI(N)=A 
7180 NEXT H 

7190 FOR N=l TO 9 STEP 1 
7200 A=»0 
15 7210 FOR K»l TO 9 STEP 1 

7220 LET G=I(M) 

7230 B=:G*C0S(2*PI*{M-S)*{N-5>y'9) 
7240 A=A+B 
7250 NEXT M 
20 7260 LET MIR{N)»A 
7270 NEXT M 

7280 FOR N=l TO 9 STEP 1 
7290 A«0 

73 00 FOR H-1 TO 9 STEP 1 
25 7310 LET G^I(M) 

7320 B«G*SIM(2*PT* (M-5) * (M-5) /9) 
7330 A=A+B 
734 0 NEXT M 
7350 LET KII(N)=A 
30 7360 NEXT N 

7365 MAT MII=(-1)*MII 
73 7 S MAT MR«MRR+Mri 
7385 MAT MI=MrR+KRI 

74 00 END SUB 

35 7410 SOB IFTZ(R(J ,I().MR() ,MI()) 

7420 DIM MRR(9) 

7430 DIM MRI(9) 

7440 DIM MIR (9) 

7450 DIM MI1(9) 
40 7460 l.X>R M«l TO 9 STEP 1 

7470 A=0 
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74 80 FOR M«l TO 9 STEP 1 
74 SO LET G=R(M) 

7500 B=G*C0S<2*PI* (H-5) * ( -N+X ) /9 ) 
7510 A=A+B 
5 7520 NEXT M 

7530 LET MRR(W)sA 
7540 mXT N 

7550 FOR N-:l TO 9 STEP 1 
7560 A=0 
10 7570 FOR M=l TO 9 STEP 1 
7500 LET G«R(M) 

7590 BcG*^SIN(2«PI*(M-5) *(-N+l)/9> 
7^00 A=A+B 
7610 NEXT M 
15 7620 LET MRI(N)«A 
7630 NEXT N 

7640 FOR N«i TO 9 STEP 1 

7641 AsO 

764 2 FOR M=l TO 9 STEP 1 
20 7650 LET G=X(M) 

7660 D=G*C0$(2*PI* (M-S) * (-Wr+l)/9) 

7670 A=A<-B 

768 0 NEXT M 

7690 LET MIR{N)=:A 
25 7700 HEXT N 

7710 FOR N=l TO 9 STEP 1 

7720 A=0 

7730 FOR M^l TO 9 STEP X' 

774 0 LET G=I(M) 
30 7750 B=G*SIN(2*PI*(M-5)*{-N+l)/9) 

7760 A=AfB 

7770 NEXT M 

7780 LET MII(>I)=A 

7790 NEXT N 
35 7795 MAT MII«(-1)*MII 

7800 MAT MR=:MRR+MII 

7810 MAT MI^MRl+MIR 

7880 END SUB 

$ 

40 
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In an embodiment, the matnx inversion and/or the reiterative algoritinns arc used 
in combination with the Fourier Transform Algorithm. For example, an NMR image 
created by the Fourier Transfonn Algorithm may used as input for the first reiteration of 
the reiterative algorithm. 

5 

The Nyquist Theorem With Hie Detennination Of The Spatial Resolution 

The derivation ofEq. (41) demonstrates that the system fiinction behaves as a 
filter of the spectrum of the bulk magnetization function (NMU image). It is well known 

10 in the art of signal processing tliat if such a filter passes all frequencies for which an 
input function has significant cucrgy, then the input function can be recovered 
completely from samples of the filtered function taken at the Nyquist rate. This premise 
embodies the Nyquist Sampling Theorem. The spectrum of the system function (Eq. (1)) 
is derived in Al^PENDIX V and shown in FIGURE Ic. This function is a band-pass for 

15 all frequencies of the bulk magnetization function where and A, are comparable. 
Thus, the bulk magnetisation function can be recovered by sampling tlie continuous RF 
field function given by Eq. (40) at the Nyquist rate, twice the highest fi-cqucncy of the 
bulk magnetization function, in each spatial dimension over the sample space for which 
the function has appreciable energy. Sampling operations otl^er than tlic present 

20 operation and Uie negligible ciror encountered by not sampling over the entire sample 
space are discussed in McC. Siebcrt (11] and the references therein disclosed which are 
all incoiporated herein by refcreuce. In the absence of noise, the spectrum of the bulk 
magnetization function can be completely recovered if the detector spacing frequency is 
equal to the Nyquist rate which is twice the highest ftcquency of the bulk magnetization 

25 function, and this represents tlie limit of resolution. However, the density of the detector 
. spacing is limited by noise. The three-dimensional bulk magnetization map is a 
reconstruction from independent recordings at independent detector spatial locations 
relative to the voxels of the image space where two detector signals are independent if 
they are sufficiently spatially displaced from each other such that the difference in signal 

30 between the two detectors due to a test voxel is above the noise level The resolution 
based on signal-to-noise arguments is discussed in the Contrast and Limits of Resolution 
Section. 
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Contrast and Limits of Resolution 

5 The ability to visualize a structure in a noisc-frce environment depends, among 

other factors, on the local contrast C, which is defined as 

(75) 

where / is the average background intensity and A/ is the intensity variation in the 
region of interest. The main source of >MR image (also called magnetic resonance 

10 images (MR!)) contrast is and which depend on tissue types. The contrast may be 
increased by using RF pulse sequences to polarize tlie protons which encode 
(spin/lattice) and 7^ (spin/spin) relaxation time information in the data of the secondary 
RF field in a manner straightforward to tliose skilled in the NMR art. 

Contrast, is not a fundamental limit on visualization since it can be artificially 

1 5 enhanced by, for example, subtracting part of the background or raising the intensity 
pattern to some power. Nojse» however, represents a fundamental limitation on the 
ability to visualise structures. The signal-to-noise ratio, a basic measure of visualization, 
is determined by the ratio of the desired Intensity variations to the random intensity 
variatioris, which are governed by the statistical properties of the system. The signal-to- 

20 noise nitio (SNR) is defijied as 

Sm-^- = C^ (76) 

where 0, is the standard deviation of the background intensity representing the rms 
value of the intensity fluctuations. The noise properties of the 4D-MRI imager involve 

additive noise only principally from theimal (y ) noise in the RF measurement circuits of 

25 the antennas of the detector array and from the fluctuations of the primary magnetic 
field. 

A feature of superconducting magnets is their extreme stability. The object to be 
imaged is magnetized with a hi^ily stable magnet such as a superconducting magnet In 
this case, a magnetic field r^ability of 10"* % over a month's time is feasible. Small 
30 antennas measure the RF signals as point samples without significant decrease in the 
signal to noise ratio relative to large ant^mas by using impedance matching while 
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minimizing resistive losses by usbg superconducting reactance elements, for example. 
In an embodiment, cross talk between antennas is ameliorated or eliminated by time 
multiplexing tlie signal detection over the array of antennas. Bxtemai sources of RF 
noise can be ameliorated by placing the 4D-MRI scanner in a shielded room (Faraday 
5 cage). 

The quality of the image (i.e. the signal to noise ratio of the image) can be 
increased by repeating the reconstruction over multiple time points wherein each data set 
of a given time point is the set of matrices of the intensity variation over the sample 
space of the RF field of the bulk magnetization M of each voxel having the 
10 corresponding magnetic moment at that synchronized (common) time point 

Sources of detector enror are random, and the noise averages out as the number of 
detectors increases. Typically, the noise is suppressed by a factor of the inverse square 
root of the number of detectors. Thus, the effective limit of RF field detection is 
increased by a factor of the square root of the number of detectors for a constant SNR. 
15 The resolution depends on the extent that the field of the ring, shell, or sphea^ of 

dipoles differs from that of a single dipole at the center. 'Vhe plot of the three cases of the 
field of a ring, shell, and a sphere of dipoles each of radius R and magnetic moment 
m = 10' Gem' given by Eq. (1.14), Eq. (0.11), and Eq. (1V.16) of APPENDIX I, 
APPENDIX n, and APPENDIX IV, respectively, as a function of radius R where the 
20 position of the center of the ring, shell, and sphere relative to the detector is the point 
(0,0, 1 0) is given in FIGURES 2, 4, and 6, respectively. The plot of the three cases of the 
field of a ring, shell, and sphere of dipoles of radius R = 0.2 cm and magnetic moment 
w = 1 0' Gem' given by Eq, (1 1 4), Eq. (II. 1 7), and Eq. (IV. 16) of APPENDIX I, 
APPENDIX n, and APPENDIX IV, respectively, as a function of the distance between 
25 the detector at the origin and tlie center of the ring, shell, and sphere at the points (0,0, 
z=A cmtoz-\5 cm)\s shovm in FIGURES 3, 5, and 7, respectively. From FIGURES 
2-7, it can be appreciated that the detector must be able lo resolve three to four 
significant figures in order to reconstruct a map of 0.2 cm resolution with a field depth of 
1 5 cm. 

30 

Finite Detector Length 

The system fimction, h , of the Reconstruction Algorithm Section is the impulse 
response for a point detector, llic following analysis vAW concern the impulse response 
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for an NMR detector which has finite dimensions. Consider a detector array 401 as 
shovjn in FIGURE 1 2 comprising multiple parallel planes 402 wherein each plane has a 
plurality of antennae coils 403 wherein each has an area A, The signal at any coil 403 is 
given by the integral of Eq. (40) over the area of the coil. This is the impulse response 
S which is the system function whicli replaces h for a finite area detector. The Fourier 
Transform of this system function contains an argument of a product of the detector area 
A and the spatial frequency variables. Reconstruction could be performed as previously 
described in the Reconstruction Algorithm Section where this system function is 
substituted for the system function for a point detector. In the limit of zero area, the 
10 measurement is that of a point sample. Thus, another approach is to use the linearity of 
the superposition of RF signal to gain a higher signal to noise ratio advantage by taking 
the difference of relatively large signals as opposed to performing the measurements 
with miniature antennas. In this case, the set of n coils each of area A of each plane 
comprise a grid of n blocks each of area a«A formed by the absence of common area 
1 5 overlap of one or more of the coils. The signal of each element of the grid of area a can 
be solved by the set of linear difference equations of the signals of coils each of area A 
that corresponds to the noncommon area for the overl^ptng coils. The solution of the 
signals due to the grid elements can be obtain using a computer by matrix inversion. 

The resulting values represent the average signal for each grid center location. 
The effect of this data processing operation on the spectrum can be modeled as a sample 
and hold, where the voltages at the centers of the grid elements arc sampled by 
multiplying by a picket fence of delta functions spaced s units i^art which are convolved 
with d , a square wave function in the x-dtrection and the y-direclion of width s units 
where the coordinates of Eq. (1) and FIGURE 9 are used. In the frequency domain, this 
data processing operation causes the spectrum of the signal function s to be multiplied 
by Dy the Fourier Transform of the square wave function of width s uriits, to fonn 
function 5* . If this multiplication does not multiply 5, the Fourier Traiisfonn of the 
signal function, s , by zero for any frequency less than its bandwidth, then S can be 
recovered from S* by multiplying iS* with the inverse of the Fourier transform of the 
sample and hold square wave ftinction, a two dimensional sine function for the x and y- 
directions. This analysis applies to all axes in which direction the detectors have finite 
length. Furthermore, z-sampliiig is achieved by translating tlic array in the z-direction by 
interval distances at which points discrete signals are recorded or by using multiple 
parallel plane detector arrays spsced at tlie sampling interval along the z-axis. However, 



wo 02/16956 PCTAJS01/259S4 

83 

if the signals are not sampled at discrete z-points, but each sample point is Ihe integral 
resultant of the signal acquired continuously over a z-displacement of q units wiiich is 
much greater than the dimension of the detector in the z-dircction, then the 
corresponding sample and hold square wave has a width of q units. 

5 

It will thus be seen that the invention efficiently attains the objects set forth 
above, among those made apparent from the preceding description. Since certain 
changes may made in tlie above constructions witliout departing finom the scope of tiie 
invention, it is intended that all matter contained in the above description or shown in the 
10 accompanying drawings be Interpreted as illustrative and not in a limiting sense. h\ 

particular, other methods such as electronic and optical methods of detection of magnetic 
resonance are within tlie scope of the present invention. 



While the clainaed invention has been described in detail and with refereiKe to 
1 5 specific embodiments thereof, it will be apparent to one of ordinary skill in the art that 
various changes and modifications can be made to the claimed invention without 
departing from the spirit and scope thereof. 
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APPENDIX I 

Integration of a RING of Dipoles 

Derivation of the Field Produced by a Ring of Magnetic Dipoles 

The z -coroponent of the magnetic field due to a dipole or a loop of current of 
radius R with dipole moment m = inl^ at the origin is given as follows: 



2z —X —y 
10 ^'^''^ItI 1 — '7^ 



a.i) 



The i -component of the magnetic field at the position (x,y,z) due to a ring of dipoles of 

radius R with dipole density centered at the origin is 

27tR 



^ ■(2(z-z'y-(x-yy-(K-y/>^^ 



wherein (xV/,z') is a variable which correspoiids to the position of each individual dipole 
15 of the ring. The relationship between Cartesian coordinates and cylindrical coordinates 
with z' = 0 is 



jc^-Z^cos^ (i-3) 
y = Rsin^ 

Substitution of Eq. (1.3) into Eq. (1.2) gives 

' " i - Rcosij^f Rsm + J' 

20 Multiplying out terras gives 

g '^ (2z' (y^ -f jg' cos^ ^ - IxRco^f)- -f/?^ sin' j> - lyRsmi^^J^^ 
' 2;t j ((x^ + /^^cos'<^^~2x/?cos(^)^(y^'^- J^'sin'<^^--2;//^sin^^^^ 

Substitution of Eq. (1.3) in the denominator and using the associative relationship and tJie 
trigonometric identity, cos^a +sin' a = I in the denomijiator gives 

^1^(2 2' - (r^ + R^ cos' ^ - 2jc/?cos^ )- + R'' sin^ - 2y/?sin^))/^ 
2 ;r (r' + - 2R{xcos^ + ;;sin 0)) 
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The denominator can be s^proximated usmg the following relationship 

(a+bYsa'+nba'-' (1.7) 

where 

a = r^ + ie (1.8) 

5 b = -2R(xcios^+ysia4') (1.9) 
and 

n = -5/2 (I.IO) 
This case gives the far field where a«b and/or r»Jt. Using Eqs. (I.7-I.10) gives 
(r' -2R{xcos<l> +;«in^))"'" 

1 5R , ^ . Oil) 

10 Substitution of Eq. (I.l I) into Eq. (16) gives 



w ^ 1 5/2 



»« x' A'cos'^ 2.xigcos^ 



2jr J! 



r 5Ax*cos<^ SjZx^ysin (j> S/?'xcos' ^ 1 
J| 5j?Vcos^^sin^ lOx^/?^ cos' ^ 10xyA'cos<^sm^ , 



d4> 



cos^ 5/?y sin^ cos ^sin^i^ 1 

2ji h S/fVsin'0 10/{^xvcos^sin^ 10/;VsiJi'^ i ^ 

The integration of Eq. (1.12) gives 

2;r)>' ;r/?' 10?ij?y , 



(112) 



(1.13) 
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Collecting terms gives 



JRingofDipoles: 



,f 1 



(1.14) 



10 • , 

Eq. (L 1 4) demonstrates that the magnetic field of a ring of dipoles is not equal to 
that of a single dipole at the origin. The ring radius, R , appears in ilie denominator of 
• each temr The first term is the magnetic field of a dipole at the origin only when the 
1 5 variable corresponding to tlie radius of the ring, R , is zero. The second antl third terms 
arc additional perturbations of the field of a dipole at the origin whose magnitude is a 
function of the radius of the ring. The second and third tenns vanish only when the 
radius of the ring is zero. A ring of dipoies has a field that is cylindrically syromctrical 
A shell, a cyhnder, and a sphere of dipoles are the only other cases which have this 
symmetry. A cylinder is a linear combinatioii of rings. Thus, the uniqueness of the 
dipole field is. demonstrated by showing that it is different fi*om that of a ring, a shell, and 
a sphere. The present result thai the field of a dipole is different for thai of a ring of 
dipoles as weli as the same result in the cases of a shell and a sphere of dipoles shown in 
APPENDIX II and APPENDIX IV, respectively, demonstrate that dipole field is unique. 
All other fields arc a linear combination of dipoles. Thus, tlie dipole is a basis element 
for the reconstruction of a NMR image. The resolution depends on the extent fliat the 
field of the ring of dipoles differs from that of a single dipole at the origin. The plot of 
the field of a ring of dipoles of radius R and magnetic moment m ^ 10* Gcm^ given by 
Eq. (L 14) as a function of radius R is given in FIGURli 2. The position of the center of 
the riffig relative to the detector is the point (0,0,10). Tlie plot of the field of a ring of 
dipoles of radius R = 0.2 cm and magnetic moment m = 10'* Gcm^ given by Eq. {1. 14) 
as a function of Uie distance between the detector at the origin and the center of the ring 
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at the poinis (0,0» z = 4 cm fo z = 1 5 cm ) is shown in FIGURE 3. From FIGURES 2 and 
3 it can be appreciated that the detector must be able to resolve three to four significant 
figures in order to reconstruct a map of 0.2 cm resolution with a field dcptJi of 15 cm. 
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APPENDIX n 

Integration of a Spherical Shell of Dipoles 

Derivation of the Field Produced by a Shell of Magnetic Dipoles 

The z -component of the inagnedc field due to a dipole or a loop of current of 
radius R with dipole moment m = in^ at the origin is given as follows: 

/ 2 . 2 2\5'2 1 (nj) 



10 B =/w 



The r -component of the magnetic field at the position (x^y^z) due to a shell of dipoles of 
radius R with dipole density centered at the origin is 

wherein (x^y^z") is a variable which corresponds to the position of each mdividual dipole 
15 of the shell. The relationship between Cartesian coordinates and spherical coordinates is 

x=psin0cos6^ a:'= /?sin^cos0 (113) 

;/= psin^sinO }f= Rsin^sinO 

Substitution of Eq. (U.3) into Eq, (11.2) gives. 

B ^^JH. 

(II.4) 

20 Multiplying out tenns gives 

U Tf (2z^ "4z/?cos^ -f 27g'cos'^-;c' -f 2.x:;?sin<^cose 

47tR^ jJ(x^-2x/esin^cosd4/?^sin^^cos^0+y-2y«sin^sinO 

-R^s'in^ <p cos' e -/ + 2y/? sin ^sing - sin^ ^ sin^ ^)/?^ sin ^0^0 
sin* ^ sin^ 0 + - 2z/?cos ^ + cos^ 

(11.5) 

Multiplying out the /f' sin ^ terra gives 



y^ ^(z-Rcos<j>y ~(x -Rsm(j}cosOf -(y^ iRsin <^ sing/ ^ sin j^cf^ifO 
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'^^litFjj (x^-2x/?sin^cosO+i?^sm^({>cos^0+y -2>/i?sm^sin0 

sin^ ^sin' 0 + - 2zi? cos^ + /J^ cos^ ^ 

(U.6) 

Substilution of Eq. (11.3} in the denominator gives 

rfi } (2 jg^z^sin ^ - AzR^ sin j)COs^ + 2^?"* cos^ <^ sin<^ - siji^ -f 2xR^ sb^ ^ cos Q 
'^4jiR^ II (p^ "ZxRsm^cosO R^ sin^ ^cos^ 9 - 2yRsin<}ts'mO 

-R"* sin' ^cos^ 0 - y sm<^ + 2>'if' sin' ^sm 0 - sin' ^sin^ 0)i^</e 
+/?^siu^ ^sin^ 0 -2zRco$^+ cos^ ^^^^ 

5 (n.7) 

Using the associative relationship and using the trigonometric identity, 
cos^ a h sin^ a - 1 in the denominator gives 

7} (2ig^z^sin ^ - 4zjg^sin^cos^ 2/?^ cos^ ^ sin^ - sirn^ A-lxR^ sin^ ^cosQ 
^' II (p' -2.vi^sin^cos0 

•-^^sin^ i^cos^ 0--//esin + 2>'jg' sin' ^sin 0-R' sin^ ^sin' 0)f^cfQ 
-2>'/t5in ^sin0 - 2z/?cos0y'^ 

(118) 

10 The denominator can be approximated using the following relationship 

(aibY^aUnba''' (11.9) 

where 

a^p'^R' (11-10) 
6 = -2x/esin^cos0-^2yjRsin^sin0-2z/?cosi^ • (K.Il) 

15 and 

« = -5/2 (n.i2) 
ITiis case gives tJie far field where a«b and/or p»R. Using Eqs. (11,8-11.12) gives 

{p^ ^R^ -2xi?sin<^cos0 -2y/?sin<^sin0 - 22/?cos(J)J^'^ 

1 57? , .v(II'13) 
+71 jrjTi (jrsin^cos0 +>?sin ^sin0 -f zcos^) 



Substitution of Eq. (11.13) into Eq. (U.8) gives 



2nn 

m 



2« «.=7fbfJ 
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2it*cos^^sin<i 10/?' cos' ^ sin . , ^ . , ♦ ^ ,\ 
+ ^ jpy^^ 4- jg^y^ ' ^ ^ cos 0 + ^sin sin 0 + z cos (p) 

~ ^^-f /g^) ~ (p^-hJ^y^^ cos 0 + ysin ^ sm 0 + 2 cos 

— n2\^^ / , ^,0/2 — (xsm^cos0+>'sm0sin0 4.zcos^) 

^p2 ^^^2y^^ " ^2^j^2y^^ ^ COS a + ySllKj^ Sin G 4- Z COS ^} 



(x sin ^ cos 0 + ysm ^ sin 0 + z cos if) 



^-—^ — — 71 — ^ .7/2 siM<^ cos e j^sin ^ sm 9 +zcos^) 



to The integral of a sum is equal to the sum of the integrals. 



'~ Anlt [j ^1 lef^ ^ (P^* +ysm4>smO hzcosi^JJrf^Q 



4zi?^siii^cosi 20ziJ*$ia^cos^ / . , ^ • . ,v i\ 
~F1 — ^,^5/3 +~71 — .oV>2 (^^^" ^<^os^ +ysm<^sinO + zcos^J 



d^de 



^fV 2i?^cos^^sin(A 10i?^cos^ ^sin^ , . . ^ - , - iS ,1 

**" il "~71 — ";x s /T ' ^ ~71 r37r-^(xsm^cos0 +;/sm^sine +zcos^m^y 

ooL \P "^^ ) (P ) 



J j 71 — '^^^7^7 — — Vcsin^cos0 + j^sm^sme + zcos^) |cf^9 

oo|_\P + J (P 

';7 2xJ?'sin^<^cosO lOx/2* sm^<^>cosO . . , „ . . - « i^ 
15 +] I h^Tl — + — 71 — ^TTtT? — (xsm^cosO+j;sin^sinO + zcos0j 

ooL (P "^'^ J \P ) . 



■(x sin^ cos<? + ysin (ftsinO+z cos^) 
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(jc sin ^ cos 0 + ysin ^sin 0 + z cos ^) 



i [' (p^' +it^y^^ (p^4- ig^y^ ' "^^ 9 sin 0 + z cos ^ ) |/0 



(n.l5) 



The integration of Eq. (11.15) gives 



m 



+ - 



AnR* 



AnR 



' 1 



(11.16) 



10 Collecting terms gives 



Shell of Dipoles: 

15 



20 Eq, QLH) demonstrates that tlie magnetic field of a shell of dipoles is not equal 

to lliat of a single dipole at tl)c origin, llie shelf radius, R , appears in the denominator of 
tlie first term. The first term is the magnetic field of a dipole at the origin only when the 
variable corresponding to the radius of the shell, R , is zero. The second term is an 
additional perturbation of the field of a dipole at tlie origin whose magnitude is a 

25 function of the radius of the shell. The second terra vanishes only when tlie radius of the 
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shell is zero. Thus, the dipole is a basis element for the reconstruction of a NMR image. 
The resolution depends on the extent that the field of the shell of dipoios differs from tliat 
of a single dipole at the origin. The plot of the field of a shell of dipolcs of radius R and 
magnetic moment w = 10* Gcfn given by Eq. (11.17) as a function of radius R is given 

5 in FIGURE 4, The position of the center of the shell relative to the detector is the point 
(0,0, 10). The plot of the field of a shell of dipoles of radius = 0.2 and magnetic 
moment m - 10^ Gan^ given by Eq, (nJ7) as a function of the distance between the 
detector at the origin and the center of the shell at the points (0,0, z) is shown in FFGURC 
5. From FIGURES 4 and 5 it can be appreciated that the detector must be able to resolve 

10 three to four significant figures in order to reconstruct a map of 0.2 cm resolution with a 
field depth of 15 cm. 
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APPENDIX in 



Proof tliat the Field Produced by a Shell of Magnetic Dipoles is Different from that of a 
5 Single Dipole 

Consider Eq. (11.5). 

^ m ^? f (22^ - 42/? cos (f> 4- 2R^ cos^ ^ - -f 2x/?sin<^ cos 0 
'^Ajzlf ilix^- 2xR sinjf cos 6 $in^ (pcos^O + - lyRsux ^sin 6 

-R^sm^<l>cos^ 0 -/ + 2yRsm^sm9 - sin^ (^sin^e)/?^ sin^^e/O 
-^R^ sin' sin^ 0 + - IzRcos^ + J?^ cos^ 

(fll.l) 

1 0 The conditions for this integral to equal that of the field of a dipole at the origin are 
-AzRcos^ + 2^?" cos' ^ + 2xRsm <pcos0-I^sin ^cos^ e+2yRsin^ sinO - sin^ ^sin^ 0=0 

(in.2) 

and 

-2xR^\n ^ cos 0 + R^ sin' ^ cos' 0 - 2yRsin ^ sin 0 + R^ sin' ^ sin' 0 - 2z/?cos <^ + cos* ^ = 0 
15 

(fflj) 

Thus, Eq, (01.2) must equal Eq. (1113), 
-42/?cos^ + 2i?^ cos^ -f 2x/?sin (/>cos0 - sin^ ^cos^ 0^2yRsm(^ sinO - i?^ sin^ ^sin^ 0 
^ -2x/?sin <f> COS0 + /?^sin^ <^ cos^ 0 - lyRsin^sinO -f R^ sin' ^sin^ 0 - 2zRcos^ + R^ cos^ ^ 

20 (ni.4) 
Therefore, 



-22/?cos^ + cos^ ^ - 2;?^ sin^ ^ cos^ 0 + 4jcsin^ cosO + 4yJ?sin^siii 0 - 21^ sin^ #sin^ 0=0 

011.5) 

25 Using Wig associative relationship and the trigonometric identity, cos^a +sin^ a = 1 as 
well as factoring out an 7? in Eq. (IJI.5) gives 

-2zRcos^ -f cos' ^ - 2R^ sin' ^ + 4;csin ^cos0 + AyRsin^ sin 0 = 0 (III.6) 
-2zcos ^i-3R cos' ^ - 2/? + 4jc sin^ cos 0 + 4>'sin ^ sin 0 =^ 0 (ID.?) 

30 
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Pot x^O, yjtQ^z^O.Jlq. (m,7) is true only ifR^O which proves that Ihe field of a 
shell of dipolcs centered on the origin is different froni that of a single dipole at the 
origin. 

Eq. (111.5) can not be intc^ted in closed form; however, the integral can be 
5 approximated so that the detector tolerances for a given image resolution can be 
determined. The denonrinator of Eq. (I1I.5) can be approximated using the following 
relationship 

(a+by^a^-^nba"" (m,8) 

wliere 

10 a^p'vR^ (ni,9) 

b = -IxRsiu ^cos 0 - 2yR sin ^ sin 0 - 22/2 cos ^ (ID. 10) 

and 

/i = -5/2 (ni.li) 

The result given by Eq. (11.17) is 



15 B =-^^4;^' 

4;r^ [ 



pr \} {ni.i2) 



Higher order terms of the approximation given by Eq. (111.8) would contain mixed 
products of the coordinate variables which would increase the deviation of the draivation 
of Uie field a shell of dipoles from that of a single dipole at the origin. 
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APPENDIX IV 
Integration of a Sphere of Dipoles 

5 

Derivation of the Field Produced by a Sphere of Magnetic Dipoles 

The z -component of the magnetic field due to a sphere of dipoles is derived 
using the equation for a shell of dipoles with the substitution of the dipole density of the 

10 sphere — -r where the radius of the sphere is R . The field of a shell of dipoles is given 
47ric 

approximately by Eq. (D. 1 7) of the Derivation of the Field Produced by a Shell of 
Magnetic Dipoles Secdon 

where the radius of Uie shell is r (The parameter r replaces the constant /? of Eq. (IIJ7) 
15 as the radius of the shell). For small r , the first term of Eq. (TV, 1) dominates and is 

used to approximately calculate the i -component of the magnetic field due to a sphere 

of 2 -oriented dipoles as follows: 

A sphere of dipoles is equivalent to (he integral over concentric shells of dipoles 

each of radius r where ^^r<R. Thus, tlie 5 -component of the magnetic field due to a 
20 sphere of dipoles is the integral of the field of the shells given by Eq. (IV.l). 



r^dr (IV.2) 



Let 



- - tanO; dr = pscc' OdO {IV 3) 

P 

and using the trigonometric identity 
25 l-l-tan'0-scc^9 (IV.4) 

Eq. (rv.2) becomes 

3 

Substit^ition in of 
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tanO=^ (IV.6) 

cose 



and 



sec0=-~ (IV.7) 

COS0 



into Eq. (TV^S) gives 



5 B, = (2z' - - y'}^^ ] -^sin' OcosOdO (rv.8) 

Let 

u = sin 9; du = cosOrfO (IV .9) 
Substitution of Eq. (TV.9) into Eq. (IV.8) followed by integration gives 

B, = (2^^-^-/)4f^77 (IV.IO) 

3 



10 Substitution of Eq. (1V.9) into Eq, (IV.IO) gives 

m 4;rsin^ 



3 

From Eq. (IV.3) 



r 



- = tanO; :sind - . ^ , (rV.12) 

fa 



Substitution of Eq. (rv.l2) into Eq. (IV.ll) gives 
)J!L. 



15 B,=(2z--.--/)^ ^ ; ^^,, |; (1V.13) 

i] 



Evaluation at the limits of the integral gives 
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(rvj5) 



Muldplying out the term in the denominator gives 



Sphere of Dipoles: 




{IV.16) 



10 



Eq. (IV, 16) demonstrates that the magnetic field of a sphere of dipoles is not 
equal to thai of a single dipolc at the origin. The shell radius, /?, appears in the 

15 denominator of the fSxsi tenn. The first temi is the magnetic field of a dipole at the origin 
only when the variable corresponding to the radius of the shell, R , is zero. The second 
icnn is an additional perturbation of the field of a dipole at the origin whose magnitude is 
a function of the radius of the sphere. The second tenn vanishes onJy when the radius of 
the sphere is zero. Thus, the dipole is a basis element for the reconstrucb'on of a NMR 

20 image. The resolution depends on the extent that die field of the sphere of dipoles diffei^ 
fipom that of a single dipole at the origin. The plot of (ho field of a sphere of dipoles of 
radius R and magnetic moment m = 10* Gcrn given by Eq. (IV J 6) as a function of 
radius R is given in FIGURE 6. The position of the center of the sphere relative to the 
detector is the point (0,0,10). The plot of the field of a sphere of dipoles of radius 

25 J? - 0.2 cm and magnetic moment w = 10^ Gcm^ given by Eq. (IV, 16) as a fimclion of 
the distance between the detector at the origin and the center of the sphere at the points 
(0,0, z) is shown in FIGURE 7. From FIGURES 6 and 7 it can be appreciated that the 
detector must be able to resolve three to four significant figures in order to reconstruct a 
map of 0.2 cm resolution with a field depth of 1 5 cm. 
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Fourier Transform of the System Function 
5 The system function, h(p^^,z), in cylindrical coordinates is 

[x +y +2 J [p +z ) 
The spacetimc Fourier transform in three dimensions in cylindrical coordinales, 
H(k^^^,k^), is given [9] as follows: 

fIiK>'^,K) = J J jA(p,*,z)exp^/2;r[A^pcos(a> - ^)+ k^zf^dpd^dz (V.2) 
1 0 With circular symmetry [9] 

The Fourier transfonn of the system function is given by the subslituiion of Eq. (V. I) 
into Eq. 0^-3), 

15 Consider the integral of Eq. (V.4) with respect to dp only. Factorization of h(p, ^, z) 
gives 

Consider the definite integral 

^ U'^z'r' 2T[u + l] 
20 Jind the modified Bessel function of the third kind relationship. 

3 

Tlie first factor of Eq. (V,5) is the same fcnm as Eq. (V.6) with v = 0; ii = , thus, 

(V.8) 
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where K,^,2[k^]= ^jnl^^] (EQ- (V.7)). The second factor of Eq. (V.5) can be made 
into Uie sanie form as Eq. {V.6) using Oie Bessel function of the first kind recurrence 
relationship 

J,.M + J,Jx]^^J,M (V.9) 
5 Consider tlie second factor of the integral of Eq. {V.5). 

Eq. (V.9) with o = 1 is 

JAxhJ2bh-/ti^] (v.ii) 

X 

10 Let 

x=k^p CVJ3) 
Substitution of Eq. CV.13) into Eq. (VJ2) is 

Substitution of Eq. (V.IO) into Eq. (V.I4) is 

(V.15) 

The first factor of the right hand side of Eq. {V.15) is the same forai as Eq. (V.6) with 
v= 1; « = ^,thus, 

(V.16) 

where /L.„j[fcpZ] = K„^[ki,2] (Eq.(V.7)). ITie second factor of theriglit liand side of Eq. 

3 

(V. 1 5) is the same form as Eq. (V.6) with v = 2; a = — , thus, 
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71Z /CL 



(V.17) 



Combining the parts of the integration with respect to dp of Eq. (V.4) by adding Eq. 
(V.8), Eq. (V.16), and Eq. (V.17) gives 



J- [ r[5/2] "'^ " ^ n5/2j ' 



(V.18) 



The modified Bessel function of the third kind formulae is 
Sjibstitiitionof Eq. (V.13) into Eq. (V. 19) with u = I is 

Substitution of Eq. (V .13) into Eq. (V.19) with u = 0 is 

Substitution of Eq. (V.20) and Eq. (V.21) into Eq. (V.18) is 



(V.19) 



(V.20) 



(V.21) 



£1 



r{5/2] [ 2V r[5/2} (2"')r[5/ 2jX2vJ J 



15 



CV.22) 



^lr(5/2] ' r{5/2]2 r(5/2j n5/2p " J 



20 Collecting terms gives 

f- rc'" f rrf3i 1 ,1 -tA*'.)" 

^-^5/231"' ■'12 X J 
With r[3] =2 and r[5/2]=3 /4n'". Eq. CV.24) is 



{V.23) 
(V.24) 

(V.25) 
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^^^3/2fc,e-''V•■''Vz (V.26) 



4Kk^^e'''e-''''Jz (V.28) 
4nk. fe''"'*'''"a: (V.29) 



P 



5 Integration of Eq. (V.29) with respect to dz gives 



•■•1:} 



(V.30) 



Multiplication of Eq. (V.31)by 



L-A+*pJ 



10 gives 



The system flinction (Eq. (V.l)) is an evai function; dins, tlie spacetinie Fourier 
. transfonn in tliree dimensions in cylindrical coordinates, Hik^X ) > given by taking the 
real pail ofEq. {VJ3) [8]. 

15 mA] = -^ . . . ■ ^''^ 

The spaceiime Fourier transfonn in three dimensions in C^esian coordinates, H{k^,k,), 



IS 



i^f^^ (V.35) 



where the relationship between Ihc wavcnumbers and the spatial Cartesian coordinates is 
20 as follows: 
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, 2k 2jr 



{V.38) 
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APPENDIX VI 



Derivation o(S = HF® U(k, ) frora Eq. (55). 



where the Faurier transform of u{z)~lforz ^0 andui^)^\forz<0 f8] is 

10 and where // is given by Eq. (52) and F is given by Eq, (49). The convolution integral 
of the second term of Eq. (VI. 1) is 

Collecting exponential lemns gives 

15 (vn.4) 

Expansion of the denominator of the second term gives 

(V1.5) 

20 The factorization of the denominator of the second tcnn, kI - 2A,ic, 4- -^k^, using the 
quadratic formula is 
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(VI.8) 



2 

= ' (VI.6) 

2 

Substitution of Eq. {VL6) into Eq. (YI.5) is 

5 (V1.7) 
The expansion of Eq. (VI.7) by the method ofpartial fractions is 

The /actors A, B, C of Eq. (Vl.g) are dctennined as fbllows: 

-2k,K, +k^+kiy % Jc. - A, -Jk^y C[k,^, -k, +Jk^y 1 (VL9) 
10 Ut K, = 0. Then, Eq. (VI.9) gives 

Thus, 

Let = A, + yAp. Then Eq. (VI.9) is 

Tlie tenn in iP vanishes. Substitution of Eq. (VL 1 1) into Eq. (VI. 12) gives 

^h^p + -'^^ -2ArJ -VKk.+kl +kiyC^jkJc^ -2kiy 1 (V1.13) 

where the lerm m A also vanishes. Solving for C gives 
20 Let If , = A. - jk^ . Then, Eq. (VI.9) is 



(VL12) 
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The term in C vanishes. Substitution of Eq. (VL 1 1) into Eq, (VI. 1 5) gives 

where the tenn in A aJso vanishes. Solving for B gives 
1 



5 = 



[-2A;-2yAi*,] 



Substitution of i4 , B , and C into the convolution integral (Eq. (VI,8) gives: 



»«^0 Jl,»-/y2 Ji.— /|/2 



The first convolution integral is of the form given by Mc Seibert f8]. 



— c:> sgnz^ ; sgnz^ = - 1 where z„ < 0 



10 Thus. 



4^ Ve»A>'- 



— I where <0 



A change of variabie in the second coiivohition integral given by letting 
K, = -A, -jk^ ; dK, = -dk, 

in 



15 



4< ] e*''-'' 



IS 



4nk' 



[-2A^-2Afc,].i 4-2*.] 



J .f 01.1 Where <0 



(Vr.l5) 



(VI.16) 



(VI.17) 



(VL18) 



(VI.19) 



(VI.20) 



(VI.21) 



(VI.22) 



(V1.23) 
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A change of variable in the third convolution integral given by letting 
5 K^=-k,+jk . dK^'^-dk, (yi.26) 



in 



IS 



10 



f ; — 2 ^ f (V1.27) 

— ^ ^ f ~ T-\-dK] where z. <0 (VI.28) 

[-2kl^2jk,k^]i j[-2k.] I '^J 



i: \ — dk, (V1.29) 

where z.<0 (VI.30) 



Combining Eqs. (VI.20), (yi.25), and (VI.30) gives the convolution of Eq. (VI.l). 

(VI.31) 
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APPENDIX Vn 

DcrivalioB of the Inverse Transfonn of Eq. (68) to Give Inverse Transfoim 1, Eq. (69). 

S 



L**».«o ,.'-t,n L % 



(VII.I) 



\2n ] ] e''i'-^J,[k^p},dk,e^'dk, ] 



whei« 2. <0 (VIL2) 
Consider Ihe definite integral: 

],-'jX„Y"^„?fM^:a><, (m3) 

0 [a +i J n:' 

In the case of the first integral of Eq. (Vn.2), the parameters of Eq. (Vn.S) are 

v = 0 

1 = k, 

,5 (vn.4) 

Thus, 
where 



(vn.6) 
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Consider the following relationship of tiie modified Bessei function of the tliird kind: 

,ri -jfC^O^r -;fl 

In the case tliat v = 0, Eq. (VTr.7) is 
5 Thus, 

(vn.!0) 

Willi a change of variable as follows: 
10 First 
Integral: Second Integral: 

' p p (vn.u) 

and an expansion of the complex exponential factor, Eq. (VII. 1 0) becomes 

= -] [cosA,2j+>sm[A:,z,]fcostt,z,]-jsm[*^,J,[A^jA, {VU.12) 

^ 0 

1 5 'Hie (lansform is given by Bateman [12]. 



-1 



Thus, 



{Vn.l3) 



0 l^fl ^-P J 

The Fourier trans fonn relationship between a function f[x] and its Fourier translbrm 
20 g[y] given by McC. Sicbcrt [8] is 
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^x]c^--~g[k] (V11.15) 
J dk 

Consider the following Inverse Fourier transfonn integral from Eq, (\^1].2) 

0 

Eq. (VU. 16) is equivalent to the following Inverse Fourier transfonn 
5 (^"•^''> 



FromEq.(Vn.l5)with 



Eq.CV1116)is 



(V11.19) 



1 0 Taking the imaginary part because Eq, (VE 1 9) has odd symmetry gives 

In the case of the third and fourth terms of Eq. (VII.2), consider tJie foUowing Fourier 
transform given by Gray [13]: 



0 

15 Thus, 



0 |?„ + p J 

The Inverse Fourier transform with ri^pecl to is given by McC. Siebert [8]. In the case 
of the first two terms of Eq. (yil.2), 

je^'^-'dk^^Slz] (VU.23) 

-09 

20 fa the case of the third and fourth tenns of Eq. (Vn.2), 
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Combining transforms from Eqs. (Vn.6), (Vn,20), and (V1123) in the case of the first 
and second tenns of Eq. (VII.2) and combining transforms from Eqs. (VH. 14) (\aL22) 
and (VII.24) in the case of the lliird and fourth terras of Eq. CV1I.2) gives: 



^^-0 Jt.-*I2 



y e 

Alt, 



{^f,'0 y,^-hn 



I 



27t8[z 



t.'*^>"r"(^+f>']"j 

1* 



where <0 



SubstiUition of Eq. (48) into Eq. (VIL25) gives 



(Vn.25) 

l^X.:=^0 Jt.=-^/2 L '^^ JJ 

where <0 

(VIL26) 

Taking the derivative given by tlie doublet function gives 

(<>,.»o ^...-.rj^a x„-4(2 L p JJ 

+V* w^z* 



15 



Inverse Transform 1 



(V!1.27) 



WO02yi695<; 



PCTAISOI/25954 



111 

REFERENCES 
The following citations are incoiporated by reference. 

5 1, MiJls, R., Magnetic Susceptibility Imaging (MSI), U.S. Patent No. 5.073,858 (1991). 

2. Reynolds, G. O., DeVelis^ J, B.. Parrent. G. B., Thompson, B.J., Tho New Ph y sical 

Optics Notebook, SPIE Optical Engineering Prcs5» (1990). 

3. Paiz, S., Cardiovasc Interven Radiol, (1986), 8:25, pp. 225-237. 

4. Jackson, J. D., Classical Electrodynaniics . Second Edition, John Wiley & Sons, New 
10 York, (1962), pp. 391-394. 

5. Ogawa, S.. Lee, T, Nayak, A. S.. Glynn, P., Magnetic Resonance in Medicine, Vol. 

14. (1990), pp. 68-78, 

6. San\'inski, R. E., "Superconducting Instruments", Cryogenics, Dec. 1977. pp. 671- 

679. 

15 7. Hounsfield, G. United Slates Patent No. 4,322,684, March 30. 1982. 

8. Sicbcrt, W., McC, Cii-cuiis. Signals, and Systems . The MTT Press. Cambridge, 

Massachusetts, (1986), p. 399.2. Siebcrt, W., McC, Circuits, S ig nals, and Syste ms.- 
The MTT Prcss, Cambridge, Massachusetts, (1986), pp. 415^16. 

9. Braceweli, R. R, The Fourier Transform and Its Applications , McGraw-Hill Book 
20 Ccrapajiy, New York, (1978), pp. 252-253, 

10. Sicbcrt, W.» McC. Circuits. Signals, and Systems. The MiT Press. Cambridge, 
Massachusetts, (1986), p. 574. 

1 1. Sieberl, W., McC, Circuits, Signals, and Svstcim . The MIT Press, Cambridge, 
Massachusetts, (1986), pp. 435-439 

25 12. Baleman, H., Tables of Integral Transforms , Vol. IH, McGraw-Hill, New York, 
(1954), p. 149, 

13. Gray, A., Mathews, G. B., A Treatise on Bcssel Functions and Their Applications to 
Physics, MacMillian and Co., Limited. London, (1952), p. p. 65. 

14. Mills. R.. Resonant Magnetic Susceptibility Imaging (ReMSI), U.S. Patent 
30 Application No. 09/191,454 filed November 12. 1998. 



wo 02/16956 



PCT/USOI/25954 



112 

CLAIMS: 

1 , Apparatus for providing a multidimensional image of inanimate or animate objects 
comprising: 

5 a magnetization source for producing a magnetic field to magndize a volume of 

an object to be imaged, the volume comprising a plurality of voxels; 

a radiation source for applying a radiation field to the object to be imaged when 
the object is disposed in the magnetic field to produce a secondary radiation field 
emanating from the object; 
10 an output signal detector for producing output signals in response to said 

secondary radiation at a plurality of spatial locations outside of tlie object as a 
function of time; 

a first processor for detennining a plurality of Fourier components, each having 
the same frequency, an intensity and a phase angle; 
! 5 a second processor for associating the Fourier components due to each voxel of 

the object to be imaged by phase to form a set corresponding to a spatial variation of 
intensity of the secondary radiation due to each voxel at said plurality of spatial 
locations outside of the object; 

a third processor for convening each said set of components into a voxel 
20 location; and 

an image processor for producing an image based on the secondary radiation 
intensity from each voxel and tlie location of each voxel. 

2. Tht apparatus according to claim 1, further comprising a display for displaying the 
25 image. 

3. The apparatus according to claim 2, characterized in that said display comprises at 
least one selected from the group consisting of printers, catliode ray tube displays^ 
liquid crystal displays, plasma screens, three dimensional modelers, laser monitors, 

30 projections monitors and holographic displays. 

4. The apparatus according to any one of the preceding claims, characterized in that one 
processor handles the processing fiinctions of at least two of said first processor, 
second processor, third processor and image processor 
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5, The apparatus according to any one of the preceding claims, characterized in tliat one 
processor handles the processing functions of at least three of said first processor, 
second processor, third processor and image processor 

5 6. The apparatus according to any one of the preceding claims, characterized in that the 
third processor comprises a reconstruction processor. 

7. The apparatus according to any one of the preceding claims, characterized in that the 
image is a magnetic resonance image. 

10 

8. The apparatus according to any one of the preceding claims, characterized in that the 
magnetization source comprises a superconducting magnet. 

9. The apparatus according to any one of the preceding claims, characterised in that the 
1 5 radiation source comprises a radio frequency generator and transmitter. 

10. The apparatus according to claim 9, characterized in that the radio frequency 
g^erator and transmitter comprises an antennae to excite protons in a magnetized 
volume. 

20 

1 1 . Tlie apparatus according to claim 10, characterized in that the ajUenna comprises two 
saddle coils driven in parallel to provide a rotating radio frequency field. 

12. The apparatus according to any one of the preceding claims, characterized in that the 
25 source of radiation produces a component of secondary radiation which is transverse 

to a direction of magnetization of a volume. 

13. The apparatus according to claim 12, characterized in that the secondary radiation 
which is transverse to the direction of the magnetization of the volume comprises 

30 radio frequency dipole radiation. 



14. The apparatus according to claim 13, characterized in tliat the radio frequency dipole 
radiation comprises near field radiation. 
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15. The apparatus according to any one of the preceding clmms, characterized in that the 
output signal delector comprises a plurality of detectors, 

16. The apparatus according to claim 15, characterized in that each of the detectors 
5 comprises an associated antenna. 

17. The J^paratus according to claim 2 6, characterized in that the plurality of detectors 
comprises a three dimensional array of antemias. 

10 IS. The apparatus according to claim 1 7, characterized in that each antenna of the three 
dimensional array of antennas responds to a secondary radio frequency field as a 
function of time. 

19. The apparatus according to claim IS, characterized in that the plurality of detectors 
] 5 are located at spatial locations outside of an object to be imaged. 

20- The apparatus according to any one of the preceding claims, characterized in that the 
output signal detector produces output signals in response to secondary radiation at a 
plurality of spatial locations as a function of time samples at a Nyquist rate in time 
over Ihc proton free induction decay. 

2 1 . The apparatus according to any one of the preceding claims, characterized in tliat tlic 
first processor iurthcr comprises an analog to digital converter to digitize signals of - 
secondary radiation. 

22. The apparatus according to any one of the preceding claims, characterized in that tlie 
first processor comprises a time Fourier transform processor to convert sccomlary 
signals from tlK signal output detector over time into Fourier components, each 
having the same frequency, an intensity and a phase angle. 

23. The apparatus according to any one of the preceding claims, characterized in that tlie 
plurality of Fourier components each having (he same frequency, an intensity and a 
phase angle are such that the phase of each component is unique. 
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24. The apparatus according to claim 23, characterized in that the unique phase of each 
component is provided by a combination of the angle 0 suspended between each 
spatial location at which the secondary field is detected and a radial vector, a vector 
from each voxel to each spatial location, and the angle ^ due to a separation distance 

5 r between each voxel and each spatial location given by a wavenumbcr of detected 
field k tinaes r . 

25. The apparatus according to any one of the preceding claims, characterized in that 
said first processor uses a first component having a phase angle, calculates the phase 

10 angle as a function of spatial position of a first spatial position relative to any other 
spatial location and identifies a component at each spatial location having the 
calculated phase angle. 

26. The apparatus according to claim 25, characterized in that a separate detector 
15 samples each separate spatial location. 

27. The apparatus according to claim 26, characterized in that each detector comprises a 
point detector. 

20 28. The apparatus according to any one of the preceding claims, characterized in lhat the 
source of radiation comprises at least one radio frequency generator selected from the 
group consisting of klystrons, backward wave oscillators, Guim diodes, and traveling 
wave tube amplifiers. 

25 29. The apparatus according to any one of the preceding claims, characterized in that the 
third processor uses a spatial variation of the secondary radiation over a sample 
volume comprising a plurality of spatial locations at which the radiation is detected 
to detcmiine a location of the voxels. 

30 30. Tlie apparatus according to claim 29, characterized in that the location of each voxel 
is determined independently by the third processor. 
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3 1 . The apparatus according to claim 2, wherein the image processor and display can 
rotate and display an image from a plurality of desired perspectives. 

32. The apparatus according to claim 2, characterized in that the image processor and 
5 display can display a reduced dimensional image. 

33. The apparatus according to any one of the preceding claims, characterized in that the 
apparatus provides a two dimensional or three dimensional nuclear magnetic 
resonance image of inanimate or animate objects. 

10 

34. The apparatus according to claim 33, characterized in that the output signal detector 
comprises a nuclear magnetic resonance detector for detecting a time dqicndcnt 
•Iransva^se radio frequency field component as a function of time of over a three 
dimensional sample volume outside of an object. 

15 

35. The apparatus according to claim 34, characterized in that secondary radiation 
emanating from a magnetized object comprises tlie transverse, time hamionic nuclear 
magnetic resonance magnetic dipole radiation. 

20 36. The apparatus according to claim 34, characterized in that the time dependent 

transverse radio frequency field component comprises a harmonic oscillation of each 
radio frequency dipole that is equivalent to a dipole rotating in a transverse plane. 

37. The apparatus according to claim 34, characterized in that the apparatus can generate 
25 an image from a three-dimensional map of a transverse resonant radio frequency 
magnetic flux external to an object and said nuclear magnetic i^onance detector 
com{»rises a plurality of detectors at spatial locations, and wherein 

a nuclear magnetic resonance signal of each voxel at a given detector gives rise to 
a Fourier component with a unique phase angle relative to a Fourier 
30 component of other voxels of the object at that detector; 

a set of Fourier components that correspond to the nuclear magnetic 
resonance signal of a given voxel over the detectors is detennined; and 

an intensity variation of a transverse radio frequency field over space is 
used to determine a coordinate location of each voxel. 



wo 02/16956 



PCT/USOl/25954 



in 

38. The apparatus according to claim 37, characterized in that the each detector is 
selectively responsive to a component of the secondary field. 

5 39. The apparatus according to claim 38, characterized in that a transverse radio 

frequency field is a near field which is a dipole that serves as a basis element to form 
a unique reoonstniction. 

40. The apparatus according to claim 39, characterized in that the each detector is 
10 responsive to a component of radio frequency in a transverse plane, 

41. The apparatus according to claim 40, characteri2ed in that each detector is 
responsive to a component of the radio frequency dipole field parallel to the radio 
frequency dipole. 

^ 15 

42. The apparatus according to claim 41, characterized in that each detector is 
responsive to the component of the radio frequency dipole field parallel to the radio 
frequency dipole as a function of time and detects a maximum intensity of the signal 
during each time hannonic cycle, 

20 

43. The apparatus according to any one of the preceding claims, characterized in that the 
third processor determines a location of a voxel from signals recorded of a transverse 
dipole radio frequency field. 

25 44. The apparatus according to claim 43, characterized in tliat a geometric system 
function corresponding to a dipole determines a spatial intensity variations of tlic 
radio frequency field. 

45. 'I1ie apparatus according to claim 44, characterized in that the system function is a 
30 band-pass when = k^, the frequency of spatial variation of the recorded secondary 
radiation in the rho direction is equal or comparable to that in the z direction of the 
three dimensional space comprising said plurality of spatial locations. 



wo 02/16956 PCTAjSOl/25954 

118 

46. The apparatus according to claim 45, characterized in that each voxel is 
reconstructed independently in pgffallcl with other voxels using each unique set of 
components associated by phase such that the scan time is no greater than the nuclear 
free induction decay time. 

5 

47. The apparatus according to claim 46, characterized in that all voxels comprising an 
image arc reconstructed in parallel. 

48. The ^paratus according to claim 46, characterized in that the apparatus can generate 
10 higli resolution, three-dimensional, reaMime anatomical images and images based on 

physiological parameters with little or no deterioration from motion artifact, 

49. Ttkc apparatus according to claim 41 » characterized in that at each point in time, each 
radio frequency dipole rotating in the transverse plane is directed at an angle 6 

1 5 relative to the direction of detection of each detector. 

50. The E^paratus according to claim 49, characterized in that the phase angle 6 of the 
radio frequency dipole relative to the direction of detection axis of the detector gives 
rise to a phase angle tenn u"'^. 

51. The apparatus according to claim 50, characterized in that at each point in time, each 
radio frequency dipole rotating in the transverse plane with a separation distance r 
between each voxel and cacli spatial location at a detector gives rise to an angle (p 
given by a wavcnuinbcr of the detected secondary radiation k times r directed at the 
angle 0 relative to a direction of detection of each detector. 

52. The apparatus according to claim 51, characterized in tliat the phase angle ^ of the 
radio frequency dipole relative to a direction of detection axis of a detector gives rise 
to a phase angle term c**'* directed at angle 6 relative to the direction of detection of 
each detector and wherein the angle 0 gives rise to a phase angle term e"^, 

53. The apparatus according to claim 52, characterized in thai a sum of the phase angles, 
kr and 6 , is unique for each voxel at each detector. 
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54. The apparatus according to claim 37, characterized in that the unique phase of each 
component is provided by a combination of 

angle 0 suspended between each spatial location at which a secondary nuclear 
5 magnetic resonance field is detected and a radial vector from each voxel to each 
spatial location^ and 

angle (fi due to a separation distance r between each voKel and each spatial 
location given by the wavenumbcr of delected nuclear magnetic resonance field k 
times r . 

10 

55. The apparatus according to claim 37, further comprising a processor, characterized 
in that a set of Fourier components that correspond to a nuclear magnetic resonance 
signal of a given voxel over detectors is detennined by the processor using a position 
of a first detector relative to a different detector to calculated the phase angle of the 

1 5 different detector relative to the first detector. 

56. Tlie apparatus according to claim 55, characterized in that the processor detennincs 
the set of Fourier compontaits that correspond to the nuclear magnetic resonance 
signal of a given voxel over the detectors by using the position of each first detector 
relative to a different detector to calculated the phase angle of the different detector 
relative to the first detectors; and the process is repeated over all of the detectors to 
give the set of intensities of the nuclear magnetic resonance signal over the plurality 
of detectors due to each voxel. 

57. The apparatus according to claim 56, characterized in that the set of intensities of the 
nuclear magnetic resonance signal over the phmility of detectors due to each voxel is 
processed to determine a coordinate location of each voxel. 

58. The apparatus according to any one of the preceding claims, further comprises coils 
for providing 

, a steady primary field, 
G,, a field gradient in an x-axis direction, 
Gy, a field gradient in a y-axis direction, 
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a radio frequency field, and 

C7y, a field gradient in a z-axis direction. 



59. The apparatus according to claim 58^ associated amplifiers for driving the coils and 
5 associated control circuits for controlling the coils which provide T, and pulse 

sequences to alter an intensity of the secondary radio frequency field as a iFbaction of 
time to provide corresponding signals to be processed. 

60. The apparatus according to claim 9, characterized in that nuclei present in object are 
10 subjected to an additional rotating radio frequency field H, , which is synchronous 

with precession and their magnetic moments, and magnetization M is caused to 
process about H, and rotate away from tlie primary field H^^ by an angle ^ in a 
coordinate frame which rotates at the Larmor frequency, the precession about II, 
continues as long as ll^ exists, and a j&nal value of which depends on the strength 
15 of Hp which determines the precession rate, aiid the time for which it is turned on is 
such that iff = 90® that a dipole is perpendicular to a detector array that is responsive 
to the nuclear magnetic resonance radiation at a plurality of ^atial locations external 
to the object. 

20 61. Tlie apparatus according fo any one of the preceding claims, characterized in that the 
magnetizing source is constructed to provide a magnetic field which permeates the 
object to be imaged and is confined to a volume to be imaged, and the confmed 
magnetic field limits the source of signal to the volume of interest, characterized in 
that the volume to be reconstructed is limited to the magnetized volume which sets a 

25 limit to required computations and reduces end efTccts of signal originating outside 
of edges of the output signal detector. 

62. T1ie apparatus according to claim 61 ^ characterized in that the apparatus is 
constructed to provide a magnetizing field which limits an imaged volume of aii 
30 object by providing a range of Larmor frequencies and data is rejected at a threshold 
Larmor frequency. 
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63. The apparatus according to any one of the preceding claims, characterized in that the 
appaiatus is constnicted to provide a magnetizing field which limits an imaged 
volume of an object by providing a range of Larmor frequencies and data is rejected 
at a threshold Larmor frequency. 

5 

64. The apparatus according to any one of the preceding claims, characterized in that the 
output signal detector comprises detector elements disposed on a first axis and 
having corresponding output sigaals. 

10 65. The apparatus according to claim 64, characterized in that the output signal detector 
includes detector elements disposed on a second axis and having a corresponding 
output signal, and the detectors on the second axis are responsive to the 
corresponding oriented field components of tlic secondary radiation. 

15 66. The apparatus according to any one of the preceding claims, characterized in that 
tlic third processor comprises at least one reconstruction processor selected from the 
group consisting of a detector processcK* for processing the detector output signal 
with an inverse matrix, a reiterative processor, a time Fourier transform processor, a 
sine function filter processor for detecting the detector output signals, a space and 

20 time Fourier transfonn processor, and a compensating processor for compensating 
for a distance from the detector along a detector responsive axis. 

67. Tlie apparatus according to any one of the preceding claims, characterised in that the 
signal output detector ftirther comprises a line or plane of detectors aad a drive 

25 mechanism that moves the detectors to sample an external secondary radiation over a 
three-dimensional volume. 

68. The apparatus according to claim 67, characterized in that the secondary radiation is 
a radio frequency magnetic field due to nuclear magnetic resonance of nuclei in a 

30 magnetic field produced by the magnetization source and excited by a radiation field 
from the source of radiation: 



69. 



The apparatus according to claim 68, characterized in that the nuclear magnetic 
resonance active nuclei are protons. 
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70. The apparatus according to clsum 19, characterized in that each antenna of tlie llircc 
dimensional array of antennas are impedance matched. 

5 71 . The apparatus according to claim 70, characterized in tiiBi each antenna of the three 
dimensional array of antennas further comprises supcrcotiducting reactance elements 
which minimize resistive losses. 

72. The a5>paratus according to claim 19, characterized in that cross talk between 

10 antennas is ameliorated or eliminated by time multiplexing the signal detection over 
the array of antennas. 

73. The apparatus accordbg to claim 19» characterized in that the signal output detector 
comprises a synchronous detector wherein tlic radio frequency field is sampled 

15 synchronously so that a phase at any given detector may be related to that at any 
other detector, 

74. The apparatus according to claim 19, characterized in that the signal output detector 
comprises a synchronous detector wherein the radio frequency field is sampled at 

20 known times so that a phase at any given detector may be related (o thai at any othei- 
detector. 

75. The apparatus according to claim 1 9, characterized in that each antenna of the three 
dimensional array of antennas responds selectively to a component of tlie secondary 

25 radio frequency field at a spatial location outside of the body as a function of time. 

76. The apparatus according to claim 1 9, characterized in tliat each antoma of the tlirec 
dimensional array of antennas responds selectively to a parallel component of the 
secondary radio frequency field at a spatial location outside of an object as a function 

30 of time. 

77. The apparatus according to claim 19, characterized in (hat each anteima of the three 
dimensional airay of antennas responds selectively to a perpendicular component of 
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the secondary radio frequency field at a spatial location outside of an object as a 
function of time. 



78. Apparatus for providing a multidimensional image of inanimate or animate objects 
5 comprising: 

means for producing a magnetic field to magnetize a volume of an object to be 
imaged, the volume comprising a plurality of voxels; 

means for applying a radiation field to the object to be imaged when the object is 
disposed in the magnetic field to produce a secondary radiation tield emanating from 
ID the object; 

means for producing output signals in response to said secondary radiation ai a 
plurality of spatial locations outside of the object as a function of time; 

first processor means for determining a plurality of Fourier components, each 
having the same ia-equency, an intensity and a phase angle; 
15 second processor means for associating the Fourier components due to each voxel 

of the object to be imaged by phase to form a set corresponding to a spatial variation 
of intensity of the secondary radiation due to each, voxel at said plurality of spatial 
locations outside of the object; 

third processor means for converting each said set of components into a voxel 
20 location; and 

image processor means for producing an image based on the secondary radiation 
intensity from each voxel and the location of each voxel. 

79. An apparatus for providing a two dimensional or tburec-dimensional nuclear 
25 magnetic resonance image of an object comprising: 

means for generating an unage iGrom a three-dimensional map of a transverse 
resonant radio frequency ma^etic flux externa! to an object having a volume to be 
imaged, the volume comprising a plurality of voxels; and 

processor means for determining the location of each voxel which is a source of 
30 the radio frequency field through spatial variations of an intensity of the nuclear 
magnetic resonance field of a set of components associated by phase. 

80. An apparatus for providing a two dimensional or tlireo-dimensional nuclear 
magnetic resonance image of an object comprising: 
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a raagaetizing source for magnetizing a volume an object in a first direction, the 
volume comprising a plurality of voxels; 

means for causing nuclear magnetic resonance magnetization (o rotate into a 
plane transverse to the first direction; 

an output signal detector for detecting the nuclear magnetic resonance field 
rotated into a transverse plane external to the object at a plurality of spatial 
locations; 

a first processor for determining a set of components of spacial variation of the 
nuclear magnetic resonance field over a sample volume by association of the 
components by phase; 

. a second processor for detcnnining a location of each voxel which is a source of 
the radio frequency Held through each set of components representative of spatial 
variations of an intensity of the nuclear magnetic resonance field over ihe sample 
volume, and 

an image processor for displaying tlie voxel locations representative of 
anatomical images of the object based on selected physiological 

8L An apparatus for providing a two dimensional or three-dimensional nuclear 
magnetic resonance image of an object comprising: 

means for magnetizing a volume an object in a first direction, (he volume 
• comprising a plurality of voxels; 

means for causing luiclcar magnetic resonance magnetization to rotate into a 
plane transverse to the first direction; 

means for detecting the nuclear magnetic resonance field rotated into a transverse 
plane external (o the object over a sample volume comprising a plurality of spatial 
locations; 

processor means for detcnnining a set of components of spatial vaiiaiion of the 
nuclear magnetic resonance field over a sample volume by association of the 
components by phase; 

processor nieajis for detcmuning a location of each voxel which is a source of the 
radio frequency field through each set of components representative of spatial 
variations of an ijQtensity of the nuclear magnetic resonance field over Che sample 
volume, and 
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processor means for displaying the voxel locations representative of anatomical 
images of the object based on selected pliysiological pardmeters. 

82. An apparatus for providing a two dimensional or three-dimensional nuclear 
5 magnetic resonance image of an object comprising: 

means for generating an image from a three-dimensional map of a transverse 
resonant radio frequency magnetic flux external to an object comprising a volume 
to be imaged, the volume comprising a plurality of voxels; and 
processor means for determining locations of voxels which is a source of a radio 
10 frequency field through spatial variations of an intensity of a nuclear magnetic 
resonance field of a set of components associated by phase. 

83 . A mclliod for providing a multidimensional image of inanimate and animate objects 
comprising the steps of: 

1 5 placing an object having a voliunc in a magnetic field, Ihe volume comprising a 

plurality of voxels; 

applying a first source of radiation to the object; 

recording a secondary source of radiation emanating fi-om the object at a plurality 
of spatial locations; 

20 forming components of the spatial variation of the secondary source of radiation 

external to the object due to voxels of the object; 

determining the location of voxels from Uie intensity variation of the components; 
. generating and displaying the image from the locations and the intensity of the 
secondary radiation of each voxel. 

25 

84. The method according to claim 83, characterized in that the secondary recorded 
radiation is a three-dimensional transverse resonant radio firequency magnetic flux 
external to the object recorded as a function of time. 

30 85. The method according to claim 84, characterised in that the three-dimensional 

transverse resonant radio frequency magnetic flux external to tlie object is nuclear 
magnetic resonance radiation. 



wo 02/16956 PCT/XJSOl/25954 

126 

86. The method according to any one of claims 84 or 85, further comprising recording 
the secondary source of radiation emanating &om the object at a plurality of spatiaJ 
locaticms using a plurality ofdctectors at each spatial location. 

5 87. The method according to any one of claims 83-86, characterized in that the recorded 
signals from each voxel contains components each having a unique phase angle 
relative to other components from different voxels. 

88. The method according to any one of claims 83-87, characterized in that a time 
10 dependent signal at each detector may be transformed into a series of components 

having intensity and phase data. 

89. The method according to any one of claims 83-88, characterized in that each set of 
components of the nuclear magnetic resonance signal over a sample volmne due to a 

1 5 given voxel is determined frofn phase data and detector positions. 

90. The mctiiod according to any one of claims 83-89, characterized in that an intwisity 
variation of a transverse ladio frequency field over the sample volimic is used to 
detcnninc the coordinate location of each voxel. 

91. The method according to any one of claims 83-90, characterized in that the radio 
frequency field is a near field wliich is a dipole that serves as a basis element to fonn 
a unique reconstruciion. 

92. The method according to any one of claims 83-91, characteri;&ed in that each volume 
element is reconstructed independently in parallel with all other volume elements 
such that a scan time is no greater than a nuclear free induction decay time. 

93. The method according to any one of claims 83-92, characterized in that said object is 
at least one selected from the group consisting of a human body, an animal, and 
tissue. 



94. A metliod for providing a multidimensional image of inanimate and animate objects 
comprising die steps of: 
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magnetizing a volume of an object to be imaged, t!ie volume comprising a 
plurality of voxels; 

applying a radiation field to the object to be imaged when disposed in the 
magnetic field to produce a secondary radiation field emanating &om the object; 
5 producing output signals from a detector in response to said secondary radiation 

at a plurality of spatial locations outside of the object as a function of time; 

determining a plurality of Fourier components each having the same frequency, 
an intensity and a phase angle; 

associating the Fourier components due to each voxel of the object to be imaged 
10 by phase to form a set corresponding to tlic spatial variation of intensity of the 

secondary radiation due to each voxel at the plurality of spatial locations outside of 
the object; 

converting each set of components into a voxel location; and 
producing an image based on an intensity of the secondary radiation from each 
1 5 voxel, the location of each voxel, and plotting the superposition of results for each 
voxel. 

95. A method for providing a two dimensional or three-dimensional nuclear ma^ietic 
resonance image of an object comprising the steps of: 

20 generating an image from a three-dimensional map of a u-ansverse resonant radio 

frequency magnetic flux external to a vohune of an object being scanned^ the 
volume comprising a pliuality of voxels, characterized in that a nuclear magnetic 
resonance signal of each voxel at aj\y given detector gives rise to a Fourier 
component with a unique phase angle relative to the Fourier component of any other 
voxel of the object at a detector; 

a set of Fourier components that correspond to a nuclear magnetic resonance 
signal of a given voxel over detectors is detennined; 

an intensity variation of the transverse radio frequency field over space is used to 
determine a coordinate location of each voxel; and 

a superposition of results for each voxel is plotted to provide a total image. 

96. A method of providing a multidimensional image of inanimate and animate objects 
comprising the steps of: 
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aligning magnetic moments of nuclei of a volume of an object to be imaged using 
a primary field, the volume comprising a plurality of voxels; 

further aligning the magnetic moments by a radio frequency pulse or series of 
pulses; 

5 recording free induction decay signals; 

Fourier transforming time dependent nuclear magnetic resonance signals to give 
intensity and phase of each Fourier component, characterized in that nuclear 
magnetic resonance signal of each voxel at any given detector gives rise to a Fourier 
component witli a unique phase angle relative to the Fourier component of any 
1 0 other voxel at that detector; 

determining a matrix of Fourier components that correspond to the nuclear 
magnetic resonance signal of a given voxel over tiic detectors; 
determining matrices for a plurality of voxels; 

determining a coordinate location of each voxel from measurements of spatial 
1 5 variations of a transverse radio frequ^cy field of each conesponding matrix; and 
processing and displaying the position of each voxel with a representation of the 
intensity of the nuclear magnetic resonance signal emitted by each voxel. 

97. A mell^od according to any one of claims 83-96, characterized in that a strength and 
20 duration of rotating radio frequency field that is resonant with the protons of the 

magnetized volume and is oriented perpendicularly to a direction of the magnetizing 
field is applied such that final precession angle of magnetization is 90 °, ^ 90^ , 

such that a radio frequency dipolc is transverse to a pritnary magnetizing field and 
perpendicular to a radio frequency magnetic field detector. 

98. The method according to any of claim 83-97, characterized in that nuclear magnetic 
resonance pulse sequences which provide signals for a 7J or 7^ image are applied, 

99. Tlte method according to claim 98, characterized in that the nuclear magnetic 
resonance pulse sequences which provide the signals for a 7J or image comprises 
applying a 90 ^ pulse followed by a series of 180 ** pulses. 
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100. The method according to claim 99, charactOTzed in that the step of applying a 90 
^ pulse followed by a series of 1 80 pulses comprises the Carr-Purce!l-Meiboom- 
GiU sequence. 

5 101, The method according to any one of claims 96-100, characterized in that the step 
of determining the matrix of Fourier components that correspond to the nuclear 
magnetic resonance signal of a given voxel over the detectors is achieved by using a 
first component having a phase angle and calculating the phase angle as a function of 
spatial position of a first detector relative to any other detector and identiiying a 
1 0 component at each detector having the calculated phase angle. 

102. The method according to any one of claims 96-100, characterized in that each 
matrix of components associated by phase comprises tiie intensity variation over the 
sample volume of the radio frequency field of bulk magnetization M of each voxel. 



1 03. The method according to claim 96, fiirtlier comprising reconstructing a 
multidimensional image representation according to spatial variations of the detected 
radio frequency magnetic field using at least one of a matrix inversion algorithm, a 
reiterative algorithm or Fourier Iransfomi algorithm. 

104. The method according to claim 103. characterized in that a spatial Fourier 
transform algorithm is used. 



105. The method according to claim 104, characterized in that the spatial Fourier 
25 transform algorithm is perfonned on each set of components over tJie detector anay 

to map each bulk magnetization M corresponding to a voxel to a spatial location of 
tlie imaged object 

106. The method according to claims 1 04 or 1 05, characterized in that the spatial 
30 Fourier transfonn comprises the steps of : 

recording a nuclear magnetic resonance signal at discrete points in the sampling 
space at position x, y, z as a function of time; 

designating each point, x, y, z, as an element in matrix A ; 
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FouriCT transforming time dqpendcnt signals to give an intensity and phase of 
each component wherein the nuclear magnetic resonance signal of each voxel at any 
given detector gives rise to a Fourier component with a unique phase angle relative 
to the Fourier component of any other voxel of the object at that detector; 
5 determining a matrix of Fourier components that cotrespond to the nuclear 

magnetic resonance signal of a given voxel over the detectors; 

determining matrices A„ for all of the voxels; 

delemiining a coordinate location of each voxel using the measurements of the 
spatial variations of the transverse radio frequency field of a given matrix; and 
10 plotting a superposition of results for each voxel to provide a total image. 

107. The method according to claim 106, characterized in that the step of determining 
the coordinate location of each voxel using the measurements of the spatial variations 
of the tiansvme radio frequency iield of a given matrix comprises the steps of 

1 5 a.) discrete Fourier transforming each matrix A„ to obtain a matrix ; 

b. ) multiplying each element of each matrix B„ by a corresponding inverse value 
or redprocal value of the Fourier transform of a system function which gives the 
spatial variation of the radio frequency field emitted by each voxel as a function of 
distance from the voxel to each detector, evaluated at the same frequency as the 

20 element of the matrix to form matrix C„ ; 

c. ) generating matrix D„ by taking the discrete inverse Fourier transform of 
matrix Q,; 

d. ) multiplying each element of each matrix by the distance squared along the 
z-axis to which the element corresponds to generate the position of the bulk 
magnetization M of voxel n which corrects the limitation of the sample volume to 
2>0; 

c.) repeating steps a«d for all matrices A^; and 

f ) plotting superposition of results for each voxel to provide a total image. 

108. The method according to claims 106 or 107, characterized in that at each 
sampled external spatial location, x, y, and z in Cartesian coordinates, a magnitude 
B' of the transverse radio frequency magnetic flux emitted by the corresponding z- 
axis aligned magnetic moment of each voxel is given by 
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where the magnetic momeat corresponds to the voxel's bulk magnetization M . 

109 . The method according to any one o f claims 106-108, fiirther compri si ng the step 
5 ofdetermining a matrix of Fourier components that correspond to the nuclear 

magnetic resonance signal of a given voxel over the detectors by using a first 
component having a phase angle and calculating the phase angle as a function of 
spatial position of a first detector relative to any other detector and identifying a 
compoacnt at each detector having (he calculated phase angle, 

10 

1 10. The method according to claim 109, characterized in that each matrix of 
components associated by phase comprises the intensity variation over the sample 
volume of the radio frequency field of the bulk magnetization M of each voxel. 

15 111. The method according to claims 106 or 1 07, characterized in that a S3rstem 
function as a Amotion of the detector location relative to each voxel in Cartesian 
coordinates, x, y, z, is by 



1 12. The method according to any one of claims 106- 111, further comprising the step 
of correcting a point spread of each rcconstnictcd voxel by assigning one voxel 
above a certain threshold with the bulk magnetization M and the other voxels are 
assigned a zero value. 

1 1 3. TIic method according to any one of claims 1 06- 1 1 2, further comprising the step 
of correctmg a finite size of each detector by multiplying each numeric value of 
matrix B„ by the conesponding matrix element of the inverse or reciprocal of the 
Fourier transform of a waveform representative of a detector geometry. 

1 14. The method according to claim 106-1 13, characterized in that each matrix 
element of the inverse or reciprocal of tlie Fourier transform of a waveform 
representative of the detector geometry is given by: 
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115. The method according to anyone of claims 106, 107, or 1 1 1-114, characterized 
in that each volume element is reconstructed independently in parallel with all other 

• 5 volume elements such that a scan time is no greater than a nuclear free induction 
decay time. 

1 16. The method according to any one of claims 106, 107, or 1 11-1 15 wherein the 
Fourier transform of the system fimctioii is given by: 

117. The method according to any one of claims 1 06- 1 1 6, characterized in that each 
matrix clement for tlie step of multiplying each element of each matrix by the 

. distance squared along the x-axis to which the element coiresponds to generate the 
15 position of the bulk magnetization M of voxel n which corrects the limitation of the 
sample volume to z > 0 is given by: 




20 118. The method according to claim 1 03, characterized in that the step of 

reconstructing the multidimensional image uses a matrix inversion reconstruction 
algorithm. 

1 19, The metfiod according to claim 118, characterized in that the matrix inversion 
25 reconstruction algorithm comprises the steps of : 

using the geometric system function to determine the spatial intensity variation of 
the transverse radio frequency field over a detector array to form a corresponding 
matrix; 

inverting the corresponding matrix; and 
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multiplying the signal over the detector array by the inverted matrix to give voxeJ 
sources. 

120. The metliod according to claim 1 19, further comprising the steps of 
5 a.) determinijig a coefficient for each voxel mathematically or by calibration 

which when multiplied by tiie bulk magnetization M of each voxel is that voxel's 
contribution to tlie signal at a given detector with the correspoading imique phase at 
each detector; 

b.) repeating step a for every detector; 
1 0 c.) using the coefficients to determine a matrix wliich when multiplied by a 

colunm vector of the bulk magnetization M values of the voxels, gives the signals at 
the detectors; 

d. ) inverting the matrix and storing the result in memory; 

e. ) recording signals as a function of time over the detector array; 

1 5 f.) Fourier transfomiing the signals as a function of time to give the Fourier 

components each having an amplitude and a unique phase; and 

g.) multjpl>'ing the components by the inverse matrix to generate a bulk 
magnetization M map which is displayed. 

20 121. The method according to claim 1 20, characterized in that the coefncient for each 
voxel determined mathematically is given by 

where is a detector signal at Csfftesian coordinate x, y, z, is a constant 

characteristic of each detector, Z^, and ly are the total lengths of the detector array 
25 in the x, y, and z directions. ^, Aj» and ky arc spacing of elements of the detector 
array in the x, y, and z directions, and #f„ n^, andn^ are integer indices of the 
elem^ts of the detector array in the x, y, and z directions. 



30 



122. 



Tlic mctliod according to claim 1 1 8, further comprising the steps of 
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detcmuniiig a matrix of Fourier components that correspond to the nuclear 
magnetic resonance signal of a given voxel over the detectors; 

determining corresponding matrices for a plurality of voxels; 

determining a coordinate location of each voxel from measureiiicnts of spatial 
5 variations of the transverse radio frequency field of a given matrix M^herein each 

matrix of components associated by phase comprises the intensity variation over the 
sample volume of the radio frequency field of tlie bulk magnetization M of each 
voxel; and 

multiplying each matrix of components by the inverse matrix of system fimction 
10 to generate tlie bulk magnetization M map. 

123. The method according to claim 122, further comprising the steps of 

a.) correcting a point spread of a reconstructed voxel by assigning one voxel 
above a certain tlireshold with the bulk magnetization M ; 
1 5 b.) assigning other voxels with a zero value; 

c. ) repeating steps a-b for other voxels; and 

d. ) plotting and di^laying the superposition of the separate msp$ for each voxel 
which comprises the total bulk magnetization map as a nuclear magnetic resonance 
image. 

20 

124. The method according to claim 103, characterized in that said multidimensional 
image is reconstructed using a reiterative reconstruction algoritlmi. 

125. The method according to claim 1 24, characterized in that the reiterative 
25 reconstruction algorithm comprises the steps of 

a. ) using the geometric system function to dctemiine a system of linear equations 
which gives the intensity, spatial variation, and phase of radio frequency field over a 
sample volume wherein the system of linear equations gives the signal and phase 
Ccom each voxel at each detector based on the bulk magnetization M value of each 

30 voxel and the position of the voxel relative to the dctcctor> 

b. ) Foim'er transforming the signal as a function of time to give Fourier 
components each having an amplitude and a unique phase; 

c. ) determining weighting coefficients based on tlie linear equations; 

d. ) determining coefBcicnts for each voxel mathematically or by calibration; 
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e. ) estimating bulk magnetization M for each voxel; 

f. ) calculating signals at each phase at each detector; 

gO comparing a result of multiplying the bulk magnetization M value of each 
voxel times its weighting coefficient and its calculated phase at a given detector to a 
5 measured signal and phase; 

h. ) making a correction to M of each voxel to^ve a second, or recomputed, 
estimate for M of each voxel; 

i. ) computing the signal value from this second estimate and makiiig corrections 
according to steps f-h; 

10 j.) repeating steps f-i until the correction for each reiteration approaches a 

predefined limit which serves to indicate that the reconstruction is within reasonable 
limits of error; 

k,) repeating steps c-j of the procedure for other detectors; 
and 

1 5 I.) plotting and displaying a bulk magnetization map. 

1 26. The method according to claim 125, wherein the coefficient for each voxel 
determined mathematically is given by 

20 where y[x,y,z] is tlie detector signal at Cartesian coordinate x, y, z, is a constant 
characteristic of each detector, /„ i^, and are tlie total lengths of the detector array 
in the x, y, and z directions, k^, k^, and are the spacing of the elements of the 
detector array in the x, y, and z directions, and n,, ri^, and are the integer indices of 
the elements of the detector array in tlie x. y, and z directions. 

25 

127. A method of providing a multidimensional image of inanimate and animate 
objects comprising tiie steps of : 

a, ) recording secondary radiation signals corresponding to a transverse radio 
frequency magnetic field in response to a first radiation applied to a magnetized 

30 object over a sample volume comprising a plurality of voxels; 

b. ) fomiing a matrix of signals at a plurality of spatial positions each 
corresponding to a detector element of a detector array; 
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c. ) transfonning the matrix numeric values according to a Fourier transfonn to 
give an intensity and phase of each component wherein a nuclear magnetic 
resonance signal of each voxel at any given detector gives rise to a Fourier 
component with a unique phase angle relative to the Fourier component of any other 

5 voxel of the object at that detector, 

d. ) determining a matrix of Fourier components that correspond to the nuclear 
magnetic resonance signal of a given voxel over the detectors; 

e. ) multiplying each numeric value of the tnuisform array by a value which Is an 
inverse or reciprocal of the Fourier transform of a waveform corresponding to aji 

1 0 operation perfonned in step b; 

f. ) multiplying each numeric value of the discrete spectrum matrix by a value 
corresponding to an inverse Fourier transfonn of a function rqjrescntative of the 
steps of exciting and detecting and providing a system corrected matrix; 

g. ) inverse-transfonning the system corrected matrix according to a 

15 multidimensional discrete inverse Fourier transform providing a voxel location; 

h. ) correcting each element in the transformed matrix for any asymmetries in the 
dimensions of the sample volume; 

i. ) rei>caling steps c-h for other matrices, each corresponding to a distinct voxel; 
and 

20 j.) superimposing the results from steps a-h, plotting and displaying the image. 



128. The method according to claim 127, characterized in tliat at each sampled 

external spatial location, x, y, and z in Cartesian coordinates, the magnitude B of the 
transverse radio frequency magnetic jflux emitted by the corresponding z-axis aligned 
25 magnetic moment of each voxel is given by 

where the magnetic moment corresponds to the voxel's bulk magnetization M . 



30 !29. The method according to claims 127 or 128, characterized in that the step of 
determining the matrix of Fourier components that correspond to the nuclear 
jnagnetic resonance signal of a given voxel over the detectors is achieved by using a 
lirst component having a phase angle and calculatmg the phase angle as a function of 
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Spatial position of a first detector relative to any other detector and idenrtfying the 
component at each detector having the calculated phase angle. 

130. The method according to claim 129, characterized in that each matrix of 

5 components associated by phase comprises an intensity variation over the sample 
volume of the radio frequency field of the bulk magnetization M of each voxel. 

131. The metiiod according to any one of claims 1 27- 1 3 0, characterized in that a 
system fimclion as a function of (he detector location relative to each voxel in 

10 Cartesian coordinates, x, y, 2> is given by 



(X 4-Z ) 

132. TItc method according to any one of claims 1 27- 1 3 1 » further comprising the step 
of correcting a point spread of each reconstructed voxel by assigning a voxel above a 
15 certain threshold with the bulk magnetization M and other voxels arc assigned a zero 
value. 



1 33. The method according to any one of claims 1 27-132, further comprising stop e of 
correcting a finite size of each detector by multiplying each numeric value of tlie 

20 matrix formed in stej) d by the corresponding matrix element of the inverse of the 
Fourier transfomi of a waveform representative of a detector geometry. 

1 34. The method according to claim 1 33, characterized in that each matrix element of 
an inverse of the Fourier transform of a waveform representative of the detector 

25 geometry is given by formula 

1 

sinc(A„A^,AJ 



30 



135. The method according to any one of claims 127-135, characterized in tlial each 
volume element is reconstructed indq)endently in parallel with all other volume 
elements such that a scan time is no greater than a nuclear free induction decay time. 
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136. The method according to any one of claims 127-135, characterized in that the 
step f of multiplying each numeric value of the discrete spectium matrix by a value 
conxsponding to the reciprocal of the Fourier transform of a ftinclion representative 
of the steps of exciting and detecting and providing a system corrected matrix 

5 wherein the discrete va lues of tlie Fourier transform of the system function matrix are 
given by: 

137. The method according to any one of claims 127-136, characterized in that the 
10 step h of correcting each element in the transformed matrix for any asymmetries in 

dimensions of the sample volume wherein the discrete values of are given by: 

138. The method according lo claim lOX further comprising the steps of using the 
15 matrix inversion algoritJun or the reiterative algoritlun in combination with tlie 

Fourier Transform Algorithm. 

139. The method according to claim 103, further comprising the steps of using the 
matrix inversion algorithm and the reiterative algoritlun in combination witli tlie 

20 Fourier Transform Algoritlun. 

1 40. The method according to claim 103, further comprising using the Fourier 
Transform Algoritlun as input for the first reiteration of a reiterative algoridim. 

25 141- The method according to claim 1 40, further comprising the steps of: 

using the Fourier Transform Algorithm to calculate a series of fn^st bulk 
magnetization images for each matrix of components associated by phase; 

detecting at least one peadcs or edge of eacli image by using a cui-off level for the 
bulk magnetization assigned to a voxel or the value of the derivative oftlic bulk 
30 magnetization as a fimction of voxel position^ respectively; 
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setting the bulk itiagnerizatjon of the voxel or voxels corresponding to the peak or 
edge equal to the peak bulk ma^etizalion and assigning the remaining voxels a zero 
value; 

calcuiating a signal at each detector using a system function; 
5 comparing each calculated signal to the actual signal and determining a 

difference or error; 

projecting the error over the voxels using the Fourier Transform Algorithm, 
Matrix Inversion Algorithm, or Reiterative Algorithm; 

adding back projected errors to a preceding corrected series of bulk 
10 magnetization maps; and 

performing the steps reileratively until it stabilizes. 

142. The method according to claim 96, characterized in that a 7*, image is produced 
by a applying at least one pulse sequence that inverts the magnetization and records 

15 the relaxation. 

1 43. The method according to cJaim 1 42, further comprising producing a T^ image by 
the steps of: 

switching on radio fircqucncy receivers to follow a decay following a nuclear 
20 excitation comprising a H, pulse; 

examining (/), a time dependent bulk magnetization in the direction of a 
primary field H^,, at a time after an inverting pulse by applying another Hj pulse 
equivalent to a rotation by 90*" after wailing the time following the initial 
inversion, 

25 characterized in that the 90** pulse puts the z magnetization A/, (/^) into the 
transverse plane for observation, and 

changing the waiting time to allow for observation of MXO at different times 
during relaxation. 

30 144. Tlie mctliod according to ciajjn 96, characterized hi that a 7\ image is produced 
by a applying at least one pulse sequence that Hips a magnetization vector into a 
transverse plajie and records a transverse relaxation by producing at least one spin- 
echo. 
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145. The method according to claim 144, further comprising producing a image by 
the steps of 

flipping the magnetization initially ah'gncd along the z„ -axis into the transverse 
5 plane by applying a H , pulse; 

recovering the relaxation from a FID by applying a 180° pulse at time 

along the yR-axis where the relaxation occurs in time 

recording the 7^ relaxation after an additional time characterized in that the 

total lime elapsed af^er the 90*^ Hj pulse is It^ = which is the spin-echo time and 
10 the dephasing due to the static inhomogeneity of II^ is canceled out, whereby the 

relaxation is due to processes that create relaxation; 

applying a plurality of 1 80^ pulses to form a plurality of spin echoes; and 
determining the relaxation from tlie envelop of the maximum amplitude of the 

spin echoes. 

15 

146. The method according to claim 145» further comprising a pulse sequence to give 
the data for a image known as the Carr-Purcell-Meiboom-Giil sequence 
comprises applying a 90** pulse along the x,^-axis followed by a series of 1 80*^ pulses 
along the y -axis at times + Int^ where n is an integer including zero. 

20 

147. A method of providing a multidimensional nuclear magnetic resonance image of 
an animate or inanimate object comprising the steps of: 

placing a volume of an object having nuclei in a static magnetic field H^^ , the 
volume comprising a plurality of voxels, and then subjected tliem to an additional 
25 rotating radio frequency field H, , ^vhich is synchronous with (heir precession, such 
that the bulk magnetization M of each voxel of the object processes about H, and 
rotates away from the primary field Ho by an angle ^„ ^ ; 

measuring of the intensity of tJie radio frequency signal over time and space 
following r, and/or encoding pulses wherein the magnitude of M is a maximum 
30 initially and decays with time by emission of the same mullipolarity radiation that it 
absorbed and by transfer of energy to tlie surrounding lattice and the iiitensity of the 
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radiation is a function of M and the coordinate position relative to the radio 
frequency emitting voxel; 

Fourier transforming the signal as a function of lime at a given detector position 
to give the components each having an amplitude and a unique phase; 
5 detennining each set of components of the nuclear magnetic resonance signal 

over the sample volume due to a given voxel from the phase data and the detector 
positions; 

determining the location of each voxel through the spatial variations of (he 
intensity of the transverse nuclear magnetic resonance field of the set of components 
1 0 associated by phase; and 

supciimposing the results of the determination of the location of M of each voxel 
and the plotting and displaying the resulting image, 

148. The method according to claim 96, characterized in that tlie plurality of Fourier 
1 5 components each having the same frequency, an intensity and a phase aiigle each 

have a unique phase. 

149. The method according to claim 148, characterized in that the unique phase of 
each component is provided by a combination of 

20 the angle 0 suspended between each spatial location at which tlie secondary field 

is detected and a radial vector from each voxel to each spatial location; and 

the angle ^ due to a separation distance r betv/een each voxel and each spatial 
location given by a wavenumbcr of a detected Geld k tunes r . 

25 1 50. The method according to claim 87, characterized in tliat the radio frequency field 
is sampled synchronously so that tlie phase at any given detector may be related to 
that at any other detector. 

151 . The method according to claim 87, characterized in tliat the radio frequency field 
30 is sampled at known times so that the phase at any given detector may be related (o 
that at any other detector. 
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1 52. The method according to claim 96, characterized in that the processor which 
determines a plurality of Fourier components each having the same frequency, an 
intensity and a phase angle uses a jSrst component having a phase angle and 
calculates (he phase angle as a fimction of spatial position of the first spatial position 

5 relative to any other spatial location and identifies the component at each spatial 
location having the calculated phase angle. 

153. The method according to claim 152, characterized in that the voxels that are on 
Ihe same axis with a given detector or a plurality of detectors that align vertically 

1 0 with one or more voxels is determined. 

1 54. The method according to claim 153, characterized in that the phase angle is a 
fimction of the distance r of tiie detector from the voxel. 

15 155, The method according to claim 154, oharactcrized in that the components may 
be identified by the presence of at least one component of a phase given by ^ + AAr 
where ^, is the phase angle of the first component, is the distance of the second 
detector relative to the first of the detector, and k is the wavenumber of the nuclear 
magnetic resonance radiation emanating from the object wherein this detector(s) is 

20 defined as having a zero angle 0 , 

156. The method according to claim 155, characterized in that a phase angle ^ofsL 
second component at a second detector aligned on the detection axis corresponding 
to a first component at a first detector having a phase angle is given by 

where Ar is the distance of the second detector relative to the first of the detector 
and k is the wavenumber of the nuclear magnetic resonance radiation emanating 
from the object. 

157, The method according to claim 156, characterized in Uiat, from the components 
identified as coaxiai with a detector, the pfiases of the components from the voxel arc 
detemained at the other detectors using the relative position of other detectors of the 
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sample volume relative to the first detector and ^^ of the component at the first 
detector. 



1 58. The method according to claim 157, characterized in that the phase ^ due to a 
5 separation distance between a voxel and a detector of r is given by the wavenimiber 
of the radio frequency field k times 

where the wavenumber k is given by 

-f 

10 the phase angel ^ is given by 

2n 



where Ax, Ay, and Az is the separation distance between a voxel and a detector in 
the X, y, and z directions, respectively, and 

in the case that a two dimensional slice is imaged in the Cransvcrse plane, the 
15 phase angel ^ is given by 



159. The method according to claim 1 58, characterized in that all voxels which have 
(lie same radial distance from a given detector wiU have the same phase 0. 

•20 

160. The method according to claim 159, characterized in that at each case that two 
detectors have the same phase ^, since the radial distances r of two voxels to a 
given detector are equivalent, a unique angle 0 suspended between the direction of 
the detector and the vector along r exists. 

25 

161 . The method according to claim 160, characterized in that the radio frequency 
dipoles of all voxels are time synchronous. 

162. The method according to claim 96, characteri;!:ed in that the nuclear magnetic 
30 resonance pulsing and detection may be synchronized. 
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163, The inetliod according to claim 96, characterized in (hat at least one standard 
dipole may be used as reference phase to set the phases of the radio frequency 
dipolcs of the voxels. 

5 1 64. The method according to claim 163, characterized in that all voxels that are not 
collincar with the axis of a detector posses an angle 0 suspended between the 
direction of the detector and the radial vector from the dipole to die detector. 

1 65. The method according to claim 1 64, characterized in that since the radio 

10 frequency dipolcs are time synchronous, at each point in time, the angle corresponds 
to a phase angle 9 of the radio frequency signal from each voxel at each detector. 

166. The method according to claim 165, characterized in that in the case tliat a two 
dimensional slice is imaged in the transverse plane, the phase angle 0 is given by 

1^ , „, Ax . -,| Ax . ^ kAx 
15 ^ =sm — = sm , ■■ ■■ =sm 



f ' ' — i>Ul t ■.-■■i;- 



167. The method according to claim 166, characterized in that the total unique phase 
angle of the radio frequency signal from each voxel at each detector is given as 
the sum of 0 and 0 given by 
20 0r=tf^ + 0. 
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